Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 10.3% A97P (GCA→CCA)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5890GC10.3% 120.8 / 6.1 39A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base C (2/2);  total (20/19)
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.50e-01

AATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTT  >  NC_000913/252294‑252888
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
aaTAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAACTGAGAATAATCCATTTCCGCCACAGTTGCGCAAATCCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGCCAGGGGGCGTTAGATGACAAAGCGTgt                                                                                                                                                                                                                                                                                                        >  1:361304/1‑300 (MQ=255)
     ggATTCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                                                                                                                                                                                                                                   <  2:485836/301‑1 (MQ=255)
         gCGGGTATTCAAAACAAAACTTATTCGACTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                                                                                                                                                                                                                               <  2:574692/301‑2 (MQ=255)
              tATTCAAATCAAAACTTATTCGACTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaacg                                                                                                                                                                                                                                                                                          <  1:488565/301‑7 (MQ=255)
                     aacaaaacTTATTCGCCTGCAACGTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                                                                                                                                                                                                                   <  2:361788/301‑1 (MQ=255)
                        aaaacTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTACGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTCCCAATTGCCATTCTGAAACCt                                                                                                                                                                                                                                                                                >  1:218176/1‑301 (MQ=255)
                                           cTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCAAAAACTGGATc                                                                                                                                                                                                                                                             >  2:504959/1‑301 (MQ=255)
                                            ttACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                                                                                                                                                                            <  2:465249/301‑1 (MQ=255)
                                               aagcagAGGAACTTGATGCGTTAACGGCGGATTTTATGTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                                                                                                                                                                                                                                         <  1:247807/300‑1 (MQ=255)
                                                agcagAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                                                                                                                                                                                                                        <  2:439868/301‑1 (MQ=255)
                                                                 ggTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                                                                                                                                                                                                                                       <  1:217895/300‑1 (MQ=255)
                                                                               ttgtttttctttatgtggttcGGTGTTTTTCCAGATATATTTGGTCGCGATGCACTCTACGACGACGCCTTTTCCTGGCCATTAAGCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCAGTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGGAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                                                                                                                                                                                                                         <  1:485569/281‑1 (MQ=255)
                                                                                cgctttccctacacgacgctcttccgatctcAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTGCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                                                                                                                                                                        <  1:439238/271‑1 (MQ=255)
                                                                                                    ggTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtt                                                                                                                                                                                                    <  2:327115/301‑1 (MQ=255)
                                                                                                              cAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCCTTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGatcggaagagcacacgtctgaactccaa                                                                                                                                                                                          >  2:439802/1‑274 (MQ=255)
                                                                                                                                       ctACGACGACTCCTTTACCTGGCCATTAATCAAATTTCAGCTCGTTGCTCATATTCATCTGGCAAATGACAATAATCCATCTCCGCCACAGTTGCGCCAATCCAGCAGAACGAATGACGAAGCGCATTTGCTATATTGTCAAGGGGCGTTTGCTGAGCAAGCATGCATGCTCATTGCCATTCGGAAACCTGAACCTCATAACCTGGCTCGAACTCACAACCAAACACATAAACATGCGAAAATGCCAGAAGCGTATCGCATGCATTTTAAAAATAAAATAATAAATAACATAACAAAAAGa                                                                                                                                                                 >  1:428144/1‑301 (MQ=255)
                                                                                                                                            ctgacTCTTTTACCGGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAAGTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTa                                                                                                                                                            <  1:179507/299‑1 (MQ=255)
                                                                                                                                                   cTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAAAAAGAAACTAACAGCCt                                                                                                                                                     >  2:451666/1‑301 (MQ=255)
                                                                                                                                                           ggtcATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGATAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGAttt                                                                                                                                             <  1:40138/298‑1 (MQ=255)
                                                                                                                                                                                gAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAACTGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCCACTATGCGCTAtttt                                                                                                                        >  1:236546/1‑301 (MQ=255)
                                                                                                                                                                                  gTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCa                                                                                                                      >  2:259154/1‑301 (MQ=255)
                                                                                                                                                                                    tGCTTATTTTCATCTGGCAAATGAGAATAATCCATTTCCTGTACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGTGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCAtt                                                                                                                    <  1:451102/301‑1 (MQ=255)
                                                                                                                                                                                    tGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCAtt                                                                                                                    <  1:258886/301‑1 (MQ=255)
                                                                                                                                                                                          tatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGAAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTAACAGCCTGGCGATTTTCAGCAACTATGCGATATTTTAATTAGAGCg                                                                                                              >  2:232025/1‑301 (MQ=255)
                                                                                                                                                                                                 atGGCAAATGAGAATAATGCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGa                                                                                                       <  1:549804/300‑1 (MQ=255)
                                                                                                                                                                                                  tGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGAt                                                                                                      >  2:161228/1‑301 (MQ=255)
                                                                                                                                                                                                         tggattattctcATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGACGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGTTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGcc                                                                                               >  2:295455/12‑301 (MQ=255)
                                                                                                                                                                                                                aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAAAGGTTACGATAACAGCCAGTCAga                                                                                        >  2:81755/1‑300 (MQ=255)
                                                                                                                                                                                                                    cATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGTTGACCACCAGTCAGCatta                                                                                    >  2:438224/1‑301 (MQ=255)
                                                                                                                                                                                                                                 tttGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCAcaacaa                                                                       <  2:236812/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCGGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCAAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAAactcc                                       >  1:543201/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                  tcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGc                                      <  2:153082/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ctttccctacacgacgctcttccgatctaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGga                              <  1:314131/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGTTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTcgcg                    <  1:295051/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCgaga                   >  2:141488/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGaa                  >  2:495901/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      agaatggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTCTATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCCGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATCATTCACCACAACAAATTTCCGCCTGTGCGCAGATTAACGAATCTCGCAGGAAAAAAAAACTCGCGaa                  >  1:482254/10‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             tgatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAAACCGCGAAATCACaa           >  1:642515/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACaa           >  2:223254/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                 gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACCCGCAAAATCACCAgtgt       >  1:23199/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  agcaaccaTGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGtgtg      >  2:565892/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAgtgagagt   >  2:498957/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                                     aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGt   <  2:643392/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGgaagatcggaagagcacacgtctgaactca  >  2:314535/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                           
AATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTT  >  NC_000913/252294‑252888

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