Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,546,502 T→C 8.9% V10V (GTA→GTG)  rlmB ← 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,546,5020TC 8.9% 191.0 / 7.0 67V10V (GTA→GTG) rlmB23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base C (3/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ATCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGCGC  >  NC_002947/5546353‑5546645
                                                                                                                                                     |                                                                                                                                               
atcTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCgtgt                                                                                                                                                 <  1:1423561/150‑1 (MQ=255)
       ggCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGt                                                                                                                                          >  1:795772/1‑150 (MQ=255)
       ggCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGt                                                                                                                                          >  1:1295875/1‑150 (MQ=255)
          gcTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCGTCCCACAGCCAGAGATGCCCATCCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                       <  2:1416069/150‑4 (MQ=255)
            tCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtt                                                                                                                                     >  1:162762/1‑150 (MQ=255)
            tCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgcagg                                                                                                                                     <  1:1021690/150‑6 (MQ=255)
               gCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAttttt                                                                                                                                  <  1:1045152/150‑1 (MQ=255)
                 CTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTc                                                                                                                                >  2:1352611‑M1/1‑150 (MQ=255)
                   GGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGATGCAAAAATGCCTGCACCGCGTGTACGCCGTAGATTTTTTCCa                                                                                                                              >  2:1155506‑M1/1‑150 (MQ=255)
                             cACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAt                                                                                                                    <  1:445725/150‑1 (MQ=255)
                                gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTcc                                                                                                                               >  2:281178/1‑136 (MQ=255)
                                gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTcc                                                                                                                               <  1:281178/136‑1 (MQ=255)
                                gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGaa                                                                                                                 >  2:1269846/1‑150 (MQ=255)
                                gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGaa                                                                                                                 >  2:740493/1‑150 (MQ=255)
                                   gTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTg                                                                                                              <  1:1269846/150‑1 (MQ=255)
                                   gTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTg                                                                                                              <  1:304123/150‑1 (MQ=255)
                                   gTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                       <  1:626293/125‑4 (MQ=255)
                                   gTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                       >  2:626293/1‑122 (MQ=255)
                                    ttttCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGc                                                                                                             <  1:118591/150‑1 (MQ=255)
                                        ccGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCtt                                                                                                         <  2:737835/150‑1 (MQ=255)
                                         cgcgGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTg                                                                                                        >  2:1244205/1‑150 (MQ=255)
                                          gcgGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGt                                                                                                       <  2:893548/150‑1 (MQ=255)
                                             gccagcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGa                                                                                                                  <  2:481607/136‑1 (MQ=255)
                                             gccagcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGa                                                                                                                  >  1:481607/1‑136 (MQ=255)
                                                agcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGc                                                                                                 >  2:499954/1‑150 (MQ=255)
                                                 gcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCg                                                                                                <  2:495401/150‑1 (MQ=255)
                                                  cgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGCTCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGt                                                                                               <  1:1021199/150‑1 (MQ=255)
                                                   gccagCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTg                                                                                              >  2:1227197‑M2/6‑150 (MQ=255)
                                                     cagCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                            <  2:531127‑M2/147‑1 (MQ=255)
                                                         aaGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                            <  1:621092/114‑1 (MQ=255)
                                                         aaGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                            >  2:621092/1‑114 (MQ=255)
                                                         aaGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAgcgc                                                                                        >  1:1363498/1‑150 (MQ=255)
                                                             tttCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGa                                                                                    >  1:316309/1‑150 (MQ=255)
                                                               tCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCt                                                                                                                           <  1:109284/109‑1 (MQ=255)
                                                               tCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCt                                                                                                                           >  2:109284/1‑109 (MQ=255)
                                                                      gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                           >  2:1160885/1‑150 (MQ=255)
                                                                      gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                           >  2:452021/1‑150 (MQ=255)
                                                                      gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCGGTCTCATGAATTGCGCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                           <  2:162762/150‑1 (MQ=255)
                                                                         gCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGcc                                                                        >  1:230760/1‑150 (MQ=255)
                                                                             gCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTa                                                                    <  1:1160885/150‑1 (MQ=255)
                                                                             gCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTa                                                                    <  1:452021/150‑1 (MQ=255)
                                                                             gCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTa                                                                    <  2:1363498/150‑1 (MQ=255)
                                                                               tGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACg                                                                  >  1:289732/1‑150 (MQ=255)
                                                                                    ccTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCt                                                             <  2:795772/150‑1 (MQ=255)
                                                                                      ttCCGACAGCCAGATCTGCTTGACGGGCTTTGGATGGTGCTGCAGAAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCAAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGGGCAGGACTTCGGCGGGCCTTTACGGTGGTTg                                                           >  2:548690/1‑150 (MQ=255)
                                                                                           aCAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCgg                                                      <  2:316309/150‑1 (MQ=255)
                                                                                                cAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGc                                                 <  1:34795/150‑1 (MQ=255)
                                                                                                    tCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGa                                             <  2:230760/150‑1 (MQ=255)
                                                                                                          tGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                       <  2:289732/150‑1 (MQ=255)
                                                                                                              ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGtt                                                    >  2:539796/1‑133 (MQ=255)
                                                                                                              ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGtt                                                    <  1:539796/133‑1 (MQ=255)
                                                                                                              ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                       >  1:570124/1‑146 (MQ=255)
                                                                                                              ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                       <  2:570124/146‑1 (MQ=255)
                                                                                                              ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCa                                   >  1:747193/1‑150 (MQ=255)
                                                                                                                  tttGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCtt                               <  1:1227197/150‑1 (MQ=255)
                                                                                                                     ggATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCAGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTtct                            <  1:1244205/150‑1 (MQ=255)
                                                                                                                      gATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTtctc                           <  2:1295875/150‑1 (MQ=255)
                                                                                                                              tgcAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgc                   <  2:26167/150‑1 (MQ=255)
                                                                                                                              tgcAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTCCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgc                   <  1:740493/150‑1 (MQ=255)
                                                                                                                               gcagcaATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgct                  <  1:728453/150‑1 (MQ=255)
                                                                                                                                    aaTGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAg             >  1:1006047/1‑150 (MQ=255)
                                                                                                                                         cTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGATGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActg        >  1:895013/1‑150 (MQ=255)
                                                                                                                                          tctaCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgc       >  1:357223/4‑150 (MQ=255)
                                                                                                                                           gCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgct      >  1:276131/1‑150 (MQ=255)
                                                                                                                                           gCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgct      <  2:276131/150‑1 (MQ=255)
                                                                                                                                            taCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgctg     >  2:123874/2‑150 (MQ=255)
                                                                                                                                            cACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgctg     >  2:23246/1‑150 (MQ=255)
                                                                                                                                               ggCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTgcgc  <  2:1324174/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGCGC  >  NC_002947/5546353‑5546645

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