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

AGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACC  >  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)
                                                         atggaagggCGAACCCGGGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTa                                                                                                                                                                               <  1:913948/232‑1 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGATAAGCCAGTCACGCATATGGACGATGACACGCTGGAACCGTATGTCCGTCACCATATCGCTGCCAGCGAACCGCAAACCGAAGGGGCTTTTCCCTGGCAGGGCGGCGAACCAACCCTACTCGGGCTTGCGTTTTACCCCAGGGCCGGAgggcccccggc                                                                                                                                                                          >  2:473217/1‑232 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTCTGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGACTTTTACCTTGCAGGGCGGCGAACCAACGCTACTCGTGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGc                                                                                                                                                                          >  1:293654/1‑239 (MQ=255)
                                                              aTGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGAAACACGAATGGGCTTTTACCTGGCAGGGCGGCGCACCAACGCAACTCGGGCTGGCGTTTTACCGCCGGGCCGGAGCGCAACAGGc                                                                                                                                                                          >  2:83034/1‑239 (MQ=255)
                                                               tgggcaaaCCGGGGGGTTCCAATTGTATCCTGAATTGTACCTTTGTTTTTTATCGCAAAAAACAACCCCTTACCCGCGAAAACCCAGTACCTAATAGGGCCGATGACCCGCTGGAGGCTGATGTCCGTCACTATACCGCTGCCAGCGAACCGAAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGTTCCCGGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCg                                                                                                                                                                         <  2:164626/234‑1 (MQ=255)
                                                                ggCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTCGCAGGCCGGCGAACCAACGCTACCCGGGCTGCGGTTTTACCGCCGTGCCGAAGGCCCACAGGGGa                                                                                                                                                                        >  2:902036/1‑239 (MQ=255)
                                                                          aGCGGATCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTAACGCATATGGACGATGACACGCTGGAATCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGGACCAACGCCACTCGGGCTGGGGTTTTACCGCCGTGCCGGCGGTCAACAGGCGAAAAATGGGGc                                                                                                                                                              >  2:907911/1‑239 (MQ=255)
                                                                           gCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCATCGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTCCAGGCGAAATATggg                                                                                                                                                                >  1:158885/1‑235 (MQ=255)
                                                                                 ctgATTTTAATCTGCAGTGTCACTACTTTTTTTATCACAAAACACAACCCCTTACCCGCGAAAACCCAGTCACGCATGGAGATGATGACCGGCTGGAAGCTTGTGTCCTGCACTATATCGCGGCAAGCAACCCGCAAACGGAAGTGGCTTTTCCTGGGCGGGGCGGCGAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAgg                                                                                                                                                       <  2:1052054/237‑1 (MQ=255)
                                                                                    gaccaaaTCTGATCGGGTACTTCTTTTTTTATCTCAAAAACCAAGCCCTTCACCACGAAAAGCCAGTCACCAATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTCCTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGa                                                                                                                                                      <  2:1051445/232‑1 (MQ=255)
                                                                                     ttGTAATCTGAACTGGGACTACTGTTTTTATCCCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                                    <  1:90703/238‑1 (MQ=255)
                                                                                     ttGAAATCTGAACTGTGACTCCTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACCCGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAg                                                                                                                                                    <  1:176840/238‑1 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAACCAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGGATGTCCGTCACTATCTCGCTGCCAGCGAACCTCAAAACGAATTGGCTTTTCCCTGGACGGGCGGCGAACCAACACTACCCGGGATGGCGTGTTACCGCCGTGCCGAGCAGCTACAGACCAAACTTGGGGCTGGCAGTAAggaaaa                                                                                                                                               >  2:361304/1‑234 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGGCCTACAGGCGACAGATGGGGCTGGCAGGAAGGTAAg                                                                                                                                               >  1:756656/1‑239 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTCCAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                                               >  1:134723/1‑239 (MQ=255)
                                                                                         aaTCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                                                               >  1:11000/1‑239 (MQ=255)
                                                                                                    tgACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTAAAGGCGAAATATGGGGCTGGCAGGAAGATAAGTAACAgggc                                                                                                                                      >  1:1051445/1‑234 (MQ=255)
                                                                                                      actactGTTTTTATCTCGAAAAACAATCCCTATACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGCAATATGGTGCTGGCAGGAAGATAAGTAACAACTTCCAg                                                                                                                                  >  1:321045/1‑239 (MQ=255)
                                                                                                        tcctGTTTTAATCTCGAAAAACAATCCTTTTCCCGCGAAAAGCCAGTCACGCATAGGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGACCCGAAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGTTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAAc                                                                                                                                           <  1:460422/226‑1 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGAAAGATAAGGAACAGCTTCCaacct                                                                                                                               >  1:141029/1‑235 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAAAATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACt                                                                                                                               >  1:54570/1‑239 (MQ=255)
                                                                                                         actGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCCCCGTGCCGTAGCGCTAAAGGGCAAATAGGGGGCTGGGGGGGAGAAAAAGAACAGGTTCCAgtaa                                                                                                                               >  2:1016045/1‑236 (MQ=255)
                                                                                                              ttttATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGCCCTCCAGGCGAAATTGGGGGCGGGCAGCAAGATAGATAACAGCTTcccaacaaaag                                                                                                                           >  2:936633/1‑229 (MQ=255)
                                                                                                                ttacccagAAAACCAATCCCTTTACCGCGAAAAGCCAGTACCGCATTGGGACGATGACACGCTGGAAGCGTATTTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCg                                                                                                                        <  1:1014679/232‑1 (MQ=255)
                                                                                                                    ctcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGAAAAAAGTTCCAAACAAAAGGCGggcg                                                                                                                    >  2:323812/1‑238 (MQ=255)
                                                                                                                     tcGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGGGCTGGCAGGGAGAAAAGAAAAAGAGTCCAGACTAAAggggggggg                                                                                                                   >  2:1098846/1‑232 (MQ=255)
                                                                                                                         aaaaCACTCCCTTTCCCGAAAAAACCCAGTCCCGCATATGGAGGATGACCCGCGGAAAGGTTTTGTCCGTCATTATATCGCTGCCAGCAAACCGCAAAACGAAGTGGCTTTTACCGGGCAGGGCGGCGAACCAACGCTACGCGGGGTGGCGTTTGACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTgg                                                                                                               <  2:566285/239‑1 (MQ=255)
                                                                                                                                     ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATg                                                                                                     >  1:566285/1‑237 (MQ=255)
                                                                                                                                      tACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGGGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATTTGGGTCGGGCAGGAAGATAAGTAAAAGCTACAAGCAAAACGGCGAGATGCCGGAGGAGGAAggggg                                                                                                  >  2:493439/1‑237 (MQ=255)
                                                                                                                                           cgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGt                                                                                             >  1:637424/1‑239 (MQ=255)
                                                                                                                                            gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGTTGACGAATGGGGCGCGtt                                                                                            >  1:361870/1‑239 (MQ=255)
                                                                                                                                                     gTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                                   >  1:288179/1‑239 (MQ=255)
                                                                                                                                                      tCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTCCTGCTGGATGACGAATGGTGCGCGTTTTTCCCggcc                                                                                  >  1:1069415/1‑237 (MQ=255)
                                                                                                                                                               tGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTTACAGCTTCCCGACTAACGGCGGTCTGCTGGATGACGAAAGGGGCGCGTTTTTCCCgggagctcctttt                                                                         >  2:469494/1‑228 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATATGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGGGCGCGTTTCCCGCGGAGCATCATTTTCttgttg                                                                  >  1:424073/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgttg                                                                  >  1:305942/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATACTTTTCttgttg                                                                  >  1:938586/1‑239 (MQ=255)
                                                                                                                                                                      tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGGAACAGCTTCCCGACTAAAGGCGTGGTGCTGGATGACGAATGGTGCGCGGTTCCCGCGGGGCATCATTTTCtttttg                                                                  >  2:1061349/1‑239 (MQ=255)
                                                                                                                                                                                        tgtcctgctttattTCGCGCCCCGCGACCCTAAAACCGGAGTGGCTTTTACTGAGCAGGCGGGCGAACCACGGCTACCGGGGTGGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                                                <  2:468053/229‑1 (MQ=255)
                                                                                                                                                                                        tGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGGGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTTACGCTGGATgtc                                                 >  1:830146/1‑236 (MQ=255)
                                                                                                                                                                                        tGTCCGTCACCATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGGGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGCAGATACGGAACAGCTTCCAGACTAACGGCGGGCTGCTGACTGGCGAAGGGGGCGCGTTTCTCGCGGGGCATCATTTTTTGGTTGGGTTATcagtggggggg                                                 >  2:972932/1‑228 (MQ=255)
                                                                                                                                                                                              tCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTCACGGCGTGCTGCTGGTTGACGAATGGGGGGCGTTTCTCGCGGGGCATACTTTTTTTGTTTGTTTTTCCATGGATGGCCCGcccg                                          >  2:372959/1‑237 (MQ=255)
                                                                                                                                                                                                    tatCGCCGCCAGCAAACCGAAAAACGAAGTGGCTTTTCCCTGGCAGGGCGGCGAACCAACGCTACTGGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:323812/239‑1 (MQ=255)
                                                                                                                                                                                                    tagcggtcccACCGCACCGAAAACCGTGGTGGCTTTTCCCTGGCGGGGCGGCGACCCAACGTTACTCGGGCGGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:133531/231‑1 (MQ=255)
                                                                                                                                                                                                    taccGCTGCCGGCAAACGGAAAAACGAATTGGTTTTTCCCGGGCGGGGCGGCGAACCAACGCTACTCGGGCTGGCTTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCGCGCCTGAGATcc                                    <  1:936633/236‑1 (MQ=255)
                                                                                                                                                                                                    aaccGCTCCCAGCGAACCGCAAACCGAGGGTGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc                                    <  1:101926/236‑1 (MQ=255)
                                                                                                                                                                                                        gCTGCAAGCGAACCGCAAAGCGAAGGGGCTTGTCCCTGGCAGGGCGGCGAACCAATGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa                                <  1:907911/239‑1 (MQ=255)
                                                                                                                                                                                                        gCTCCCGGCAAACCGCAAACCGAAGTGCTTTTTCCCTGGAAGGGCGGCAAACCAACGCTACGCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa                                <  1:942093/239‑1 (MQ=255)
                                                                                                                                                                                                         cTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTaa                                                                                                                                            >  1:188662/1‑130 (MQ=255)
                                                                                                                                                                                                                   aCCGCAAAACGAAGTGGCTTTTACCTGGAAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc                      <  1:181482/238‑1 (MQ=255)
                                                                                                                                                                                                                   aCCGAAAACCGAGGGGGTTTTTCCTTGGCAGGGGGGCGAACCAACGCTACCCGGGCGGGCTTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCATACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc                      <  1:93504/238‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCTTTTCCCTGGCAGGGCGGCAAACCAACGCTTCTGGGCCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTgg        <  1:469494/237‑1 (MQ=255)
                                                                                                                                                                                                                                      tttCCCTGGGCGGGCGGCTAACCACGCCTACGCGGGGTGGCTTTTTACCGCCGTGCCGTAGGGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACCGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTTCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc  <  2:39800/239‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
AGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACC  >  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: