Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,114,917 T→C 7.4% intergenic (‑106/+73) rlmF ← / ← valS ribosomal RNA large subunit methyltransferase F/valine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,114,9170TC 7.4% 687.5 / 30.1 232intergenic (‑106/+73)rlmF/valSribosomal RNA large subunit methyltransferase F/valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (122/89);  new base C (8/9);  total (133/99)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

Only 100 of 235 total aligned reads displayed.

GGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTT  >  NC_002947/1114769‑1115066
                                                                                                                                                     |                                                                                                                                                     
ggtggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGtt                                                                                                                                                      <  1:1083335/150‑1 (MQ=255)
 gtggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTg                                                                                                                                                     >  2:5575373/1‑150 (MQ=255)
 gtggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTg                                                                                                                                                     <  2:5266378/150‑1 (MQ=255)
   ggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGgc                                                                                                                                                   <  1:6234169/150‑1 (MQ=255)
    gCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGgcg                                                                                                                                                  <  1:8455455/150‑1 (MQ=255)
       gTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCg                                                                                                                                               <  2:105041/150‑1 (MQ=255)
         gtccgaacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTc                                                                                                                                             >  2:5043050‑M1/98‑150 (MQ=255)
            gcgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctg                                                                                                                                          <  1:3627297‑M2/150‑6 (MQ=255)
            gcgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGCAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctg                                                                                                                                          <  2:4179861‑M2/150‑6 (MQ=255)
             gaacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCtgt                                                                                                                                         >  1:2919809‑M1/94‑150 (MQ=255)
             cgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgt                                                                                                                                         >  1:6689921‑M2/1‑144 (MQ=255)
             cgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgt                                                                                                                                         >  1:7819854‑M2/1‑144 (MQ=255)
              aacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCtgtg                                                                                                                                        >  2:8959290‑M1/93‑150 (MQ=255)
               acaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTgg                                                                                                                                       >  2:2884944‑M1/92‑150 (MQ=255)
               acaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTgg                                                                                                                                       >  1:629337‑M1/92‑150 (MQ=255)
               acaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTgg                                                                                                                                       >  1:8599186‑M1/92‑150 (MQ=255)
                 tGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtggga                                                                                                                                     >  2:7488614‑M2/1‑140 (MQ=255)
                   cAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagc                                                                                                                                   >  1:6975035‑M2/1‑138 (MQ=255)
                      ggtgggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggc                                                                                                                                <  1:7739148‑M2/150‑16 (MQ=255)
                         gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggcctt                                                                                                                             >  2:814651‑M2/1‑132 (MQ=255)
                         gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggcctt                                                                                                                             >  1:2699453‑M2/1‑132 (MQ=255)
                         gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTCGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggcctt                                                                                                                             >  1:6836159‑M2/1‑132 (MQ=255)
                           gtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgt                                                                                                                           <  1:1679481‑M2/150‑21 (MQ=255)
                           gtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggccttgt                                                                                                                           >  1:6636697‑M2/1‑130 (MQ=255)
                            tttGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtg                                                                                                                          <  2:3096936‑M2/150‑22 (MQ=255)
                                   ggggTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGaa                                                                                                                   >  1:635581/1‑150 (MQ=255)
                                    gggccctGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGACGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTTTGGGAGCGGCCTTGTGTCGCGaaa                                                                                                                  <  1:3879698/144‑1 (MQ=255)
                                    gggTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTCGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggccttgtgtcgcaaaa                                                                                                                  >  2:2031539‑M2/1‑121 (MQ=255)
                                       tCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGAGCCGTGCTGTGGGAGCGGCCTTGTGTAGCGAAAtgg                                                                                                               >  1:2246644/1‑150 (MQ=255)
                                            gCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTa                                                                                                          >  2:1648258/1‑150 (MQ=255)
                                              aGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAATGGCCGTACa                                                                                                        >  1:5637064/1‑150 (MQ=255)
                                                   ccAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGc                                                                                                   <  2:2805979/150‑1 (MQ=255)
                                                    cAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCGCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGcc                                                                                                  >  1:2655093/1‑150 (MQ=255)
                                                     aGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGccc                                                                                                 >  2:9085853/1‑150 (MQ=255)
                                                     aGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGccc                                                                                                 >  1:6058995/1‑150 (MQ=255)
                                                          tataCGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACg                                                                                            >  2:7033177/1‑150 (MQ=255)
                                                                   ggTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTa                                                                                   <  2:5056368/150‑1 (MQ=255)
                                                                     tGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGt                                                                                 <  2:2538259/150‑1 (MQ=255)
                                                                     gggCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGt                                                                                 <  2:7157233/149‑1 (MQ=255)
                                                                       gCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTTCAGCGGCCCCAACGGTCTCTCTAGTGa                                                                               >  2:8310902/1‑150 (MQ=255)
                                                                         ttGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAtt                                                                             >  2:1408319/1‑150 (MQ=255)
                                                                             aGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCGCAGGAGGGTGGCGGCTGGTAGGGTGGCGCCGTCCTGTGGGAGCGGCCGTGGGTCGCGAAAGGGCCGTCCAGCGGCCCCACCGGTCGCTCGAGGGagtcccg                                                                         >  2:9035697/1‑144 (MQ=37)
                                                                               gcaaaagggCTGCGCAGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAgcg                                                                       >  1:1868010/7‑150 (MQ=38)
                                                                                 ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCGGTGGGAGCGGCCGGGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCGAGTGATTACAGCGCt                                                                     >  2:290487/1‑150 (MQ=255)
                                                                                   tGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGc                                                                   >  2:8891729/1‑150 (MQ=255)
                                                                                   aagggCTGCGCAGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGc                                                                   >  1:4859970/3‑150 (MQ=255)
                                                                                         tgctgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCGGGGGGGAGCGGCCGTGTGTCGCGTGTGGGGCGTACGGCGGCCCCCGCGGTCTCTCTCGTGATTTCAGCGCGGCacgccg                                                             >  2:8626884/1‑145 (MQ=255)
                                                                                         tgctgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcg                                                             >  2:6042625/1‑150 (MQ=255)
                                                                                           ctgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcgcg                                                           >  2:7112421/1‑150 (MQ=255)
                                                                                            tgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCgcgc                                                          >  2:4418128/1‑150 (MQ=255)
                                                                                                 aGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCAtgct                                                     >  1:653457/1‑150 (MQ=255)
                                                                                                    ccTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGGTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACCGCGGCCCCAACGGTCTCTCTCGTGATTACAGAGCTGCAATACGCGCATgctgct                                                  >  1:7809057/1‑150 (MQ=38)
                                                                                                     ctttcGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                                 <  2:4700298‑M1/145‑1 (MQ=255)
                                                                                                     cTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                                 <  1:4341083/150‑1 (MQ=255)
                                                                                                      ttttGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCg                                                >  1:7170912/1‑150 (MQ=255)
                                                                                                         tGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTg                                             <  1:2754039/150‑1 (MQ=255)
                                                                                                          GCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCGTGAAGGTTGGCGCCGGCCTGTGGGGGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACCGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGa                                            <  1:3043121‑M1/150‑1 (MQ=255)
                                                                                                          GCGACACAAGGCCGCGCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCACTCCGCGCCTGCTGctcgggga                                            <  1:8705042‑M1/150‑4 (MQ=255)
                                                                                                           CGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTgac                                           <  2:7624254‑M1/150‑2 (MQ=255)
                                                                                                            GACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                          <  2:2406140‑M1/150‑1 (MQ=255)
                                                                                                            GACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                          <  1:5948924‑M1/150‑1 (MQ=255)
                                                                                                              CACAAGGCCGCTCCCACAGGAGGGTGGCGCCGTGAAGGTTGGCGCCGTCCGGTGGGCGCGGCCGGGGGTCGCGCAAGGGCCGTACAGCGGCCCCAACGGTCGCGCTAGTGATTACAGCGCGGCAATCCGCGCCTGCGGCTCGGTGaaggt                                        <  2:6276727‑M1/150‑3 (MQ=255)
                                                                                                               ACACGGCCGCTCCCACAGGTGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                       >  2:2655093‑M1/1‑150 (MQ=255)
                                                                                                               ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                       >  1:9085853‑M1/1‑150 (MQ=255)
                                                                                                               ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                       <  1:4974354‑M1/150‑1 (MQ=255)
                                                                                                               ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                       >  1:5277277‑M1/1‑150 (MQ=255)
                                                                                                               ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCGGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                       <  1:3568934‑M1/150‑1 (MQ=255)
                                                                                                                   GGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCa                                   >  2:3934555‑M1/1‑150 (MQ=255)
                                                                                                                     CCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGaa                                           >  1:8298931‑M1/1‑140 (MQ=255)
                                                                                                                         tCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgca                             <  1:7026483‑M2/150‑115 (MQ=255)
                                                                                                                           CCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGt                           <  2:5691588‑M1/150‑1 (MQ=255)
                                                                                                                            cacaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcgaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcaccc                          <  2:3007051‑M2/150‑118 (MQ=255)
                                                                                                                            cacaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcaccc                          <  1:4792292‑M2/150‑118 (MQ=255)
                                                                                                                            CACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTTTGGGATCGGCATTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                          >  2:3051892‑M1/1‑150 (MQ=255)
                                                                                                                            CACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                          <  1:1498618‑M1/150‑1 (MQ=255)
                                                                                                                            CACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCATGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                          >  2:7952164‑M1/1‑150 (MQ=255)
                                                                                                                             acaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccg                         <  1:4884560‑M2/150‑119 (MQ=255)
                                                                                                                             acaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccg                         <  2:6765289‑M2/150‑119 (MQ=255)
                                                                                                                             ACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTc                         <  2:1258444‑M1/150‑1 (MQ=255)
                                                                                                                             ACAGGAGGGTGGCGCCATGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTc                         <  1:2844381‑M1/150‑1 (MQ=255)
                                                                                                                              CAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAGTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCg                        >  2:6561689‑M1/1‑150 (MQ=255)
                                                                                                                               AGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCgg                       <  1:2324716‑M1/150‑1 (MQ=255)
                                                                                                                                GGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGa                      <  1:4383866‑M1/150‑1 (MQ=255)
                                                                                                                                 GAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGAc                     <  1:3940564‑M1/150‑1 (MQ=255)
                                                                                                                                   gggTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTc                   <  2:1679481/150‑1 (MQ=255)
                                                                                                                                    ggTGGCGCCTTGAAGGTTGGCGCCGtcctcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgca                  <  1:1376837‑M2/150‑126 (MQ=255)
                                                                                                                                      tGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGc                >  2:363225/1‑150 (MQ=255)
                                                                                                                                          gcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGc            >  1:520047/1‑150 (MQ=255)
                                                                                                                                          gcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGc            >  2:1868902/1‑150 (MQ=255)
                                                                                                                                          gcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGGGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGc            >  2:7699640/1‑150 (MQ=255)
                                                                                                                                            ccTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgcaaccgccacca           <  1:3275984‑M2/150‑133 (MQ=255)
                                                                                                                                             cTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCtt          >  1:1408330/1‑150 (MQ=255)
                                                                                                                                              ttGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTg         >  2:5631354/1‑150 (MQ=255)
                                                                                                                                              ttGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTg         >  1:554764/1‑150 (MQ=255)
                                                                                                                                              ttGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTg         >  2:4778657/1‑150 (MQ=255)
                                                                                                                                              ttGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTg         >  1:2109583/1‑150 (MQ=255)
                                                                                                                                                 aaGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGgcg      <  1:8891729/150‑1 (MQ=255)
                                                                                                                                                    gTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTa                                                                                                          >  2:736304/1‑47 (MQ=255)
                                                                                                                                                    gTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTa                                                                                                          <  1:736304/47‑1 (MQ=255)
                                                                                                                                                     ttGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGtt  <  1:8056541/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTT  >  NC_002947/1114769‑1115066

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