Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 7.9% coding (685/1131 nt) wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,8171.G 7.9% 540.2 / 42.1 127coding (685/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (62/55);  new base G (5/5);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 130 total aligned reads displayed.

CGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970673‑3970962
                                                                                                                                                 |                                                                                                                                                 
cGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑Tgg                                                                                                                                                >  2:763627/1‑148 (MQ=255)
         gTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCg                                                                                                                                       <  1:763627/148‑1 (MQ=255)
                    cGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCAc                                                                                                                            >  2:1074412/1‑148 (MQ=255)
                           cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAgg                                                                                                                        <  4:1009642/145‑1 (MQ=255)
                           cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAgg                                                                                                                        >  3:1009642/1‑145 (MQ=255)
                             gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                           >  1:183943/1‑124 (MQ=255)
                             gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                           <  2:183943/124‑1 (MQ=255)
                              cGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                  <  1:1074412/148‑1 (MQ=255)
                                    gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                            <  3:411070/148‑1 (MQ=255)
                                     aaCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                           <  3:840391/148‑1 (MQ=255)
                                      aCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACAt                                                                                                                                         >  2:375194/1‑117 (MQ=255)
                                         ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑Tg                                                                                                                                                 >  1:963606/1‑106 (MQ=255)
                                          gcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                      >  4:804854/1‑148 (MQ=255)
                                             gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            <  4:317038/107‑1 (MQ=255)
                                                aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCa                                                                                                                             >  2:454142/1‑119 (MQ=255)
                                                aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCa                                                                                                                             <  1:454142/119‑1 (MQ=255)
                                                     cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          <  2:301678/117‑1 (MQ=255)
                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTg                                                                                                 >  2:606509/1‑140 (MQ=255)
                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTg                                                                                                 <  1:606509/140‑1 (MQ=255)
                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                         >  4:1782/1‑148 (MQ=255)
                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            <  2:234661/97‑1 (MQ=255)
                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            >  1:234661/1‑97 (MQ=255)
                                                         tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGcca                                                                                                                                    <  2:636703/103‑1 (MQ=255)
                                                         tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            >  2:482740/1‑95 (MQ=255)
                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  3:403931/1‑112 (MQ=255)
                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  3:803321/1‑112 (MQ=255)
                                                               gACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACg                                                                                 >  1:347439/1‑148 (MQ=255)
                                                                aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacc                                                                                                                                  >  4:56312/1‑98 (MQ=255)
                                                                aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacc                                                                                                                                  <  3:56312/98‑1 (MQ=255)
                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                            <  3:946646/117‑1 (MQ=255)
                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                            >  4:946646/1‑117 (MQ=255)
                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACCCACCAGGACa                                                                                                                     <  3:1036493/108‑1 (MQ=255)
                                                                       tAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                               <  4:350939/126‑1 (MQ=255)
                                                                       tAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                               >  3:350939/1‑126 (MQ=255)
                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACAt                                                                                                                                         >  1:381735/1‑82 (MQ=255)
                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCGCCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAcc                                                                       <  1:739201/148‑1 (MQ=255)
                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            >  3:1038087/1‑79 (MQ=255)
                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            <  4:1038087/79‑1 (MQ=255)
                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAt                                                                                                              <  1:297348/107‑1 (MQ=255)
                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                          >  4:108778/1‑127 (MQ=255)
                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                          <  3:108778/127‑1 (MQ=255)
                                                                              atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  3:980559/1‑92 (MQ=255)
                                                                              atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                      <  4:185349/112‑1 (MQ=255)
                                                                               tgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                 <  3:1782/148‑1 (MQ=255)
                                                                                  ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                         >  4:522219/1‑137 (MQ=255)
                                                                                   tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  3:114782/1‑148 (MQ=255)
                                                                                   tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  4:746244/1‑148 (MQ=255)
                                                                                         gACCGGTCACAGGCGTGAGAGTTTCGGTCGTAGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                        <  1:1059260/147‑1 (MQ=255)
                                                                                         gACCGGTCACAGGCGTGAGAGTTTCGGTCGTAGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                        >  2:1059260/1‑147 (MQ=255)
                                                                                            cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                     <  2:805522/83‑1 (MQ=255)
                                                                                            cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                     >  4:64249/1‑147 (MQ=255)
                                                                                               tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                                >  2:631509/1‑117 (MQ=255)
                                                                                               tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                                <  1:631509/117‑1 (MQ=255)
                                                                                                        tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                               >  1:729713/1‑141 (MQ=255)
                                                                                                        tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                               <  2:729713/141‑1 (MQ=255)
                                                                                                          agagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGt                                      >  1:619826/1‑148 (MQ=255)
                                                                                                           gagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                     >  1:157930/1‑148 (MQ=255)
                                                                                                             gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                         >  2:234245/1‑142 (MQ=255)
                                                                                                             gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                   <  1:319124/148‑1 (MQ=255)
                                                                                                               ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                               <  2:403261/102‑1 (MQ=255)
                                                                                                               ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                               >  1:403261/1‑102 (MQ=255)
                                                                                                                  ggTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                              <  2:347439/148‑1 (MQ=255)
                                                                                                                      gTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcga                          <  3:746244/148‑1 (MQ=255)
                                                                                                                            ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCTCCa                    >  4:647307/1‑148 (MQ=255)
                                                                                                                            ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    <  4:114782/148‑1 (MQ=255)
                                                                                                                            ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAATGTCTATCCgg                                                                                                   <  2:220711/69‑1 (MQ=255)
                                                                                                                            ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAATGTCTATCCgg                                                                                                   >  1:220711/1‑69 (MQ=255)
                                                                                                                              gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAgg                  <  3:235753/148‑1 (MQ=255)
                                                                                                                               aagaaATCTGTCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                 >  1:837971/1‑148 (MQ=255)
                                                                                                                               aagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                        >  3:209656/1‑77 (MQ=255)
                                                                                                                               aagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                        <  4:209656/77‑1 (MQ=255)
                                                                                                                               aagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                 >  4:911012/1‑148 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCa                                                  <  2:1126709/114‑1 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCa                                                  >  1:1126709/1‑114 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                             <  1:297872/135‑1 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                             >  2:297872/1‑135 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                >  4:827532/1‑148 (MQ=255)
                                                                                                                                agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                >  4:173334/1‑148 (MQ=255)
                                                                                                                                 gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             <  3:510192/102‑1 (MQ=255)
                                                                                                                                 gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  4:510192/1‑102 (MQ=255)
                                                                                                                                 gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    >  4:1021048/1‑143 (MQ=255)
                                                                                                                                 gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    <  3:1021048/143‑1 (MQ=255)
                                                                                                                                 gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt               >  4:403755/1‑148 (MQ=255)
                                                                                                                                   tgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                        >  3:592322/3‑89 (MQ=255)
                                                                                                                                    aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtg                                                <  1:925456/112‑1 (MQ=255)
                                                                                                                                    aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtg                                                >  2:925456/1‑112 (MQ=255)
                                                                                                                                    aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt            <  3:911012/148‑1 (MQ=255)
                                                                                                                                     tCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt           >  3:192/1‑148 (MQ=255)
                                                                                                                                      cTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTa          <  4:952005/148‑1 (MQ=255)
                                                                                                                                          cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt      >  1:784876/1‑148 (MQ=255)
                                                                                                                                          cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCACCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATCCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt      <  3:647307/148‑1 (MQ=255)
                                                                                                                                           aCGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                         <  1:316913/112‑1 (MQ=255)
                                                                                                                                           aCGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                         >  2:316913/1‑112 (MQ=255)
                                                                                                                                            cgcgC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                 <  1:359774/87‑1 (MQ=255)
                                                                                                                                            cgcgC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                 >  2:359774/1‑87 (MQ=255)
                                                                                                                                            cgcgC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt           <  1:1025974/141‑1 (MQ=255)
                                                                                                                                            cgcgC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGttt    >  4:90102/1‑148 (MQ=255)
                                                                                                                                              cgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                              <  4:596624/120‑1 (MQ=255)
                                                                                                                                              cgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                              >  3:596624/1‑120 (MQ=255)
                                                                                                                                              cgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGt  <  3:827532/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970673‑3970962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑