Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,278,042 Δ1 bp 5.3% coding (192/1125 nt) rlmG → ribosomal RNA large subunit methyltransferase G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,278,0420A. 5.3% 1046.9 / 45.9 302coding (192/1125 nt)rlmGribosomal RNA large subunit methyltransferase G
Reads supporting (aligned to +/- strand):  ref base A (163/123);  new base . (5/11);  total (168/134)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

CGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGGTGCTGG  >  NC_002947/5277893‑5278190
                                                                                                                                                      |                                                                                                                                                    
cGCCAGCGGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCa                                                                                                                                                      <  2:3343610/150‑1 (MQ=255)
cGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCa                                                                                                                                                      <  1:6523188/150‑1 (MQ=255)
cGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCa                                                                                                                                                      <  2:1936026/150‑1 (MQ=255)
cGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCa                                                                                                                                                      <  2:4817503/150‑1 (MQ=255)
cGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCa                                                                                                                                                      <  1:2697715/150‑1 (MQ=255)
 gCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAg                                                                                                                                                     <  1:1446719/150‑1 (MQ=255)
 gCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAg                                                                                                                                                     <  1:2429010/150‑1 (MQ=255)
 gCCAGCCGGAGCAGGCCAACGACCCGGTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAg                                                                                                                                                     <  1:3531207/150‑1 (MQ=255)
     gCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCa                                                                                                                                                 >  1:684872/1‑150 (MQ=255)
     gCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCa                                                                                                                                                 <  2:3661393/150‑1 (MQ=255)
       cGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGCGGCCAgg                                                                                                                                               <  1:523577/150‑2 (MQ=255)
       cGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGc                                                                                                                                               >  2:2691463/1‑150 (MQ=255)
           gCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGCGGCCAgggcac                                                                                                                                           <  2:5338959/150‑6 (MQ=255)
            cAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGTGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAg                                                                                                                                          >  2:3901858/1‑150 (MQ=255)
            cAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAg                                                                                                                                          >  1:2573427/1‑150 (MQ=255)
            cAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAg                                                                                                                                          <  1:2845072/150‑1 (MQ=255)
              ggCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAggt                                                                                                                                        >  2:4294549/1‑150 (MQ=255)
                 cAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAggtggt                                                                                                                                     <  2:6607397/150‑1 (MQ=255)
                  caCGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCAGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTc                                                                                                                                    <  2:835194/149‑1 (MQ=255)
                  aaCGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGAGGTc                                                                                                                                    >  1:5388510/1‑150 (MQ=255)
                        ccGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCagcagc                                                                                                                              >  2:1269182/1‑150 (MQ=255)
                          gTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCgg                                                                                                                            >  1:3468505/1‑150 (MQ=255)
                              cAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGAc                                                                                                                        <  1:3223308/150‑1 (MQ=255)
                              cAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGAc                                                                                                                        <  2:6551073/150‑1 (MQ=255)
                               aGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACt                                                                                                                       >  2:5080336/1‑150 (MQ=255)
                                ggCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCGATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTc                                                                                                                      <  2:5976240/150‑1 (MQ=255)
                                ggCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTc                                                                                                                      <  1:704807/150‑1 (MQ=255)
                                 gCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCg                                                                                                                     >  1:442894/1‑150 (MQ=255)
                                  ccTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGc                                                                                                                    >  1:5679818/1‑150 (MQ=255)
                                    ttCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCAc                                                                                                                  <  1:3589001/150‑1 (MQ=255)
                                     tCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCAcc                                                                                                                 >  2:5249770/1‑150 (MQ=255)
                                        aTGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTgg                                                                                                              <  1:4041542/150‑1 (MQ=255)
                                        aTGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTgg                                                                                                              <  1:3851993/150‑1 (MQ=255)
                                            cgccgACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCa                                                                                                          >  2:3438353/1‑150 (MQ=255)
                                            cgccgACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCa                                                                                                          >  2:5599794/1‑150 (MQ=255)
                                             gccgACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCAc                                                                                                         <  1:6076258/150‑1 (MQ=255)
                                                   gaGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGcc                                                                                                   >  1:3749050/1‑150 (MQ=255)
                                                     gTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCt                                                                                                 <  1:4333750/150‑1 (MQ=255)
                                                     gTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCt                                                                                                 >  1:1042287/1‑150 (MQ=255)
                                                       aCCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAACTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTgg                                                                                               >  1:3253900/1‑150 (MQ=255)
                                                         cTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGaa                                                                                             <  2:1801626/150‑1 (MQ=255)
                                                          tGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGaaa                                                                                            <  2:6250278/150‑1 (MQ=255)
                                                             ttGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGa                                                                                         >  2:4620307/1‑150 (MQ=255)
                                                              tGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGaa                                                                                        <  2:6764988/150‑1 (MQ=255)
                                                                cgcAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAAcc                                                                                      <  2:2059101/150‑1 (MQ=255)
                                                                cgcAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAAcc                                                                                      >  1:6797960/1‑150 (MQ=255)
                                                                 gcAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCt                                                                                     >  1:592974/1‑150 (MQ=255)
                                                                 gcAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCt                                                                                     >  2:305273/1‑150 (MQ=255)
                                                                        cacGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGgcgcgc                                                                              >  1:1136105/1‑150 (MQ=255)
                                                                           gATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGaa                                                                                             >  2:4076454/1‑132 (MQ=255)
                                                                           gATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAAt                                                                           >  1:2021776/1‑150 (MQ=255)
                                                                           gATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAAt                                                                           >  1:5467414/1‑150 (MQ=255)
                                                                               aGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGcc                                                                       >  1:4823225/1‑150 (MQ=255)
                                                                                ggCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCt                                                                      >  2:2209211/1‑150 (MQ=255)
                                                                                ggCACCCGACGCCAAGTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCt                                                                      <  1:4715885/150‑1 (MQ=255)
                                                                                 gCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTg                                                                     >  2:5805286/1‑150 (MQ=255)
                                                                                 gCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTg                                                                     >  2:694532/1‑150 (MQ=255)
                                                                                      cGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGtt                                                                <  2:2764839/150‑1 (MQ=255)
                                                                                             aaTTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAgt                                                         <  1:515861/150‑1 (MQ=255)
                                                                                             aaTTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGACCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAgt                                                         <  2:3602528/150‑1 (MQ=255)
                                                                                                tGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGtgtg                                                      >  2:2967793/1‑150 (MQ=255)
                                                                                                      gTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTc                                                <  1:1255322/150‑1 (MQ=255)
                                                                                                      gTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTc                                                <  1:4078395/150‑1 (MQ=255)
                                                                                                      gTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTc                                                <  2:2870742/150‑1 (MQ=255)
                                                                                                          tGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTAc                                            >  1:2367860/1‑150 (MQ=255)
                                                                                                            gTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCg                                          <  1:4175012/150‑1 (MQ=255)
                                                                                                            gTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCg                                          <  1:18884/150‑1 (MQ=255)
                                                                                                            gTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCg                                          <  1:1702739/150‑1 (MQ=255)
                                                                                                              gCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGc                                        <  1:1977700/150‑1 (MQ=255)
                                                                                                                tCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCa                                      <  1:6462108/150‑1 (MQ=255)
                                                                                                                tCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCa                                      >  1:2894643/1‑150 (MQ=255)
                                                                                                                 cAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAg                                     >  1:1630110/1‑150 (MQ=255)
                                                                                                                  aaTGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAgg                                                                                                                                         <  2:6122024/49‑1 (MQ=255)
                                                                                                                  aaTGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGc                                    <  1:308103/150‑1 (MQ=255)
                                                                                                                  aaTGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCg                                   >  1:6020641/1‑150 (MQ=255)
                                                                                                                    tGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGa                                  >  1:6214808/1‑150 (MQ=255)
                                                                                                                     gACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGaa                                 <  1:978069/150‑1 (MQ=255)
                                                                                                                      aCAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAac                                >  1:6322403/1‑150 (MQ=255)
                                                                                                                       cAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAaca                               <  1:2323651/150‑1 (MQ=255)
                                                                                                                       cAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAaca                               <  2:544162/150‑1 (MQ=255)
                                                                                                                        aGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAacac                              >  2:4362805/1‑150 (MQ=255)
                                                                                                                          cTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTg                            >  2:184313/1‑150 (MQ=255)
                                                                                                                          cTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTg                            <  1:900478/150‑1 (MQ=255)
                                                                                                                          cTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTg                            >  2:6268558/1‑150 (MQ=255)
                                                                                                                              ggTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCaa                        <  2:1645795/150‑1 (MQ=255)
                                                                                                                              ggTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCaa                        <  2:6208732/150‑1 (MQ=255)
                                                                                                                               gTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCGATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAg                       >  2:4030588/1‑150 (MQ=255)
                                                                                                                               gTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAg                       >  1:422854/1‑150 (MQ=255)
                                                                                                                                     tGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGt                 >  2:1479571/1‑150 (MQ=255)
                                                                                                                                      ggcctgCAAGGCTTGCGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACa                                                           >  2:4536441/6‑107 (MQ=255)
                                                                                                                                      ggcctgCAAGGCTTGCGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACa                                                           <  1:4536409/102‑1 (MQ=255)
                                                                                                                                      ggcctgCAAGGCTTGCGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACa                                                           >  2:1092694/6‑107 (MQ=255)
                                                                                                                                      ggcctgCAAGGCTTGCGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCAGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGtt                >  1:3713706/6‑150 (MQ=255)
                                                                                                                                      ggCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGtt                >  2:3244428/1‑150 (MQ=255)
                                                                                                                                         cGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGa             >  2:4732372/1‑150 (MQ=255)
                                                                                                                                              aGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCgg         <  2:4158727/150‑1 (MQ=255)
                                                                                                                                               gCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCggg        <  1:6636898/150‑1 (MQ=255)
                                                                                                                                                    gCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGGTGCTg   <  2:3819114/150‑1 (MQ=255)
                                                                                                                                                    gCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGGTGCTg   <  2:399632/150‑1 (MQ=255)
                                                                                                                                                     cAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGGTGCTgg  >  1:2130113/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGCCGACGAGTACCTGCTTGCGCAGCTACACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCGCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGGCGCGCAATGGCCTGCCGTTCGACAGTGTGCCGTTCGTACCGGCCAGCGAACACTGGCAAGGCCCGTTCGACCGGGTGCTGG  >  NC_002947/5277893‑5278190

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