Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,310,770 A→G H32H (CAT→CAC yciI ← protein YciI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,310,7700AG98.1% 176.4 / ‑4.3 54H32H (CAT→CACyciIprotein YciI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/26);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

AACACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTG  >  NC_000913/1310526‑1311015
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
aaCACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGt                                                                                                                                                                                                                                                 >  2:669796/1‑251 (MQ=255)
   acGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGAAGTAAc                                                                                                                                                                                                                                              >  2:628849/1‑251 (MQ=255)
        ttcAAAAACCGTCCCCTTAAATATTGAATCTCTATTCATAACTCTTTCTTAAATGTTTTCACTGAAACGTGGTTATAGACTCCTGCCGCTACGTACGGGTCAGCATCTGCCCAGGTCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTa                                                                                                                                                                                                                                          <  1:655935/248‑1 (MQ=255)
             aaaCCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTg                                                                                                                                                                                                                                    >  2:676309/1‑251 (MQ=255)
             aaaCCGTGCCATTAAATATGGAATCTATATACAGAACCCTTTCTTAAATGGTTTCACTGAAACGTGTTAAAAGACCCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTg                                                                                                                                                                                                                                    <  2:72124/251‑1 (MQ=255)
                    ccccTTAATTATTGAATCTCTATTCAGAACCCTTTCTTAATTGGTTCCACTGAAACGGGTCCATAGACACCTGCCGCTTCGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAtg                                                                                                                                                                                                                             <  2:343624/250‑1 (MQ=255)
                        ttAAATATTGAATATCTATCCAGAACACTTTCTTAAATGGTTTCACTGAAACGCGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTTGCTCGTTACTGTCTACCGCTGGCATTGGTACAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTCAACGTGCTCAATGTGcc                                                                                                                                                                                                                         >  2:321969/1‑251 (MQ=255)
                           aaTATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCgg                                                                                                                                                                                                                       >  1:280186/1‑250 (MQ=255)
                                 accgctcttcaaatctcacTTTATTAAAGGGTTTACCTGAAACGTGTTAATAGACGCCGGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACg                                                                                                                                                                                                                <  2:633434/235‑1 (MQ=255)
                                            cAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCCTCGGCCCAGGACTGCGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGCGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGCCCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAACGGCGttt                                                                                                                                                                                                     >  1:37993/1‑251 (MQ=255)
                                                 cacTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAACCGTGCTAAATGTGCCGGACGAACg                                                                                                                                                                                                                >  1:633434/1‑235 (MQ=255)
                                                      tctaaaagggTTCCACTGAAACGTGTTCAGAGTCGCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGGTGAGCCAGTAAATCCCGCAGCCCCTGGAGCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTGTTTCGAGGCt                                                                                                                                                                                            <  2:452587/242‑1 (MQ=255)
                                                         tAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTAt                                                                                                                                                                                          <  1:625285/249‑1 (MQ=255)
                                                                tttCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGt                                                                                                                                                                                 >  1:185395/1‑251 (MQ=255)
                                                                 ttCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGt                                                                                                                                                                                 <  1:228808/250‑1 (MQ=255)
                                                                  tCACGGAACCGTTTTCATAGCCTCCTGCCGCTACGGACGGGGCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTc                                                                                                                                                                                <  2:604373/250‑1 (MQ=255)
                                                                   cACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCtt                                                                                                                                                                              >  1:555413/1‑251 (MQ=255)
                                                                      tGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGCGGCTATCAGCTTTGTCTTg                                                                                                                                                                             >  1:216703/1‑249 (MQ=255)
                                                                           cgggtTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTaa                                                                                                                                                                      <  1:17063/248‑1 (MQ=255)
                                                                           cGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTa                                                                                                                                                                       >  1:269465/1‑250 (MQ=255)
                                                                                 cATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGGGCGTAAATAAc                                                                                                                                                                 >  1:553577/1‑250 (MQ=255)
                                                                                       cTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAACTGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATacaa                                                                                                                                                           >  2:40372/1‑250 (MQ=255)
                                                                                         ccTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATacaac                                                                                                                                                          <  1:669796/249‑1 (MQ=255)
                                                                                            gCCGCTACGTACGGGTCCGCATCGGCCCAGGCCTGTGCTGCTGCCAGCGATTCAAATTCATCCATAACGGTTGCGCCAGTAAAGCCCGCAGCCTCTGGATCGTTTCTGTCTAACGCTGGCATTTTACCAGCTGTCAACAATCGACCTTCATTGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTGGTCTTGAGCTTAAATAACATACAACACg                                                                                                                                                       <  2:540301/249‑1 (MQ=255)
                                                                                                 tacttccgGGTCAGCATCGTCCCATGTCTCAGCTGTTTCCATCGATTCATATTCCGCATTAACGTTTGAGTCCTTAAACACTGCAGTCCCTGGTACGTTTCTATCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCCTCATCTTGCACTACCTGTAACCGCGTTAAATGTGCCTTACGCACGGAATGTCGTTTTTCGATGCTATCAGCTTTGTCTTGATCGTAAATAACATaccccacgcccaca                                                                                                                                                 <  2:543023/244‑13 (MQ=255)
                                                                                                                     cccAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTg                                                                                                                            <  1:53009/251‑1 (MQ=255)
                                                                                                                            cTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGAACAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTCTCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGTGAAAGTTGTAAAGt                                                                                                                      >  1:331847/1‑250 (MQ=255)
                                                                                                                             tGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTAc                                                                                                                    >  1:176110/1‑251 (MQ=255)
                                                                                                                              gAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTAc                                                                                                                    <  2:185395/250‑1 (MQ=255)
                                                                                                                               aGCTGCTTCCAGCGATTCAAATTCAGAAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGAAACTGTAAACGTGCTAAAAGTGCCGGACGAACGGAACGGCGTTTTTAGAGGCTATAAGCTATGTCTTGAGCGTAAAAAACAAAAAACCCGGGAACAACTCCTTGTTCGGGCAAGTTGTaaataaagt                                                                                                                  >  2:494468/1‑245 (MQ=255)
                                                                                                                                                  aaTTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGCCCTTCATCGTGCAGTAACTGTAACCGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTCAAGAACGTTATTTGAAAGGGCGa                                                                                                   >  2:477671/1‑247 (MQ=255)
                                                                                                                                                      cagaaaaaacagTTGAGCCAGTAAATCCCTCAGCTGCAGGATCGTTACTGGCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCGTCATCGTGCAGTAAATGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTGGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGAGCTg                                                                                            <  2:141951/239‑1 (MQ=255)
                                                                                                                                                      cAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGc                                                                                           >  1:233267/1‑251 (MQ=255)
                                                                                                                                                       aGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACGGTCTACCGCTGGCATTGGACCAGCTGGAAACAAACGACCTTTATCGTGCAGTAACTGTAAACGTGCTAAATGGGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCa                                                                                          <  2:516198/251‑1 (MQ=255)
                                                                                                                                                             aaCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGa                                                                                     <  1:40372/250‑1 (MQ=255)
                                                                                                                                                             aaCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACCCGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGa                                                                                     >  1:80242/1‑250 (MQ=255)
                                                                                                                                                                 gTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAgat                                                                                <  1:92906/251‑1 (MQ=255)
                                                                                                                                                                          ggaaaTCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCttt                                                                       <  2:329810/249‑1 (MQ=255)
                                                                                                                                                                          gTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCttt                                                                       <  1:676309/251‑1 (MQ=255)
                                                                                                                                                                            aaaTCCCGCAGCCCCTGTAATGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGt                                                                     <  2:499635/251‑1 (MQ=255)
                                                                                                                                                                             aaTCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTg                                                                                                                            >  2:85869/1‑195 (MQ=255)
                                                                                                                                                                             aaTCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTg                                                                                                                            <  1:85869/195‑1 (MQ=255)
                                                                                                                                                                                       ccccTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGc                                                          <  1:247550/251‑1 (MQ=255)
                                                                                                                                                                                            ggagcGTTCCTGTCTAACGCAGGCATTGGACCAGTTGGCAACAAACGACCGTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTGTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACTCGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTATGTCCATGCATa                                                       <  2:303732/245‑1 (MQ=255)
                                                                                                                                                                                             gATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATaa                                                      >  1:39024/1‑249 (MQ=255)
                                                                                                                                                                                                   tACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAGGTAAggg                                              >  1:143124/1‑251 (MQ=255)
                                                                                                                                                                                                                tGGCATTTGACCAGCTGTCAAAAAACGACCTTCATCGTGCAGTGACTGTAAACGTGCTAAATGTGCCGGACTAACGGACAGGCGTTTTTCGATGCTATCAGCTTTATCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGTGAAATTTGTAAAGTACGTTATTAGAACGGACGACGATCTGCAACGCAAAGATGATGTTTTTGTTAAGTCCATGCTTAAAGTAAgaataataccgcaaa                                  >  2:25526/1‑236 (MQ=255)
                                                                                                                                                                                                                 ggCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCaaa                                 <  1:668762/250‑1 (MQ=255)
                                                                                                                                                                                                                     ttGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATg                             >  1:227579/1‑250 (MQ=255)
                                                                                                                                                                                                                        gACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAACCGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAACGTTGTAAAGTACGTTATTTGAAAGGGCGCAGATCTGCAACGGAAAGATGATGGCTTTGTTAATGCCATGCATAAAGTAAGGGTAATTAAGCCAAAAATGACAt                         >  2:468873/1‑251 (MQ=255)
                                                                                                                                                                                                                                gTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGTGTAATTACGCCAAAAATGACATTTTCAc                   <  1:240894/249‑1 (MQ=255)
                                                                                                                                                                                                                                gTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCAct                  <  1:468873/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAACTGACATTTTCActgatcctg           >  2:243335/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCg       >  1:365652/1‑251 (MQ=255)
                                                                                                                                                                                                                                                catcgtGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTGTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTg  <  2:280186/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
AACACGGTTGCAAAAACCGTGCCCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCTTG  >  NC_000913/1310526‑1311015

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