Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,022,380 A→T 6.1% K55* (AAA→TAA)  tatB → TatABCE protein translocation system subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,022,3800AT 6.1% 465.9 / 16.3 132K55* (AAA→TAA) tatBTatABCE protein translocation system subunit
Reads supporting (aligned to +/- strand):  ref base A (61/63);  new base T (5/3);  total (66/66)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 139 total aligned reads displayed.

TTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAA  >  NC_000913/4022232‑4022525
                                                                                                                                                    |                                                                                                                                                 
ttttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTaa                                                                                                                                                   <  1:335798/149‑1 (MQ=255)
  ttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTaaaa                                                                                                                                                 <  1:748943/149‑1 (MQ=255)
  ttAGCGAACTGCTAGTGGTCTTCATCATCGGCCTCGTCGCTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTaaaa                                                                                                                                                 <  4:779544/149‑1 (MQ=255)
    aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACt                                                                                                                                               >  1:1060150/1‑149 (MQ=255)
    aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACt                                                                                                                                               >  2:447416/1‑149 (MQ=255)
    aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACt                                                                                                                                               >  2:307665/1‑149 (MQ=255)
     gCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTc                                                                                                                                              >  3:187777/1‑149 (MQ=255)
      cGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTcc                                                                                                                                             <  3:714349/149‑1 (MQ=255)
       gAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCa                                                                                                                                            <  2:155170/149‑1 (MQ=255)
        aagtGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGGCGCAACGACTGCTTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAg                                                                                                                                           >  2:988482/4‑149 (MQ=255)
        aaCTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAg                                                                                                                                           <  4:940366/149‑1 (MQ=255)
        aaCTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAg                                                                                                                                           <  2:486237/149‑1 (MQ=255)
             cTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGtt                                                                                                                                      <  1:1128396/149‑1 (MQ=255)
                ttGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGTGTTTCa                                                                                                                                   <  2:814235/149‑1 (MQ=255)
                 tGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAg                                                                                                                                  <  3:1062847/149‑1 (MQ=255)
                   gtgtTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGa                                                                                                                                <  1:558513/149‑1 (MQ=255)
                   gtgtTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGa                                                                                                                                <  2:547272/149‑1 (MQ=255)
                     gtTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACa                                                                                                                              <  4:644221/149‑1 (MQ=255)
                            atcGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGaa                                                                                                                       >  3:448717/1‑149 (MQ=255)
                                gCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAg                                                                                                                   >  1:930346/1‑149 (MQ=255)
                                gCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAg                                                                                                                   <  2:253728/149‑1 (MQ=255)
                                gCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGCTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGGAAAAg                                                                                                                   >  1:775312/1‑149 (MQ=255)
                                 ccTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAgg                                                                                                                  <  3:270013/149‑1 (MQ=255)
                                     gtcgtTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGa                                                                                                              >  2:699018/1‑149 (MQ=255)
                                        gtTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTc                                                                                                                           >  3:836026/1‑133 (MQ=255)
                                        gtTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTc                                                                                                                           <  4:836026/133‑1 (MQ=255)
                                        gtTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAACAGGTTGaaaa                                                                                                           <  4:942171/149‑1 (MQ=255)
                                         ttCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGa                                                                                                              >  1:654526/1‑145 (MQ=255)
                                         ttCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGa                                                                                                              <  2:654526/145‑1 (MQ=255)
                                         ttCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAg                                                                                                          >  4:737497/1‑149 (MQ=255)
                                           cTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTTTAACTCCTgg                                                                                                                                          >  4:814132/1‑115 (MQ=255)
                                           cTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGc                                                                                                        <  2:1043774/149‑1 (MQ=255)
                                            tGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCg                                                                                                       >  3:883661/1‑149 (MQ=255)
                                             ggggCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGa                                                                                                      <  2:20584/149‑1 (MQ=255)
                                              gggCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAg                                                                                                     >  4:1025714/1‑149 (MQ=255)
                                                  cGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTc                                                                                                 <  3:952127/149‑1 (MQ=255)
                                                   gCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCa                                                                                                >  4:1108916/1‑149 (MQ=255)
                                                    cAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAgg                                                                                                                                 >  2:649135/1‑115 (MQ=255)
                                                    cAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAgg                                                                                                                                 <  1:649135/115‑1 (MQ=255)
                                                        gACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCGTCACTAAc                                                                                           <  3:1055298/149‑1 (MQ=255)
                                                        gACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCt                                                                                                  <  1:808377/142‑1 (MQ=255)
                                                          cTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCt                                                                                         <  1:216722/149‑1 (MQ=255)
                                                           tGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTaaaa                                                                                                                                                 <  1:90300/92‑1 (MQ=255)
                                                           tGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTaaaa                                                                                                                                                 >  2:90300/1‑92 (MQ=255)
                                                           tGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGaaaaa                                                                                                                    <  3:797958/121‑1 (MQ=255)
                                                            gCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAGAGGCGAGCCTCACTAACCTGa                                                                                       >  4:1078407/1‑149 (MQ=255)
                                                            gCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGa                                                                                       <  3:1019220/149‑1 (MQ=255)
                                                            gCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGa                                                                                       <  2:1060150/149‑1 (MQ=255)
                                                              cTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGGAAAGGCGa                                                                                                      >  1:873078/1‑132 (MQ=255)
                                                              cTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACg                                                                                     >  4:154312/1‑149 (MQ=255)
                                                                gtgGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGcc                                                                                   >  1:234696/1‑149 (MQ=255)
                                                                gtgGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGcc                                                                                   >  4:1050487/1‑149 (MQ=255)
                                                                        aaaaaCGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGa                                                                           >  3:66729/1‑149 (MQ=255)
                                                                         aaaaCGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGaa                                                                          >  4:108308/1‑149 (MQ=255)
                                                                          aaaCGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCTGGACAGCCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGaaa                                                                         <  1:985186/149‑1 (MQ=255)
                                                                          aaaCGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGaaa                                                                         <  1:415351/149‑1 (MQ=255)
                                                                                tAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCg                                                                   >  1:612988/1‑149 (MQ=255)
                                                                                 aGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGGCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCga                                                                  <  4:448717/149‑1 (MQ=255)
                                                                                  gCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgat                                                                 >  2:708825/1‑149 (MQ=255)
                                                                                  gCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgat                                                                 <  2:234696/149‑1 (MQ=255)
                                                                                     ggCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAg                                                                                                                   <  3:478143/96‑1 (MQ=255)
                                                                                     ggCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAg                                                                                                                   >  4:478143/1‑96 (MQ=255)
                                                                                      gCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatggat                                                             <  4:66729/149‑1 (MQ=255)
                                                                                      gCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgatggat                                                             <  2:930346/149‑1 (MQ=255)
                                                                                        tGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAg                                                                                                                                        >  2:959577/1‑72 (MQ=255)
                                                                                        tGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGa                                                                                                      <  4:83016/106‑1 (MQ=255)
                                                                                        tGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGa                                                                                                      >  3:83016/1‑106 (MQ=255)
                                                                                          gATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGa                                                                                                              >  3:92545/1‑96 (MQ=255)
                                                                                               gcgcgTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTcc                                                                                                                                             <  4:250913/60‑1 (MQ=255)
                                                                                                   gttgcgttCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAgg                                                <  3:745407/149‑1 (MQ=255)
                                                                                                    ttgcgttCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGc                                               >  2:1023533/1‑149 (MQ=255)
                                                                                                        gttCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCgat                                                                 <  2:833070/127‑1 (MQ=255)
                                                                                                            aCTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGaa                                                          >  3:143177/1‑130 (MQ=255)
                                                                                                             cTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTc                                      <  4:187777/149‑1 (MQ=255)
                                                                                                              tGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCg                                     >  1:1045779/1‑149 (MQ=255)
                                                                                                                gCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                   <  1:988482/149‑1 (MQ=255)
                                                                                                                  gACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGa                                 >  2:614261/1‑149 (MQ=255)
                                                                                                                         gTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTc                          <  1:1023533/149‑1 (MQ=255)
                                                                                                                              gaacgaacTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACg                     >  2:286535/1‑149 (MQ=255)
                                                                                                                              gaacgaacTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACg                     >  4:571574/1‑149 (MQ=255)
                                                                                                                               aacgaacTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTc                                                                    <  2:570798/101‑1 (MQ=255)
                                                                                                                               aacgaacTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                   <  2:393396/134‑1 (MQ=255)
                                                                                                                               aacgaacTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGAt                                   >  1:393396/1‑134 (MQ=255)
                                                                                                                                    acTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaa               >  2:764048/1‑149 (MQ=255)
                                                                                                                                     cTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCaaa              <  4:1001440/149‑1 (MQ=255)
                                                                                                                                         tcctggAGTTTTAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGa                                    >  1:528337/5‑123 (MQ=255)
                                                                                                                                         tcctggAGTTTTAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGAt          >  1:169448/5‑149 (MQ=255)
                                                                                                                                         cccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGcc                                                                                                   >  4:1128217/1‑60 (MQ=255)
                                                                                                                                          cctggAGTTTTAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc         <  2:169448/146‑1 (MQ=255)
                                                                                                                                          cctggAGTTTTAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc         >  3:386525/4‑149 (MQ=255)
                                                                                                                                          ccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc         <  1:18006/149‑1 (MQ=255)
                                                                                                                                          ccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc         <  3:1050487/149‑1 (MQ=255)
                                                                                                                                          ccAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATc         <  3:1108916/149‑1 (MQ=255)
                                                                                                                                            tggAGTTTTAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt       >  3:1049683/2‑149 (MQ=255)
                                                                                                                                            aGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCt       >  3:300839/1‑149 (MQ=255)
                                                                                                                                                 ttAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa  >  1:679346/1‑149 (MQ=255)
                                                                                                                                                 ttAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa  <  1:699018/149‑1 (MQ=255)
                                                                                                                                                 ttAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa  >  3:247352/1‑149 (MQ=255)
                                                                                                                                                 ttAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGCTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGaaaa  >  4:1084639/1‑149 (MQ=255)
                                                                                                                                                   aaaaCTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAg                                                                                                          >  3:567218/1‑43 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCGTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGCGTTGCGTTCACTGGCGACAACGGTGCAGAACGAACTGACCCAGGAGTTAAAACTCCAGGAGTTTCAGGACAGTCTGAAAAAGGTTGAAAAGGCGAGCCTCACTAACCTGACGCCCGAACTGAAAGCGTCGATGGATGAACTACGCCAGGCCGCGGAGTCGATGAAGCGTTCCTACGTTGCAAACGATCCTGAAAA  >  NC_000913/4022232‑4022525

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