Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 635,136 C→T 11.5% E81K (GAA→AAA)  ybdM ← Spo0J family protein, ParB‑like nuclease domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913635,1360CT11.5% 224.7 / 3.2 78E81K (GAA→AAA) ybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (26/42);  new base T (7/3);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 8.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCC  >  NC_000913/634915‑635374
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ggTGATTTGAAGGCGTCCACGAGCGCGATTATGGCGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATATCCTTTTAAGCGTAATTTTAATGATGATGATCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTgcagc                                                                                                                                                                                                                                    <  1:600512/238‑1 (MQ=255)
   gATTTGGTGGCCTCCCCGCACGCGATTATGACGAATAGTCCGTGTGATGCTTTGATTACGCGTACCTTCCAGACAGGTCCCCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGt                                                                                                                                                                                                                                 <  1:474909/238‑1 (MQ=255)
        tgtgggttCCCCGAGCGCGTTTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGttttttt                                                                                                                                                                                                                           <  2:317734/232‑1 (MQ=255)
              gTCCACGAGCGCGATTATGACGAAAAGGCGGTGCGATTCGTTTAATACGCGTACCTTCCAGACAGGTCCCCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCggg                                                                                                                                                                                                                      <  1:442921/238‑1 (MQ=255)
              ctccccggggcggtttATGACGGATTGGCGGGGGGATGGGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCggg                                                                                                                                                                                                                      <  2:98502/224‑1 (MQ=255)
              ctCCACGGGCGCTATTATGACGAATAGTCGGGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCggg                                                                                                                                                                                                                      <  1:347364/237‑1 (MQ=255)
                      gcgcGATTAAGGCGAATAGTCGCGGCGAGGCGTTGATTACGCGTACCTTCCAGGAAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATTGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                              <  1:412843/238‑1 (MQ=255)
                       cgggATTAAGAACAATAGTCGCGGGTATGGGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGGGTAATTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                              <  1:1061632/234‑1 (MQ=255)
                        gcgATTATGACGAATAGGCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACg                                                                                                                                                                                                            <  1:368577/238‑1 (MQ=255)
                         cgATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGGAAACACATCACGAGTTCCATGTGACTTTTTTTCTCATTGTGTACCCCAAAGTGTTGGGGGAAACCCCCTCTTTTTATGTGGTTTTTCAGCAGTTTTTTTTTGGGGGGGGGCCCgtt                                                                                                                                                                                                          >  2:350974/1‑239 (MQ=255)
                             tttggcGAATTGTCGGGGGGATTCGTTGATTTCGGGTTCCTTCCAGACCGGTCCCCGGGAAATTACCTTTTAAGCGTAATTTTAATGGTGATTGGCCTTTACCAATTTCATGGCGGTGAAAACCAACCACGATTTCCAAGGCATTTTTATCCGTATGTGTTACCACAATGGGTTGGGTAAAACCATCTATTTTAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                      <  1:350974/234‑1 (MQ=255)
                             tttGACGAATAGTCGCGGCGATGCTTTGATTCCGCGTTCCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATTATGGGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                      <  1:900960/237‑1 (MQ=255)
                             ttggacgaatgtgcccggggATTCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                      <  1:431048/224‑1 (MQ=255)
                             tATGACGAATAGTCGCTGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGttt                                                                                                                                                                                                      <  1:853980/239‑1 (MQ=255)
                                         tCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCCCCATTTTCACGGCCTTTTTTTCCCTATGTGTGGACACAATTGAGTGGGTTAAAACCTCCATTTCAATTGATTTTTTCAGCAGGTTTTTTTTTGGGGGGGCCCCCTTGTTTGGGttttttgaa                                                                                                                                                                                          >  2:260880/1‑232 (MQ=255)
                                          gggggggttgggTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCAt                                                                                                                                                                                        <  2:725721/228‑1 (MQ=255)
                                          ccctggtgttctttGATTACCGGTACTTTCCAGACCGGTCACCGGCAAATAACCTTTTAAGCCTAATTTTAATGATTATGAGGCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCGCGTTGTTTGGCTTGTAATCac                                                                                                                                                                                        <  2:1007739/227‑2 (MQ=255)
                                                        gtttACGCGGATCTGCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGGGCATTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTgg                                                                                                                                                                         <  1:482581/237‑1 (MQ=255)
                                                         aTTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCAACCACAAATTTCATGGCAGTTTTTTTCGTATGTGTGACCACCATTGGGTTGTGGAAACCATCTTTTTTATTGGATTTTTTTAGCAGGTTTTTTTTGGGTGGGGGCGCCTTGTTTTGGTTTTGATCATTAGGCCTAAtttt                                                                                                                                                                          >  1:461901/1‑236 (MQ=255)
                                                                  tccttgcagaaagtgtctcGGCAAATAACCTTGTAAGCGTAATTTTAATGATTATGTGCCTTTACCTATTGCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTa                                                                                                                                                              <  1:372463/220‑1 (MQ=255)
                                                                   ccTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACCATTTTCCTGGCCTTTTTTTCCCTATGTTTGACCACAATGGGTTGGGGAAAACCATCTATTTTAATTGTTTTTTTTACAAGTTTTTTTTTTGGGGTGGGCCCGTTTTTTTGGTTGTAAATATTAGGAAAAAATTTGTTGtttttta                                                                                                                                                             >  2:372463/1‑236 (MQ=255)
                                                                   ccTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCACGATTTTCCTGGCATTTTTATCCGTATGTGTGGCCACAATGGGGTGGGTAAAACCATCTATTTCAATGGGTTTTTTCAGCAGTTTTTTTTCGGGGGGTGCCACGTTGTTTTGGTTTTTATCATTAGGCTTAAGTTTGCTGTTTTTaa                                                                                                                                                             >  1:409979/1‑238 (MQ=255)
                                                                   ccTTCCAAACAGTGCACCGGAAAATAGCCTTTTAAGCGTAATTTGAATGATGATGAGCCTTTAGCAATTGCATGGCGGTGAAAACCATCCACGATTTCCAGGGCATTTTTATCCGTAGGTGTGACCACAATGGTTTGGGTAAAACCATCTATTTCAATGTATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAAc                                                                                                                                                            <  1:246050/239‑1 (MQ=255)
                                                                        cAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCcaca                                                                                                                                                       <  1:430788/239‑1 (MQ=255)
                                                                              ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTTCATGGCATTTTTTTCCGTATGTTTGGCCACAATGGGTTGGGGAAAACCAACTATTTCAATGGGTTTTTTCAGCAGTTTTTTTTCGGGGGGTGGCCCGTTTTTTTGGTTTTAATCATTAGGGATAAATTGGCTGTTTTTTAACcaccaaaagg                                                                                                                                                 >  1:734141/1‑232 (MQ=255)
                                                                              ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCACCGTTTTCATGGCAATTTTTTTCCTATTTTTGGCCAACAATGGGTGGGGAAAACCATCTATTTTTATTGGTTTTTTTAGCAGATTTTTTTTTGGGGGGGCCCAGTGGTTTGGGTTGTGTACATTTTGCTTAAATTGGGTGGTTTTTAACCCCCAAAAca                                                                                                                                                  >  1:16907/1‑238 (MQ=255)
                                                                              ggTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCACCGCTTTTCCTTGCATTTTTTATTGTATGTGTGAACACAACAGGTTTGGTAAAAACATCCATTTTACTTTAATTTTTTACCAGTTTTTTTTTTGGGGTGGTCACGTTTTTTTGGTTTTGATATATAGGGGTAGATTGGTTGTTTTTTAACCcaacaaaag                                                                                                                                                  >  2:772577/1‑232 (MQ=255)
                                                                                ttcccGGCAAATAACCTTTTAAGCTTAATTTTAATGGTTATGAGCCGTTCCCAATTTCCTGGCGGTGAAAACCATCCCCTGTTTCCATGGCATTTTTAACCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCa                                                                                                                                                <  2:17896/235‑1 (MQ=255)
                                                                                cccccGGCAAATTACCTTTTAAGCTTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGTTTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCa                                                                                                                                                <  2:757131/235‑1 (MQ=255)
                                                                                     ggCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTTAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCAc                                                                                                                                          <  1:326015/239‑1 (MQ=255)
                                                                                          aTAACCTTTTAAGCGTAATTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAACCCTCCCAGATTTCCATGGCATTTTTATCCGTATGTGTTACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGct                                                                                                                                     <  1:491211/239‑1 (MQ=255)
                                                                                              ccTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAAACATCCCCGATTTCCCTGGCATTTTTTTCCGTATGTGTGACCACAAAAGGTTGGGGAAAACCATCCATTTTAATTGGTTTTTTTCACAGGTTTTTTTTGGGGGGGGCCCCCTTTTTTGGGTTGTTTTCTTTTGGGCAAAGGTTGGTGGTTTTTACCCACAGAAAGCAAACCAccaggccttc                                                                                                                                 >  2:725129/1‑230 (MQ=255)
                                                                                                 tttAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTTTTCGTATGTGTGACCACAATGGGGTGGGGGAAACCCACTTTTTTCATTGATTTTTTCTGCCGTTTTTTTTTGGGGGGTGCCCCCTTGTTTTGGTTTTGATCATTTGGCATAAGTTTGGTTGTTTTTATCCCCAGaaaccaaacccacgggccttcccc                                                                                                                              >  2:749444/1‑217 (MQ=255)
                                                                                                          aaTTTTAATGATGATGGGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGc                                                                                                                      <  1:294539/238‑1 (MQ=255)
                                                                                                                     tgatgaGCCTTTACCAATTTCATGGCGGTGAAAACCAGCCACGATTTCCATGGCATTTTTATCCGTATGGGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGa                                                                                                          <  1:1054434/239‑1 (MQ=255)
                                                                                                                             cTTTACCAATTTCAGGGCGGTGAAAACCGTTCACGATTTCCATGGCATTTTTATTCGTATGTTTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATg                                                                                                  <  1:379280/239‑1 (MQ=255)
                                                                                                                                  ccAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCa                                                                                                                                                                                                              <  1:573894/130‑1 (MQ=255)
                                                                                                                                       ttCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCt                                                                                                                     <  1:830094/210‑1 (MQ=255)
                                                                                                                                                           tCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCt                                                                                                                                                                       <  1:57690/143‑1 (MQ=255)
                                                                                                                                                           tCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCt                                                                                                                                                                       >  2:57690/1‑143 (MQ=255)
                                                                                                                                                                   ttCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATCCGGATTTCATTAATCGATTTGAGGGAATCATCTTCCg                                                             >  2:379504/1‑239 (MQ=255)
                                                                                                                                                                       aTGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGTTCATCTTCCGGCa                                                          >  1:247786/1‑238 (MQ=255)
                                                                                                                                                                          gCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTACTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCaa                                                         >  1:112735/1‑236 (MQ=255)
                                                                                                                                                                              ttttttccttttggttgATCATAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                                  <  2:112735/224‑1 (MQ=255)
                                                                                                                                                                              ttttttCCGTATGTTTGACCACCATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAAGCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGGAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCg                                                  <  1:1013704/239‑1 (MQ=255)
                                                                                                                                                                              ttttATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGGTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATCCGGATTTCATTAATAGCTTTGATGCGATAATATTCCGCTAAGCTggtg                                                  >  1:1058178/1‑237 (MQ=255)
                                                                                                                                                                                   tCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGaaa                                             <  1:799017/239‑1 (MQ=255)
                                                                                                                                                                                      gTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGCTAACCTGATGGATCGCCATGCGGATTTCTTTAATCGCTTTGATGCGATCACTTTCGGGCAAGCTGGAGGGCAAttcc                                         >  2:985044/1‑237 (MQ=255)
                                                                                                                                                                                       tATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTAACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGCTGCGATCATCTTCCGGCAAGCTGGCGAAAAATTGc                                         >  2:246505/1‑239 (MQ=255)
                                                                                                                                                                                       tATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCAATGGCTCTTCACGAAAGGGGTTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAAtt                                           >  1:609450/1‑237 (MQ=255)
                                                                                                                                                                                         tgtgtgACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGt                                       >  1:19164/1‑239 (MQ=255)
                                                                                                                                                                                                aCCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                                      >  2:452953/1‑202 (MQ=255)
                                                                                                                                                                                                aCCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGat                                                                      <  1:452953/202‑1 (MQ=255)
                                                                                                                                                                                                aCCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAAt                                  >  1:316732/1‑238 (MQ=255)
                                                                                                                                                                                                aCCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCAATGCGGATTTCATTAATCGCTTTGAGGCGATCATCTTCCGGCAAGCTGGGGAGAAAATGGGTTaaaa                                  >  2:1034749/1‑237 (MQ=255)
                                                                                                                                                                                                  cacaATTGTTTGGGTAAAAACATTTTTTTCAATGGATTTTTGCGGAAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCt                               <  2:16907/239‑1 (MQ=255)
                                                                                                                                                                                                             ttgggtAAAACCTTCTATTTCAATGGATTTTTGCAGCAGTTTTTTGTGGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCTACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAAt                       <  1:367486/238‑1 (MQ=255)
                                                                                                                                                                                                             ttgggtAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCAAGAAATTGCGTTAAATCCTGCGTTaa                        >  1:210650/1‑237 (MQ=255)
                                                                                                                                                                                                              ttggtTAAACCTTTTTTTTTAATTGATTTTTTGAGGAGTTTTTTTTCGGGGTGTGCCACGTTGGTTGGGTTGTAATCCTTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCg                     <  2:461901/237‑1 (MQ=255)
                                                                                                                                                                                                              tgggtAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCg                     >  1:9463/1‑239 (MQ=255)
                                                                                                                                                                                                                ggtAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgtt                   >  1:161676/1‑239 (MQ=255)
                                                                                                                                                                                                                 gtAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTGTGTCGGGTGGTGCACAGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                  <  1:985044/238‑1 (MQ=255)
                                                                                                                                                                                                                  tAAAACCTTCTATTTCAATGGTGTTTTGCAGGAGTTTTTTGTCGGGTGGTGCCACGTTGGTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttg                  <  2:118491/238‑1 (MQ=255)
                                                                                                                                                                                                                  tAAAACCAGTTGTTTCAATGGATTTTTGCAGCAGTTTTTTTGCGGGTGGTGCCAAGTTGTTAGGGGTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAAGCGCTTTGATGCGATCATCTTCCGGCAAGCGGGCGAGAAATTGCGTTAAATCCTGCGTTAAGCgttg                  <  1:361992/238‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                <  1:725129/239‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCgttgtt                <  1:195796/238‑1 (MQ=255)
                                                                                                                                                                                                                     aaCCTTCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCt                                                                                                                     <  1:463576/136‑1 (MQ=255)
                                                                                                                                                                                                                     aaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGCAAGGGGCTAACCTGATGGATCGCCATGCGGCTTTCATTAATCGCTTTGATGCGTTCATCTTCCGGCAAGCTGGCGGGAAATTGCGTTAAATCCTGCGTTAATCctttttca              >  2:244853/1‑237 (MQ=255)
                                                                                                                                                                                                                     aaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGTTGCGATCATCTTCCGGCAACCTTGCGAGAAATTGCGTTAAATCCTTGGTTAATCgttgtggg              >  2:563450/1‑236 (MQ=255)
                                                                                                                                                                                                                     aaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCAAGAAAGGGGCt                                                                                                                     >  2:463576/1‑136 (MQ=255)
                                                                                                                                                                                                                     aaCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCCCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAATTCGGTTTAATCTTGCGTTAATCGTCGTTGc              >  1:196702/1‑239 (MQ=255)
                                                                                                                                                                                                                      aCCATTTTTTTTAATGGATTTTTGCAGGAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAAGCCTTAGGCATAAGTTGGATGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATTGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCa             <  2:409979/239‑1 (MQ=255)
                                                                                                                                                                                                                      aCCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCa             <  1:749444/239‑1 (MQ=255)
                                                                                                                                                                                                                        cgacTATTTCAATTGATTCTTACATCATTTTTGTGGCGGTTGGGGCCACTTTGGTTGGGGTGTAATCATTAGGCATAAGTTGGCTGGTGTTAACCCACAGAACTCAATCCACTGGCTCTTCACGAAAGGGGCTAACCGGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTTCATa           <  1:777909/236‑1 (MQ=255)
                                                                                                                                                                                                                          tCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCTTTGTTCCATAGt         >  1:370362/1‑239 (MQ=255)
                                                                                                                                                                                                                          tCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATatt         >  1:419214/1‑237 (MQ=255)
                                                                                                                                                                                                                          tCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGt         >  1:305016/1‑239 (MQ=255)
                                                                                                                                                                                                                          tCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTAAAACGGTGTTGCAAAGt         >  2:297400/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGGGAAATTGCGTTAAATCCTGCGTTAATCGTTGTGCAAAATTATccccc  >  2:514988/1‑239 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GGTGATTTGATGGCGTCCACGAGCGCGATTATGACGAATAGTCGCGGCGATGCGTTGATTACGCGTACCTTCCAGACAGGTCACCGGCAAATAACCTTTTAAGCGTAATTTTAATGATGATGAGCCTTTACCAATTTCATGGCGGTGAAAACCATCCACGATTTCCATGGCATTTTTATCCGTATGTGTGACCACAATGGGTTGGGTAAAACCATCTATTTCAATGGATTTTTGCAGCAGTTTTTTTTCGGGTGGTGCCACGTTGTTTGGGTTGTAATCATTAGGCATAAGTTGGCTGTTTTTAACCCACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCC  >  NC_000913/634915‑635374

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