Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 13.4% W98C (TGG→TGC 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,5940GC13.4% 123.2 / 9.9 43W98C (TGG→TGCyafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base C (4/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTTGCGATC  >  NC_000913/252299‑252894
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
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)
                                                                                                                                                                                                                                                                                                       gTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTTGCGATc  <  1:498310/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       gTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGATCGAGATAACAACCAAATGCATAAAATTGGGAAAAGGGCAGAAGCGTTTAGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCAGTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTTGCGATc  <  1:161077/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                            
GGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTCGCTGGAAGGAGAAACTCGCGAAATCACAAGTGTGGGTTGCGATC  >  NC_000913/252299‑252894

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