Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,683 T→C 14.3% intergenic (‑775/+102) rpsF ← / ← rlmB 30S ribosomal protein S6/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,545,6830TC14.3% 144.4 / 16.0 63intergenic (‑775/+102)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (24/30);  new base C (5/4);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

TGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545572‑5545830
                                                                                                               |                                                                                                                                                   
tGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTAc                                                                                                               <  2:1730664/150‑1 (MQ=32)
  ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                             >  1:292815/1‑150 (MQ=33)
   gAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                            <  2:1821271/150‑1 (MQ=34)
    aaCCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAggg                                                                                                           <  2:1357666/150‑1 (MQ=35)
         tAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACg                                                                                                      >  2:828246/1‑150 (MQ=39)
         tAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACg                                                                                                      <  1:828246/150‑1 (MQ=39)
           gaaaaagaaaaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTcctacagtggtcgcgt                                                                                                    >  1:1210988‑M2/92‑144 (MQ=255)
            cAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTNCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCaa                                                                                                   >  1:546609/1‑150 (MQ=255)
              aaGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATc                                                                                                 >  1:711260/1‑150 (MQ=255)
                   ttCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                            <  2:1730970/150‑1 (MQ=255)
                   ttCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                            <  1:177896/150‑1 (MQ=255)
                       gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  <  2:587805‑M2/29‑1 (MQ=255)
                       gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  >  1:587805‑M2/80‑108 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATctcctacgggcga                                                                                                         >  1:1679281‑M2/78‑125 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATctcctacgggcga                                                                                                         <  2:1679281‑M2/54‑7 (MQ=255)
                         ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTAcac                                                                                                             >  2:812894‑M2/78‑126 (MQ=255)
                         ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTAcac                                                                                                             <  1:812894‑M2/50‑2 (MQ=255)
                          gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAgg                                                                                                                                <  1:150851‑M2/31‑1 (MQ=255)
                          gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAgg                                                                                                                                >  2:150851‑M2/77‑107 (MQ=255)
                           cTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTAATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGc                                                                                    <  1:841225/150‑1 (MQ=255)
                            tGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAgag                                                                                   >  1:1078721/1‑148 (MQ=255)
                              aCCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCggg                                                                                 <  2:1153101/150‑1 (MQ=255)
                              aCCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCggg                                                                                 <  1:1323463/150‑1 (MQ=255)
                                cTACCCGTACCAGTCGGCGAATTCTTTTCCAGACCTCGATGTGAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTTCGACATCGCAGCCCGGAAGGGCTGGGGAAGCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGtt                                                                               <  2:546609/150‑1 (MQ=17)
                                  aCCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCa                                                                             <  2:722286/150‑1 (MQ=255)
                                                    aTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGaa                                                           <  2:292815/150‑1 (MQ=255)
                                                      gCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGt                                                                                          >  1:1820871/1‑117 (MQ=255)
                                                      gCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGt                                                                                          <  2:1820871/117‑1 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGATGCCg                                                              <  2:1501845/143‑1 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGATGCCg                                                              >  1:1501845/1‑143 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                >  1:495574/1‑141 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                <  2:495574/141‑1 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       >  1:586382/1‑150 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       >  2:1905201/1‑150 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       <  1:247231/150‑1 (MQ=255)
                                                        ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       >  1:932472/1‑150 (MQ=255)
                                                          ttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCt                                                     <  2:352547/150‑1 (MQ=255)
                                                          ttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCt                                                     <  2:586382/150‑1 (MQ=255)
                                                            cTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTga                                                   >  2:105125/1‑150 (MQ=255)
                                                              aGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgaga                                                 <  2:932472/150‑1 (MQ=255)
                                                                 tctcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTc                                              >  2:1407358/1‑150 (MQ=255)
                                                                        gTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCtt                                       <  1:105125/150‑1 (MQ=255)
                                                                        gTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCtt                                       <  2:1787043/150‑1 (MQ=255)
                                                                           aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAc                                    >  1:1440463/1‑150 (MQ=255)
                                                                           aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAc                                    >  1:879909/1‑150 (MQ=255)
                                                                           aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAc                                    >  2:1757958/1‑150 (MQ=255)
                                                                             tGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCt                                  >  1:1463306/1‑150 (MQ=255)
                                                                             ggCAGATCTGACGTCATGCGCATCCTCTTTGCTGTGCGACAGCGCAGCCCGGAAGGACTGGGTGATCTCGTACAGGGCAACGCAATCACCCTGTAGGAGCTGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCt                                  <  1:1622803/149‑1 (MQ=37)
                                                                              gCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTg                                 >  1:282935/1‑150 (MQ=255)
                                                                               cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                <  1:1905201/150‑1 (MQ=255)
                                                                               cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                >  2:1584717/1‑150 (MQ=255)
                                                                                       aCGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACg                        >  1:1810954/1‑150 (MQ=255)
                                                                                        cGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGc                       <  1:656629/150‑1 (MQ=255)
                                                                                            aTGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACgg                   >  1:85466/1‑150 (MQ=255)
                                                                                             tGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGc                  <  2:282935/150‑1 (MQ=255)
                                                                                                  aTCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGAGTGGGCAATCTCCGACAGGGCAAGGCAATCACCCCCTAGCACCGGATGAAGTGATTAAGATGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCAGACGCTCGGCTTCGa             >  1:1893460/1‑150 (MQ=16)
                                                                                                    ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgaga                                                 >  2:1259298/1‑112 (MQ=255)
                                                                                                    ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgaga                                                 <  1:1259298/112‑1 (MQ=255)
                                                                                                     cTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa          <  2:1810954/150‑1 (MQ=255)
                                                                                                      tGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAg         <  2:1440463/150‑1 (MQ=255)
                                                                                                       gTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGa        <  1:140621/150‑1 (MQ=255)
                                                                                                       gTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAACCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAGAGGCTGAGAGTCAGCGCTTGACCTGTCGCTGACGCACGGCTTCGAACAGa        <  2:711260/150‑1 (MQ=255)
                                                                                                        tttGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGAc       >  2:822848/1‑150 (MQ=255)
                                                                                                             tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcc  >  1:104757/1‑150 (MQ=255)
                                                                                                               tgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAc                                    >  1:846303/1‑114 (MQ=255)
                                                                                                               tgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAc                                    <  2:846303/114‑1 (MQ=255)
                                                                                                               |                                                                                                                                                   
TGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545572‑5545830

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