Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,612 T→C 18.1% I42V (ATA→GTA)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,6120TC18.1% 100.6 / 2.7 54I42V (ATA→GTA) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base C (3/7);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAA  >  NC_000913/2596388‑2596837
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
ctttCCCCGGTTCATGTCCAGCTTCGGTTGTTGATTGGTCTTACCCATCTTTCCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAg                                                                                                                                                                                                                                 <  1:77901/226‑1 (MQ=255)
                 ccggCTTCGGTTGCTGATAGGTCTTACTCATATTCACGATCGTTTGTCCCCAGGTGTCATCCGTGATTGCCAGAATGTCAGAGTTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGTACGGTAATCAACCCGGTTAAATCAACACCAAAGTAGCCTGCAACCAGCACGTCTATCAGTAAAATCAGCCCGCCtt                                                                                                                                                                                                                <  2:560617/225‑1 (MQ=255)
                    gCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTa                                                                                                                                                                                                              >  2:382298/1‑227 (MQ=255)
                    gCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAACCAACCCCCCTTTAc                                                                                                                                                                                                             >  1:525768/1‑228 (MQ=255)
                    gCTTCGGTTGCTGATAGGTCTCACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCACACCCCCTTAAc                                                                                                                                                                                                             >  1:39494/1‑228 (MQ=255)
                          gTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTACCAGTAAAATCAGCCCCCCTTTACCCCtt                                                                                                                                                                                                        >  2:158065/1‑227 (MQ=255)
                                           cccATCTGCTGGATCTGTTTGCCCCAGTTGTCTTCCGTGTTTGCAAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATGTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAACCCGGCTGCCCGGTCATCAACCAGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGa                                                                                                                                                                                      <  2:215739/228‑1 (MQ=255)
                                            ccATCTTCTCGAACTGGTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAc                                                                                                                                                                                     <  2:39494/228‑1 (MQ=255)
                                                cttctCGAACTGTTGTTCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCCGCACGACTACAAGTAAAATAAGCCCGCCTTTACCGCTTGTCAGGCGAAAACCGGGACCGcc                                                                                                                                                                                 >  2:553107/1‑228 (MQ=255)
                                                cttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATAAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcgac                                                                                                                                                                                  >  2:287778/1‑224 (MQ=255)
                                                              gtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACAATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacca                                                                                                                                                                    >  2:606520/1‑227 (MQ=255)
                                                                    aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATATTAGCCTGCAACCAGCACGACTATCAGTAAAATCATCCCGCCTTTACCGCTTGGCCGGCGCAAACCGGGACGGCCCGAAGATGAACCACCaaagg                                                                                                                                                             >  2:556992/1‑227 (MQ=255)
                                                                    aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccaccaaga                                                                                                                                                              >  2:479897/1‑224 (MQ=255)
                                                                      gtgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGc                                                                                                                                                             >  2:352962/1‑226 (MQ=255)
                                                                            tcggtgggtggcaaatcccccGGGGGGATTTTTGGGGTTTCGGCTCAATTTGGCATAATTGACCGCGTGGATTTGTGTGGGGAAACCGGCACCCCGGTCACCAACCCGGTTAAAGCAACCCCATAGTAGCCGGCAACCAGCCCGACCACCAGAAAAACCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgcgc                                                                                                                                                     <  1:352962/209‑1 (MQ=255)
                                                                            cccGTGGTTGCAAAAAGCACCGAGGTTAATTTTGCGGTTTCGGTTTCATTGGGGCTAATTGCCCTGGTTGAGTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAAGCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgcgc                                                                                                                                                     <  2:378732/227‑1 (MQ=255)
                                                                                tGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGCCCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAAATGGACCACCAAGGCTGTTGCGCCTGt                                                                                                                                                 >  2:433773/1‑228 (MQ=255)
                                                                                   ttGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCTGCTGCCCGGTCATCAACCCGGTTAAATCAACACCCTAGAAGCCTGCAACCAGCACCACTATCAGGCAAATCAGCCCGCCTTTACCGCTTGGCCGGCGAAAACCGGGCCCGCCCAAAAATGGACCACCACAGCTGTTGCGCCTGTCtg                                                                                                                                              >  1:165697/1‑227 (MQ=255)
                                                                                                     tGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAAGTGGACCAACAGAGCTGTTGCGCCGGCCTTCACCTTTTTCACTTTCa                                                                                                                             >  1:513393/1‑227 (MQ=255)
                                                                                                           ttGGGGTTTCGTTTTAATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGcc                                                                                                                      <  2:493463/228‑1 (MQ=255)
                                                                                                                 cTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGcc                                                                                                                <  2:513393/228‑1 (MQ=255)
                                                                                                                  ttCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGcc                                                                                                                <  2:525768/227‑1 (MQ=255)
                                                                                                                   tctcttctttttGGGTTATTGGCCCGGGTTAATTTTTGTTGGGAACCCGGTTGCCCGGACACCACCCCGGTTAAATCAACCCCATAGTAGCCTGCAACCAGCACTACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCaa                                                                                                              <  1:382298/219‑1 (MQ=255)
                                                                                                                         ccgttggggATAATTGACCGCGTTGATTGTTGTTGGTAAACCGGCTGCCCGTTCATCAACCCGGTTAAATCACCACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATa                                                                                                        <  2:288704/222‑1 (MQ=255)
                                                                                                                           ctttGGACTAATTTCCCGCGTTAATTGTTCTTGGGAACCCGGCTCCCCGGTCCCCCACCCCGTTAAACCAACCCCATAGGAGCCGGCAACCAGCCCAACAACCAGTAAAACCAGCCCTCCTTTACCGCTTTGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGa                                                                                                      <  1:606520/227‑1 (MQ=255)
                                                                                                                            gttGGGATAATCGCGCGCTGTGATGGTTGTTGGGAAACCGGCGTCCCGGTCATCAACCCGATTAAAACAACCCCAGATTAGCCTGCAACCAACACAACGACAGGTAAAATCAGCGCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACGTTCACGTCGCCCTTGCCAACGCATAGa                                                                                                      <  1:158065/226‑1 (MQ=255)
                                                                                                                                gtgtaatgacccGGGTTGAGTGTTATTGGGAACACGCACACCCGTTCAAGAACCCGGTTAAAACAACCCCATTGTAGCCAGCAAGCAGCCCAACTAGCAGTAAACTCAGCCCGCCTTGACCGCGTGGCAGGCGAAAACCGGGACCGCCCATAGATTGACCACCAGAGCTGTTGCGCATGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCt                                                                                                 <  1:678655/218‑1 (MQ=255)
                                                                                                                                ggCTAATTTACCGCGTTAATTGTTTATGAGAAACCGGCGGCCCGGTCAACCACCCGGTTAAACCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGATGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCt                                                                                                 <  2:558981/228‑1 (MQ=255)
                                                                                                                                     aTTGACCGGGTTGATGTTTGTTGGGAAACCGCGCGCCCGGTCAACAACCCGGTTAAATCAACCCCATAGGAGCCTGCAACCAGCACGACAATCAGTAAAACCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTcaaca                                                                                            <  1:373965/228‑1 (MQ=255)
                                                                                                                                       tgcccGGGTTTATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACCCCATAGTAGCCTGCAACCAGCACGATTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAat                                                                                          <  1:479897/225‑1 (MQ=255)
                                                                                                                                           cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCCTAGATGGACCCCCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGCCGCCCTTGCCAACGCATAGAAACCTCCACAatatat                                                                                     >  2:18436/1‑228 (MQ=255)
                                                                                                                                           cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCACCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCAATAGATGGACAACCAGAGCTGTTGCGCCTGCCTTCAACATTGTCACTTTCCCGTCGCCCTTGCCAAGCCTTAGATACCccaccacataat                                                                                      >  2:678655/1‑222 (MQ=255)
                                                                                                                                              gTTGATTTTTTGTGGGAAACCGGAGTGCCGGTAAGCAACCCGGGTAAATCAACCCCATGTTAGCGTGCAACCAGCCAGACTACAAGTAAAACAAGCCCGGCTTGACCGCTTGGCAGGCGGAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATCCCTCAAAAATATATTCa                                                                                  <  1:68313/228‑1 (MQ=255)
                                                                                                                                               ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCATTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGGACCTGGCTTCAACATTGTCACTTTCACGTCGCCCTTGCCCACGCCTAGAAACCTCAACAATATATTTAt                                                                                 >  2:667315/1‑228 (MQ=255)
                                                                                                                                                      ttgttgGGAACACGGCTGCCCGGTCATCAACCCGGTTAAATCACAACCATAGTAGCCTGCAACCAGCACAACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGaa                                                                          <  2:683323/228‑1 (MQ=255)
                                                                                                                                                      gttttgGAAATCCGGCTGCCCGGTCCCCATCCCGTTTAAACCACCACCATAGTAGCGTGAAACCAGCACAACTACCAGTAAAACCAGCGCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGaa                                                                          <  1:359460/225‑1 (MQ=255)
                                                                                                                                                        gttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGACAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTCTCACTTTCACGACGCCCTCTCCAACGCCTAGATACCCCCACAATATCTTCATTATAGAAAt                                                                        >  1:693723/1‑228 (MQ=255)
                                                                                                                                                         ttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTg                                                                                                                                      >  2:35097/1‑166 (MQ=255)
                                                                                                                                                                 ccGGCTGCCCGGTTATCAACCCGGTTAAATCAACCCCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTGACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAgg                                                               <  2:591977/228‑1 (MQ=255)
                                                                                                                                                                    gCTGCCCGGTCATCAACCCGGTTAAATCAACACCTTAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCAGCCCTTACCGCTTGGCAGGCGAAAACCGGGAACGCCCATAGATGGACCACCAGAGCTGTTGCGCCGGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATAcc                                                                                                  >  2:251572/1‑191 (MQ=255)
                                                                                                                                                                      cgcgcagtacacaACCCGGTTAGAACAGCCCCCTAATAGCCCGCAACCAGCACGACTATCAGTAAAACCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCGCAACAATATATTCATTATAGAAATGATCGTAGGTGGtt                                                          <  2:290203/217‑1 (MQ=255)
                                                                                                                                                                       ccccGGGCAACACTCCGGTTAAACCACCCCCATGGGAGCCTGCACACGGCCCGACATACAGTAAAACCCGCCCGCCATTCCCGCTTGGTAGGCGAAATCCCGGACCGCCCCTAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGtt                                                          <  1:433773/226‑1 (MQ=255)
                                                                                                                                                                                 caACCCGGTTAACTCAACACCATAGTAGCCTGAAACCAGCACGACTATCAGTAAAACCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGAt                                               <  2:693723/228‑1 (MQ=255)
                                                                                                                                                                                       ggTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCAtcg                                                                                                       >  2:258770/1‑165 (MQ=255)
                                                                                                                                                                                              tCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAg                                  <  2:129541/228‑1 (MQ=255)
                                                                                                                                                                                                  cacCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATCGACCACCAAAGCTGTTGCGACTGTCTTCAAAATTGTCACTTTCACGTCGCCCTTTGTAACGCATAGACAACTCACAATAACATCCATTTAAGAAATGCCCGTCGGGAGATCAGAGGAAAATTACCAAAGGAAAACgag                              >  1:176136/1‑226 (MQ=255)
                                                                                                                                                                                                   cccATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGa                             <  2:114325/227‑1 (MQ=255)
                                                                                                                                                                                                            gCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATCGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCAccc                                                                                                                     >  2:254573/1‑132 (MQ=255)
                                                                                                                                                                                                                 cAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCCGAGCGGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCACCGCATCGCTACCTCACCCATATATTCATTATAGCAATGATCGTAGGTGGTTCAGAGGAACATAACCCACGGAAAACGGGAAGAAAAgggggggg               >  1:257847/1‑228 (MQ=255)
                                                                                                                                                                                                                 cAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGAGGAAATGGCggagg               >  2:520801/1‑228 (MQ=255)
                                                                                                                                                                                                                       gccccaCTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCCAAAACCGGGTCCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGa          <  2:176136/222‑1 (MQ=255)
                                                                                                                                                                                                                       gCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGCGGTTCAGAGGAAGATTACCCCAGGAAAGCGGGGTGCAAAGGCGGGGGaaaggg         >  1:241725/1‑224 (MQ=255)
                                                                                                                                                                                                                           gACTACCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCTCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGtt     <  2:165697/228‑1 (MQ=255)
                                                                                                                                                                                                                              tATCAGAAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCaa  <  2:378010/228‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                  
CATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCGCCCTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAA  >  NC_000913/2596388‑2596837

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