Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,983,165 A→C 13.6% T70P (ACC→CCC)  aslB → putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,983,1650AC13.6% 187.6 / 2.3 74T70P (ACC→CCC) aslBputative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (49/13);  new base C (2/10);  total (51/23)
Fisher's exact test for biased strand distribution p-value = 6.97e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACC  >  NC_000913/3982932‑3983400
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
ggacaggcTTCGGGCAAGGGGCGACCATTCTTCAACCGGTTCCAAGGCGGGCTTTTAATTTGAGGGGGAACCCGAGCGGTTCCGATTGTAATCTGAACTGTGACTCCTTGTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTgg                                                                                                                                                                                                                                         <  2:14820/232‑1 (MQ=255)
    tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGGGCGGTTCCGATTGTAATCTGAACTGTGACTACGGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTAACGCATATGGACGATGACACGCTGGAAGCGTAGTACCGTCACTATCTCGCTGCCGGCCAACCGAAAAAAGAAGTGGCTGTTAcacggcgcg                                                                                                                                                                                                                                     >  2:426662/1‑231 (MQ=255)
    tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGCTGACACGCTTGAACCGAATGTCCGTCACTATATCGCTGCCAGCGACCCGCAAAACGAAGTGGCTTTTAACTTGCAgg                                                                                                                                                                                                                                     >  1:517985/1‑239 (MQ=255)
    tGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTCACGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACCCTGGAAACGTATGTCCGTCCCTATATCGCGGCCAGCGAACCGCCAAACGCAGAGGCTTTAACCTggacgg                                                                                                                                                                                                                                     >  2:989114/1‑235 (MQ=255)
      aCTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGCCCGCAACTATATCGATGCCAGCGAACCGCAAAACGCAGAGGCTTTTACCTGCCAgggg                                                                                                                                                                                                                                   >  2:942093/1‑238 (MQ=255)
          cAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCACCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGgcggcg                                                                                                                                                                                                                               >  1:965734/1‑239 (MQ=255)
                 ggAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCaa                                                                                                                                                                                                                         >  1:1081255/1‑238 (MQ=255)
                  gAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGactgcc                                                                                                                                                                                                                         >  1:14820/1‑232 (MQ=255)
                  gAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAAc                                                                                                                                                                                                                        >  1:281101/1‑238 (MQ=255)
                   ggCGGCCATGCTTCAACAGGTTTCATAGCGTTCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGGGACTACTGTTTTTATCTGGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCaa                                                                                                                                                                                                                         <  2:217510/235‑1 (MQ=255)
                   aGCGACCATGCTTCAACAGGGTCCAACGCCGTCTTTTCATGTGGTGGCGAAACCGAGCGGTTCCGATTGTAATCCGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACg                                                                                                                                                                                                                       <  2:281101/238‑1 (MQ=255)
                       aCCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGGTGGCGAAACCGCGCGGTTCCGATTGTAATCTGAACTGTGGCTACTGTTTTTATCTCCAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGTACGCTGACACGACGGAACCGTATGCCCGTCACTATATTCCTGCCACCGAAACGCAAAACGAAGAGGCTTTTACCTGGCAGAGCGGCGCACAAAcacaac                                                                                                                                                                                                                   >  2:280136/1‑233 (MQ=255)
                       aCCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCAAACCAACGCAAcc                                                                                                                                                                                                                  >  2:1014679/1‑238 (MQ=255)
                             ctgcAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGCAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGCAAGGGCGGCGAACCAACGcaccccgg                                                                                                                                                                                                               >  1:135972/1‑229 (MQ=255)
                             ctgcAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCgg                                                                                                                                                                                                               >  1:695350/1‑236 (MQ=255)
                                  aCAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTCTATCCCTGCCAGCGAACCGCAAAACGAATGGGCTTTTACCTGGCAGGGGGGCGAAACAAACCAACACGGGCTggggg                                                                                                                                                                                                       >  2:176840/1‑236 (MQ=255)
                                        tCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCCCTATATCGCTGCCACCGAACCGCAAAACGACGTGGCTTTTACCTGGCATGGCGGCGAACCAACCCCACCCGGGCTGGCGTTTTAcc                                                                                                                                                                                                 >  2:93504/1‑239 (MQ=255)
                                        tCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCCCTGCCAGCGAACCGCCAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAAAGCTACCCGGGCTGGCGTTTTAc                                                                                                                                                                                                  >  1:26018/1‑238 (MQ=255)
                                            aCGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGCCCGTCACTATATCGCTGCCACCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGGGGACCAACGCAACACGGGGTGGCGTTTTAccgccg                                                                                                                                                                                             >  2:553616/1‑239 (MQ=255)
                                                 tGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGACGTGGCTTTTACCTGGCAGGGCGGCGATCCCACGCTACTCGGGCTGGCGTTTTCCCgccgtgccg                                                                                                                                                                                        >  1:608601/1‑239 (MQ=255)
                                                 tGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgtgccg                                                                                                                                                                                        >  1:39800/1‑239 (MQ=255)
                                                         atggaagggCGAACCCGGGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTa                                                                                                                                                                               <  1:913948/232‑1 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGATAAGCCAGTCACGCATATGGACGATGACACGCTGGAACCGTATGTCCGTCACCATATCGCTGCCAGCGAACCGCAAACCGAAGGGGCTTTTCCCTGGCAGGGCGGCGAACCAACCCTACTCGGGCTTGCGTTTTACCCCAGGGCCGGAgg‑gcccccggc                                                                                                                                                                          >  2:473217/1‑232 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTCTGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGACTTTTACCTTGCAGGGCGGCGAACCAACGCTACTCGTGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGc                                                                                                                                                                          >  1:293654/1‑239 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGAAACACGAATGGGCTTTTACCTGGCAGGGCGGCGCACCAACGCAACTCGGGCTGGCGTTTTACCGCCGGGCCGGAGC‑GCAACAGGc                                                                                                                                                                          >  2:83034/1‑239 (MQ=255)
                                                               tgggcaaaCCGGGGGGTTCCAATTGTATCCTGAATTGTACCTTTGTTTTTTATCGCAAAAAACAACCCCTTACCCGCGAAAACCCAGTACCTAATAGGGCCGATGACCCGCTGGAGGCTGATGTCCGTCACTATACCGCTGCCAGCGAACCGAAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGTTCCCGGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCg                                                                                                                                                                         <  2:164626/234‑1 (MQ=255)
                                                                ggCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTCGCAGGCCGGCGAACCAACGCTACCCGGGCTGCGGTTTTACCGCCGTGCCGAAGG‑CCCACAGGGGa                                                                                                                                                                        >  2:902036/1‑239 (MQ=255)
                                                                          aGCGGATCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTAACGCATATGGACGATGACACGCTGGAATCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGGACCAACGCCACTCGGGCTGGGGTTTTACCGCCGTGCCGGCGG‑TCAACAGGCGAAAAATGGGGc                                                                                                                                                              >  2:907911/1‑239 (MQ=255)
                                                                           gCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCATCGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTCCAGGCGAAATATggg                                                                                                                                                                >  1:158885/1‑235 (MQ=255)
                                                                                 ctgATTTTAATCTGCAGTGTCACTACTTTTTTTATCACAAAACACAACCCCTTACCCGCGAAAACCCAGTCACGCATGGAGATGATGACCGGCTGGAAGCTTGTGTCCTGCACTATATCGCGGCAAGCAACCCGCAAACGGAAGTGGCTTTTCCTGGGCGGGGCGGCGAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAgg                                                                                                                                                       <  2:1052054/237‑1 (MQ=255)
                                                                                    gaccaaaTCTGATCGGGTACTTCTTTTTTTATCTCAAAAACCAAGCCCTTCACCACGAAAAGCCAGTCACCAATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTCCTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGa                                                                                                                                                      <  2:1051445/232‑1 (MQ=255)
                                                                                     ttGTAATCTGAACTGGGACTACTGTTTTTATCCCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                                    <  1:90703/238‑1 (MQ=255)
                                                                                     ttGAAATCTGAACTGTGACTCCTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACCCGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                                    <  1:176840/238‑1 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGGATGTCCGTCACTATCTCGCTGCCAGCGAACCTCAAAACGAATTGGCTTTTCCCTGGACGGGCGGCGAACCAACACTACCCGGGATGGCGTGTTACCGCCGTGCCG‑AGCAGCTACAGACCAAACTTGGGGCTGGCAGTAAggaaaa                                                                                                                                               >  2:361304/1‑234 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGG‑CCTACAGGCGACAGATGGGGCTGGCAGGAAGGTAAg                                                                                                                                               >  1:756656/1‑239 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTCCAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                                               >  1:134723/1‑239 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                                               >  1:11000/1‑239 (MQ=255)
                                                                                                    tgACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTAAAGGCGAAATATGGGGCTGGCAGGAAGATAAGTAACAgggc                                                                                                                                      >  1:1051445/1‑234 (MQ=255)
                                                                                                      actactGTTTTTATCTCGAAAAACAATCCCTATACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGCAATATGGTGCTGGCAGGAAGATAAGTAACAACTTCCAg                                                                                                                                  >  1:321045/1‑239 (MQ=255)
                                                                                                        tcctGTTTTAATCTCGAAAAACAATCCTTTTCCCGCGAAAAGCCAGTCACGCATAGGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGACCCGAAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGTTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAAc                                                                                                                                           <  1:460422/226‑1 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGAAAGATAAGGAACAGCTTCCaacct                                                                                                                               >  1:141029/1‑235 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAAAATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACt                                                                                                                               >  1:54570/1‑239 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCCCCGTGCCGTAGC‑GCTAAAGGGCAAATAGGGGGCTGGGGGGGAGAAAAAGAACAGGTTCCAgtaa                                                                                                                               >  2:1016045/1‑236 (MQ=255)
                                                                                                              ttttATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGC‑CCTCCAGGCGAAATTGGGGGCGGGCAGCAAGATAGATAACAGCTTcccaacaaaag                                                                                                                           >  2:936633/1‑229 (MQ=255)
                                                                                                                ttacccagAAAACCAATCCCTTTACCGCGAAAAGCCAGTACCGCATTGGGACGATGACACGCTGGAAGCGTATTTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCg                                                                                                                        <  1:1014679/232‑1 (MQ=255)
                                                                                                                    ctcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGAAAAAAGTTCCAAACAAAAGGCGggcg                                                                                                                    >  2:323812/1‑238 (MQ=255)
                                                                                                                     tcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGGGCTGGCAGGGAGAAAAGAAAAAGAGTCCAGACTAAAggggggggg                                                                                                                   >  2:1098846/1‑232 (MQ=255)
                                                                                                                         aaaaCACTCCCTTTCCCGAAAAAACCCAGTCCCGCATATGGAGGATGACCCGCGGAAAGGTTTTGTCCGTCATTATATCGCTGCCAGCAAACCGCAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGCTACGCGGGGTGGCGTTTGACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTgg                                                                                                               <  2:566285/239‑1 (MQ=255)
                                                                                                                                     ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATg                                                                                                     >  1:566285/1‑237 (MQ=255)
                                                                                                                                      tACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGGGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATTTGGGTCGGGCAGGAAGATAAGTAAAAGCTACAAGCAAAACGGCGAGATGCCGGAGGAGGAAggggg                                                                                                  >  2:493439/1‑237 (MQ=255)
                                                                                                                                           cgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGt                                                                                             >  1:637424/1‑239 (MQ=255)
                                                                                                                                            gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGTTGACGAATGGGGCGCGtt                                                                                            >  1:361870/1‑239 (MQ=255)
                                                                                                                                                     gTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                                   >  1:288179/1‑239 (MQ=255)
                                                                                                                                                      tCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTCCTGCTGGATGACGAATGGTGCGCGTTTTTCCCggcc                                                                                  >  1:1069415/1‑237 (MQ=255)
                                                                                                                                                               tGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTTACAGCTTCCCGACTAACGGCGGTCTGCTGGATGACGAAAGGGGCGCGTTTTTCCCgggagctcctttt                                                                         >  2:469494/1‑228 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATATGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGGGCGCGTTTCCCGCGGAGCATCATTTTCttgttg                                                                  >  1:424073/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgttg                                                                  >  1:305942/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATACTTTTCttgttg                                                                  >  1:938586/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGGAACAGCTTCCCGACTAAAGGCGTGGTGCTGGATGACGAATGGTGCGCGGTTCCCGCGGGGCATCATTTTCtttttg                                                                  >  2:1061349/1‑239 (MQ=255)
                                                                                                                                                                                        tgtcctgctttattTCGCGCCCCGCGACCCTAAAACCGGAGTGGCTTTTACTGAGCAGGCGGGCGAACCACGGCTACCGGGGTGGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                                                <  2:468053/229‑1 (MQ=255)
                                                                                                                                                                                        tGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGGGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTTACGCTGGATgtc                                                 >  1:830146/1‑236 (MQ=255)
                                                                                                                                                                                        tGTCCGTCACCATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGCAGATACGGAACAGCTTCCAGACTAACGGCGGGCTGCTGACTGGCGAAGGGGGCGCGTTTCTCGCGGGGCATCATTTTTTGGTTGGGTTATcagtggggggg                                                 >  2:972932/1‑228 (MQ=255)
                                                                                                                                                                                              tCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTCACGGCGTGCTGCTGGTTGACGAATGGGGGGCGTTTCTCGCGGGGCATACTTTTTTTGTTTGTTTTTCCATGGATGGCCCGcccg                                          >  2:372959/1‑237 (MQ=255)
                                                                                                                                                                                                    tatCGCCGCCAGCAAACCGAAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTACTGGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:323812/239‑1 (MQ=255)
                                                                                                                                                                                                    tagcggtcccACCGCACCGAAAACCGTGGTGGCTTTTCCCTGGCGGGGCGGCGACCCAACGTTACTCGGGCGGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:133531/231‑1 (MQ=255)
                                                                                                                                                                                                    taccGCTGCCGGCAAACGGAAAAACGAATTGGTTTTTCCCGGGCGGGGCGGCGAACCAACGCTACTCGGGCTGGCTTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCGCGCCTGAGATcc                                    <  1:936633/236‑1 (MQ=255)
                                                                                                                                                                                                    aaccGCTCCCAGCGAACCGCAAACCGAGGGTGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:101926/236‑1 (MQ=255)
                                                                                                                                                                                                        gCTGCAAGCGAACCGCAAAGCGAAGGGGCTTGTCCCTGGCAGGGCGGCGAACCAATGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa                                <  1:907911/239‑1 (MQ=255)
                                                                                                                                                                                                        gCTCCCGGCAAACCGCAAACCGAAGTGCTTTTTCCCTGGAAGGGCGGCAAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa                                <  1:942093/239‑1 (MQ=255)
                                                                                                                                                                                                         cTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTaa                                                                                                                                            >  1:188662/1‑130 (MQ=255)
                                                                                                                                                                                                                   aCCGCAAAACGAAGTGGCTTTTACCTGGAAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc                      <  1:181482/238‑1 (MQ=255)
                                                                                                                                                                                                                   aCCGAAAACCGAGGGGGTTTTTCCTTGGCAGGGGGGCGAACCAACGCTACCCGGGCGGGCTTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCATACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc                      <  1:93504/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCTTTTCCCTGGCAGGGCGGCAAACCAACGCTTCTGGGCCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTgg        <  1:469494/237‑1 (MQ=255)
                                                                                                                                                                                                                                      tttCCCTGGGCGGGCGGCTAACCACGCCTACGCGGGGTGGCTTTTTACCGCCGTGCCGTAGG‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACCGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTTCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc  <  2:39800/239‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
AGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGC‑GCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACC  >  NC_000913/3982932‑3983400

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