Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,616,197 C→A 5.8% V145V (GTC→GTA yneK → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,616,1970CA 5.8% 229.6 / 2.9 70V145V (GTC→GTAyneKuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base A (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

CAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCACCGCCATTCTCGACAT  >  NC_000913/1615918‑1616472
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
cAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTc                                                                                                                                                                                                                                                                                     >  1:747057/1‑280 (MQ=255)
cAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTc                                                                                                                                                                                                                                                                                     >  2:1053670/1‑280 (MQ=255)
cAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAAATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTAAGACCGGCGTCg                                                                                                                                                                                                                                                                                    >  2:511981/1‑281 (MQ=255)
                      tttAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGAGCGTATTTCAGACCGGCGTCGAGCGGTTACTCtttctgtttct                                                                                                                                                                                                                                                              >  2:1047512/1‑281 (MQ=255)
                                         tATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTACCCGGCCAGGAATGGAAGACACAAAGAAAATATAGAATTTCAGGTAGTAATAAATTTAAGCCCTCTTTATTTAGATAATTTTCGCCAGGATGGAGAAATCAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTCTGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTGCTGAATGATTTTATTGAACag                                                                                                                                                                                                                                           >  2:425204/1‑280 (MQ=255)
                                                 atTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCaa                                                                                                                                                                                                                                    >  2:1011056/1‑280 (MQ=255)
                                                  ttaagagGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGCGCGAGTATTTCAGACCGGCGTGGAGCGCATACTATTTCTGTTGCGGACTGATTTTATAGAACAATTTCAAAt                                                                                                                                                                                                                                   >  2:291356/4‑280 (MQ=255)
                                                  ttCCGAGGACTTCTTTTCGGGGTAACAAAACGCTATACAATAAAAAGGTGCAATTTAAACTCACTTATCCCACCAGGAATGGAAGACCCAATGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTAGTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACTATTTCCAATg                                                                                                                                                                                                                                  <  2:915100/281‑1 (MQ=255)
                                                    cAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATTGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                <  2:724212/281‑1 (MQ=255)
                                                    cAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                <  1:1066402/281‑1 (MQ=255)
                                                    cAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                <  2:548017/281‑1 (MQ=255)
                                                       agGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCaa                                                                                                                                                                                                                             >  2:808267/1‑281 (MQ=255)
                                                       agGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTGTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGGAGTAATAAATTTAAGCCAGATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGAGCGTATTTCAGACCGGCGTCGAGCGTGTAATCTTTCTGTTTCGGAACGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                >  1:915101/1‑278 (MQ=255)
                                                                 ttCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTAGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGtt                                                                                                                                                                                                                   >  2:1131842/1‑281 (MQ=255)
                                                                  tCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGtt                                                                                                                                                                                                                   <  2:747057/280‑1 (MQ=255)
                                                                            gACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAag                                                                                                                                                                                                        <  1:843126/281‑1 (MQ=255)
                                                                                cGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCa                                                                                                                                                                                                                                     <  1:975927/248‑1 (MQ=255)
                                                                                cGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCa                                                                                                                                                                                                                                     >  2:975926/1‑248 (MQ=255)
                                                                                cGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAgc                                                                                                                                                                                                     >  2:762422/1‑280 (MQ=255)
                                                                                 gCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAgcgc                                                                                                                                                                                                   <  1:1053671/281‑1 (MQ=255)
                                                                                   tataCGATGAGAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAGGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCAt                                                                                                                                                                                                 <  1:128428/281‑1 (MQ=255)
                                                                                     taCAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATac                                                                                                                                                                                               <  1:762422/281‑1 (MQ=255)
                                                                                     taCAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTAGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATac                                                                                                                                                                                               >  2:439266/1‑281 (MQ=255)
                                                                                        aaTGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcacc                                                                                                                                                                                            <  1:808268/281‑1 (MQ=255)
                                                                                         aTGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGTAATTGAAGACACAAAGAAAATATTTAATTTCAGGTAGTAATAAATTTAAGCCCTATTGTTTTAGATAATTTTCGCCATGGTGGAGAAATTAATATTTTTGGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGGCTAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcacca                                                                                                                                                                                           <  2:996605/281‑1 (MQ=255)
                                                                                                     aGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTAGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAAccc                                                                                                                                                                               <  1:1131843/281‑1 (MQ=255)
                                                                                                      gTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCAt                                                                                                                                                                                                 >  1:691741/1‑262 (MQ=255)
                                                                                                      gTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCAt                                                                                                                                                                                                 <  2:691741/262‑1 (MQ=255)
                                                                                                                aCTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCg                                                                                                                                                                      <  1:788086/279‑1 (MQ=255)
                                                                                                                 cTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGCTCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGa                                                                                                                                                                    >  2:201695/1‑280 (MQ=255)
                                                                                                                 cTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcacc                                                                                                                                                                                            <  2:275350/256‑1 (MQ=255)
                                                                                                                 cTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcacc                                                                                                                                                                                            >  1:275350/1‑256 (MQ=255)
                                                                                                                          ccAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGTTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGACCACCATAcc                                                                                                                                                          >  2:301139/1‑281 (MQ=255)
                                                                                                                              gAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTg                                                                                                                                                      <  1:1047513/281‑1 (MQ=255)
                                                                                                                              gAATGAAAGACACAAAGAAAATATTGAATTTCAGGTAGTAAGAAATTTAAGCCCTATTTATATAGATAATTTTCGCCATGATGGAGAAATTCATATTTTGTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTAGGATCACCATACCGCTg                                                                                                                                                      <  2:311407/281‑1 (MQ=255)
                                                                                                                                         acaAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGAttactta                                                                                                                                           >  1:1119199/1‑281 (MQ=255)
                                                                                                                                                         aTTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGc                                                                                                                           >  2:108260/1‑281 (MQ=255)
                                                                                                                                                         aTTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAAATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGGATGATTTTATTGAACAATTTCCGATAATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGc                                                                                                                           <  1:291356/281‑1 (MQ=255)
                                                                                                                                                          tttCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTAGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCt                                                                                                                          <  1:439266/281‑1 (MQ=255)
                                                                                                                                                             cAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGt                                                                                                                       >  2:731711/1‑281 (MQ=255)
                                                                                                                                                                           tttAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCa                                                                                                                                                                                                  <  1:243244/192‑1 (MQ=255)
                                                                                                                                                                           tttAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCa                                                                                                                                                                                                  >  2:243244/1‑192 (MQ=255)
                                                                                                                                                                               aGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTcg                                                                                                     <  1:1011057/281‑1 (MQ=255)
                                                                                                                                                                                      tttatttaGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                >  2:995588/1‑151 (MQ=255)
                                                                                                                                                                                      tttatttaGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGAt                                                                                                                                                                                                                                <  1:995589/151‑1 (MQ=255)
                                                                                                                                                                                      tttatttaGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCaaa                                                                                                                                                                                                          >  2:499250/1‑173 (MQ=255)
                                                                                                                                                                                      tttatttaGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCaaa                                                                                                                                                                                                          <  1:499250/173‑1 (MQ=255)
                                                                                                                                                                                           ttaGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTg                                                                                          >  2:313866/1‑280 (MQ=255)
                                                                                                                                                                                             aGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCa                                                                                       <  1:731711/281‑1 (MQ=255)
                                                                                                                                                                                                           atgataGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCATGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTa                                                                         <  1:128214/281‑1 (MQ=255)
                                                                                                                                                                                                                  agaAATTAATATTTTTTGCGCTCACAACCCAAAGCCGGTCACTATGGGGCGCTTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAAt                                                                   <  1:511981/280‑1 (MQ=255)
                                                                                                                                                                                                                      aTTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTTCCCCGATTATGGAATAAct                                                               >  1:572988/1‑280 (MQ=255)
                                                                                                                                                                                                                        tAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcaccac                                                                                                                                                                                          >  1:708110/1‑155 (MQ=255)
                                                                                                                                                                                                                        tAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATAcaccac                                                                                                                                                                                          <  2:708110/155‑1 (MQ=255)
                                                                                                                                                                                                                           tatTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAgg                                                         >  2:1075780/1‑281 (MQ=255)
                                                                                                                                                                                                                                tttGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGt                                                                                                                                    <  1:709116/201‑1 (MQ=255)
                                                                                                                                                                                                                                tttGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGt                                                                                                                                    >  2:709116/1‑201 (MQ=255)
                                                                                                                                                                                                                                                       gTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCa                                                                                                                                                                                                                                     <  2:1012504/81‑1 (MQ=255)
                                                                                                                                                                                                                                                       gTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCa                                                                                                                                                                                                                                     >  1:1012505/1‑81 (MQ=255)
                                                                                                                                                                                                                                                        tCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGCTTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCCAAGCTGGTTCGTCAATCAAGAAcc                             >  2:116755/1‑280 (MQ=255)
                                                                                                                                                                                                                                                        tCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACACTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATcc                             >  1:424433/1‑280 (MQ=255)
                                                                                                                                                                                                                                                            tATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAAc                         <  1:9750/280‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgcgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTTTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATc                   >  1:712807/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                   cgcgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCAc                 <  1:148354/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgcgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCAc                 >  2:148354/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                    gcgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCAc                 >  1:820632/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                    gcgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATCATGGAATAACTCTGAGGTTCACCGATGGTTCGCCAATAAAGAACCTAAATTGATTAcc                >  1:228470/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                     cgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGt                                         >  1:26617/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                     cgTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGt                                         <  2:26617/255‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCACCGCCATTCTCGACa   >  1:245416/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                   gCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGCTCCTAACTTGATCACCGCCATTCTCGAcct  >  1:853671/1‑278 (MQ=255)
                                                                                                                                                                                                                                                                                   gCGTCGAGCGTGTACTCTTTCTGTTTCTGAACGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCACCGCCATTCTCGACAt  >  1:864626/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                   
CAACTCAGCCACAAAAAAGATATTTAATAACCTTGACAACTTATCATATATTCAGAGGACTTCTTTTCGGGGTAACGACACGCTATACAATGAAAAGGTTCAGTTTAAACTCACTTATCCCGCCAGGAATGGAAGACACAAAGAAAATATTGAATTTCAGGTAGTAATAAATTTAAGCCCTATTTATTTAGATAATTTTCGCCATGATGGAGAAATTAATATTTTTTGCGCTCCCAACCCAAAGCCTGTCACTATGGGGCGCGTATTTCAGACCGGCGTCGAGCGTGTACTCTTTCTGTTTCTGAATGATTTTATTGAACAATTTCCAATGATCAACCCTGGTGTTCCCATCAAAAGAGCGCATACACCACATATTGAACCCCTGCCTTCGGATCACCATACCGCTGCAGATTACTTACGCCAGTTTGATTTGCTTGTCCTGAATTTTATCTCTCGCGGTAATTTTGTCATACTCCCCCGATTATGGAATAACTCTGAGGTTCACAGATGGTTCGTCAATAAAGATCCTAACTTGATCACCGCCATTCTCGACAT  >  NC_000913/1615918‑1616472

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