Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,114,912 G→C 7.3% intergenic (‑101/+78) 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,9120GC 7.3% 765.6 / 37.9 233intergenic (‑101/+78)rlmF/valSribosomal RNA large subunit methyltransferase F/valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (124/92);  new base C (8/9);  total (132/101)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 234 total aligned reads displayed.

GACCCTGGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTG  >  NC_002947/1114763‑1115059
                                                                                                                                                      |                                                                                                                                                   
 aCCCTGGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGa                                                                                                                                                    <  1:2935323/150‑1 (MQ=255)
  cccTGGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGaa                                                                                                                                                   <  2:6851142/150‑1 (MQ=255)
   ccTGGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAg                                                                                                                                                  <  2:3358357/150‑1 (MQ=255)
      ggtggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGtt                                                                                                                                               <  1:1083335/150‑1 (MQ=255)
       gtggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTg                                                                                                                                              <  2:5266378/150‑1 (MQ=255)
         ggCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGgc                                                                                                                                            <  2:7937823/150‑1 (MQ=255)
             gTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCg                                                                                                                                        <  2:105041/150‑1 (MQ=255)
               gtccgaacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTc                                                                                                                                      >  1:6651850‑M1/98‑150 (MQ=255)
                  gcgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctg                                                                                                                                   <  2:8056541‑M2/150‑6 (MQ=255)
                   gaacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCtgt                                                                                                                                  >  1:5261074‑M1/94‑150 (MQ=255)
                   gaacaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCtgt                                                                                                                                  >  1:2919809‑M1/94‑150 (MQ=255)
                   cgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgt                                                                                                                                  >  2:5670891‑M2/1‑144 (MQ=255)
                   cgGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgt                                                                                                                                  >  2:6596804‑M2/1‑144 (MQ=255)
                    gggTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtg                                                                                                                                 >  2:3082558‑M2/1‑143 (MQ=255)
                    gggTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtg                                                                                                                                 >  1:6153359‑M2/1‑143 (MQ=255)
                     acaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTgg                                                                                                                                >  2:4786862‑M1/92‑150 (MQ=255)
                     acaggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTgg                                                                                                                                >  2:2885500‑M1/92‑150 (MQ=255)
                       tGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtggga                                                                                                                              >  1:6886651‑M2/1‑140 (MQ=255)
                       aggccctggccaacttcaccgagcagcatgcgcggattgcagcgctgtaatcactagagagaccgttggggccgctgtacggccctttcGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGa                                                                                                                              >  2:5465211‑M1/90‑150 (MQ=255)
                        gCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggag                                                                                                                             <  1:3861976‑M2/150‑12 (MQ=255)
                        gCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtccggtgggag                                                                                                                             <  2:8524221‑M2/150‑12 (MQ=255)
                            ggtgggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGAGACAAGGGCGCTGCCACAGGTGGGTGGCGCCGTGATGGTTGGGGGCGtcctgtgggagcggc                                                                                                                         <  1:3051892‑M2/150‑16 (MQ=255)
                               gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggcctt                                                                                                                      >  2:5053260‑M2/1‑132 (MQ=255)
                               gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggcctt                                                                                                                      >  2:9068870‑M2/1‑132 (MQ=255)
                               gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggcctt                                                                                                                      >  1:906124‑M2/1‑132 (MQ=255)
                               gggtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggcctt                                                                                                                      >  1:2699453‑M2/1‑132 (MQ=255)
                                 gtTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgt                                                                                                                    <  1:1679481‑M2/150‑21 (MQ=255)
                                      ttCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCATGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccgtgtgtcgc                                                                                                               <  2:8978710‑M2/150‑26 (MQ=255)
                                          gggccctGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGACGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTTTGGGAGCGGCCTTGTGTCGCGaaa                                                                                                           <  1:3879698/144‑1 (MQ=255)
                                          gggTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGaaa                                                                                                           <  2:2952123/150‑1 (MQ=255)
                                          gggTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTCGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggccttgtgtcgcaaaa                                                                                                           >  2:2031539‑M2/1‑121 (MQ=255)
                                           ggTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCtcctgtgggagcggccttgtgtcgcaaaag                                                                                                          >  2:5036788‑M2/1‑120 (MQ=255)
                                             tCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGAGCCGTGCTGTGGGAGCGGCCTTGTGTAGCGAAAtgg                                                                                                        >  1:2246644/1‑150 (MQ=255)
                                                  gCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTa                                                                                                   >  2:1648258/1‑150 (MQ=255)
                                                    aGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACa                                                                                                 >  1:5637254/1‑150 (MQ=255)
                                                    aGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACa                                                                                                 >  1:641237/1‑150 (MQ=255)
                                                         ccAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGc                                                                                            <  1:5303495/150‑1 (MQ=255)
                                                          cAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCGCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGcc                                                                                           >  1:2655093/1‑150 (MQ=255)
                                                                 ataCGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGCGCGGCCTTGTGTCGCGAACGGGCCGTACAGCGGCCCCAACgg                                                                                    >  1:6982062/1‑150 (MQ=255)
                                                                       gCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGtctctc                                                                              >  2:5599616/1‑150 (MQ=255)
                                                                          gTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTGGAAGGTTGGCGCCGTCCGGTGGGAGCGGCCTGGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAg                                                                           >  2:8735650/1‑150 (MQ=255)
                                                                           tGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGt                                                                          <  2:2538259/150‑1 (MQ=255)
                                                                           gggCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGt                                                                          <  2:7157233/149‑1 (MQ=255)
                                                                            ggCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTg                                                                         <  2:4401550/150‑1 (MQ=255)
                                                                            ggCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTGGCGAAAGGGCCGTACAGCGGCCCCCACGGTCTCTCGAGTg                                                                         >  1:8641460/1‑150 (MQ=255)
                                                                               ttGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAtt                                                                      >  1:6294865/1‑150 (MQ=255)
                                                                               ttGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAtt                                                                      >  2:1408319/1‑150 (MQ=255)
                                                                                     gcaaaagggCTGCGCAGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAgcg                                                                >  1:1868010/7‑150 (MQ=38)
                                                                                       ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                              >  1:6770213/1‑150 (MQ=255)
                                                                                                 ctgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcgcg                                                    >  2:7112421/1‑150 (MQ=255)
                                                                                                  tgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCgcgc                                                   >  2:4418128/1‑150 (MQ=255)
                                                                                                     gcAGCCCTTGTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGGTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCAtg                                                <  1:7987539/150‑1 (MQ=255)
                                                                                                       aGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCAtgct                                              >  2:1749176/1‑150 (MQ=255)
                                                                                                           ctttcGCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                          <  2:4700298‑M1/145‑1 (MQ=255)
                                                                                                           cTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                          >  2:2754380/1‑150 (MQ=255)
                                                                                                           cTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                          <  1:4341083/150‑1 (MQ=255)
                                                                                                           cTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTc                                          <  2:1261713/150‑1 (MQ=255)
                                                                                                            ttttGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCg                                         >  1:4861389/1‑150 (MQ=255)
                                                                                                               tGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTg                                      <  1:2754039/150‑1 (MQ=255)
                                                                                                                  GACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                   <  1:5948991‑M1/150‑1 (MQ=255)
                                                                                                                  GACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                   <  2:2406140‑M1/150‑1 (MQ=255)
                                                                                                                  GACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                   <  2:6463684‑M1/150‑1 (MQ=255)
                                                                                                                    CACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGtt                                 <  1:8221382‑M1/150‑1 (MQ=255)
                                                                                                                     ACACGGCCGCTCCCACAGGTGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                >  2:2655093‑M1/1‑150 (MQ=255)
                                                                                                                     ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                >  1:9085853‑M1/1‑150 (MQ=255)
                                                                                                                     ACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCGGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTg                                <  1:3568934‑M1/150‑1 (MQ=255)
                                                                                                                      CAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTgg                               <  1:8300112‑M1/150‑1 (MQ=255)
                                                                                                                       AAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                              <  2:5751292‑M1/150‑1 (MQ=255)
                                                                                                                       AAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                              <  1:8544600‑M1/150‑1 (MQ=255)
                                                                                                                       AAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                              <  1:7335586‑M1/150‑1 (MQ=255)
                                                                                                                         GGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCa                            >  2:3934555‑M1/1‑150 (MQ=255)
                                                                                                                           CCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGaa                                    <  2:8298931‑M1/140‑1 (MQ=255)
                                                                                                                            CGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAggg                         <  1:6587345‑M1/150‑1 (MQ=255)
                                                                                                                               tCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgca                      <  2:2817304‑M2/150‑115 (MQ=255)
                                                                                                                               tCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgca                      <  1:7026483‑M2/150‑115 (MQ=255)
                                                                                                                               tCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtccggtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcggatagggctggacactgccatgaccccgaacaaacccaccctgca                      <  2:7398708‑M2/150‑115 (MQ=255)
                                                                                                                                 ccACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacc                    <  1:5568239‑M2/150‑117 (MQ=255)
                                                                                                                                  cacaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcaccc                   <  1:6250829‑M2/150‑118 (MQ=255)
                                                                                                                                  cacaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcaccc                   <  1:4792292‑M2/150‑118 (MQ=255)
                                                                                                                                  cacaGGAGGGTGGCGCCGTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcgaaagggctgcgcagcagccccagggtcttcaagccaccgcgcgtctagggctggacactgccatgaccccgaacaagcccaccctgcaccc                   <  2:3007073‑M2/150‑118 (MQ=255)
                                                                                                                                  CACCGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGTGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                   >  1:4965246‑M1/1‑150 (MQ=255)
                                                                                                                                  CACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                   <  1:1498618‑M1/150‑1 (MQ=255)
                                                                                                                                   acaGGAGGGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccg                  <  1:5076894‑M2/150‑119 (MQ=255)
                                                                                                                                   ACAGGAGGGTGGCGCCATGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTc                  <  1:2844381‑M1/150‑1 (MQ=255)
                                                                                                                                    CAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCg                 >  2:7636921‑M1/1‑150 (MQ=255)
                                                                                                                                     aggaggGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcg                <  1:7471371‑M2/150‑121 (MQ=255)
                                                                                                                                     AGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCgg                <  1:2324716‑M1/150‑1 (MQ=255)
                                                                                                                                      GGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGa               <  1:4383866‑M1/150‑1 (MQ=255)
                                                                                                                                       gaggGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcgaaagggccgtacagcggccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgca              <  1:776428‑M2/150‑123 (MQ=255)
                                                                                                                                       gaggGTGGCGCCTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgca              <  2:6213021‑M2/150‑123 (MQ=255)
                                                                                                                                       GAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGAc              <  1:3940564‑M1/150‑1 (MQ=255)
                                                                                                                                          ggTGGCGCCTTGAAGGTTGGCGCCGtcctcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgca           <  1:1376837‑M2/150‑126 (MQ=255)
                                                                                                                                          gacggCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCg           >  2:7616215/4‑150 (MQ=255)
                                                                                                                                            tGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGc         >  2:363225/1‑150 (MQ=255)
                                                                                                                                              gcgcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCa       >  2:8892698/1‑150 (MQ=255)
                                                                                                                                                gcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGc     >  2:1868902/1‑150 (MQ=255)
                                                                                                                                                  ccTTGAAGGTTGGCGCCGtcctgtgggagcggccttgtgtcgcaaaagggctgcgcagcagccccagggtcttcaagccaccgctcgtatagggctggacactgccatgaccccgaacaaacccaccctgcacccgcgcaaccgccacca    <  1:3275984‑M2/150‑133 (MQ=255)
                                                                                                                                                  ccTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCt    >  2:706267/1‑150 (MQ=255)
                                                                                                                                                   cTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCtt   >  1:1408330/1‑150 (MQ=255)
                                                                                                                                                    ttGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTg  >  1:2109583/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GACCCTGGTGGCGGTTGCGCGGGTGCAGGGTGGGTTTGTTCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTG  >  NC_002947/1114763‑1115059

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