Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,683 T→C 8.7% 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,6830TC 8.7% 178.3 / 7.2 69intergenic (‑775/+102)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (3/3);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

TTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545571‑5545830
                                                                                                                 |                                                                                                                                                   
ttGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAAGGCTTTTGTAGATATCGATGACAGCCCGATATCGTCGACATACTAATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTa                                                                                                                >  1:896079/1‑150 (MQ=11)
   ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGCTCTTACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                             >  1:595710/1‑150 (MQ=18)
   ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                             <  2:1115718/150‑1 (MQ=33)
    gAACCATAGCAAAGTATTCGGCGTCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCATATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                            <  2:676159/150‑1 (MQ=21)
     aaCCATAGCAAAGTGTTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAggg                                                                                                           >  1:1499368/1‑150 (MQ=25)
        cATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGAAAGGGCTGGGTGATCTCCTACAGGGGaa                                                                                                        >  1:1449137/1‑150 (MQ=38)
          tAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACg                                                                                                      >  1:736871/1‑150 (MQ=39)
           aGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGAGGACCGc                                                                                                     >  2:1434156/1‑150 (MQ=25)
               aaGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATc                                                                                                 <  2:1812275/150‑1 (MQ=255)
                  tATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCTCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAATGGCTGGGTGATCTCCTACAGGGGAACGCAATCAcc                                                                                              <  1:151974/150‑1 (MQ=255)
                   aTTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAccc                                                                                             <  2:1901702/150‑1 (MQ=255)
                      cggcggCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGt                                                                                          <  1:1230445/150‑1 (MQ=255)
                       ggcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTa                                                                                         <  2:829043/150‑1 (MQ=255)
                       ggcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTa                                                                                         <  2:540290/150‑1 (MQ=255)
                       ggcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTa                                                                                         <  2:244189/150‑1 (MQ=255)
                       ggcggcTGACCTACCCGGATCACTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTTCAGGGGAACGCAATCACCCTGTa                                                                                         <  2:952355/150‑1 (MQ=255)
                        ggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  <  2:1588062‑M2/29‑1 (MQ=255)
                        ggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  >  1:1588062‑M2/81‑109 (MQ=255)
                           ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg                                                                                                                            <  2:102450‑M2/35‑1 (MQ=255)
                           ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg                                                                                                                            >  1:102450‑M2/78‑112 (MQ=255)
                           gcTGACCTACCCGGATCAGTCGGCGAATGCTTTAATTTGTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAg                                                                                     >  1:151574/1‑150 (MQ=17)
                            gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAgg                                                                                                                                >  1:1433341‑M2/77‑107 (MQ=255)
                            gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAgg                                                                                                                                <  2:1433341‑M2/31‑1 (MQ=255)
                              gACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGAAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCgg                                                                                  <  2:1449137/150‑1 (MQ=255)
                                             gTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                   >  1:1394438/1‑150 (MQ=255)
                                              tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAccc                                                                                             <  2:674295/123‑1 (MQ=255)
                                              tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAccc                                                                                             >  1:674295/1‑123 (MQ=255)
                                               cGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                 <  2:283852/150‑1 (MQ=255)
                                                ggCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                <  2:1724612/150‑1 (MQ=255)
                                                  ccaaTCCTTTTCTAGATCTCGATGTCATTGCAGATTTTACTTCATGCGCATCCTGTTTGCTGTGCGACACCTCAGCCCTGAATGGCTGGGTGTTCTCCTACAGTAGAACGCAATCACCCTTTAGGAGCTGGTTCACCCGCTAAGAGGcct                                                              <  2:896079/148‑2 (MQ=18)
                                                  cGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCg                                                              <  2:435051/150‑1 (MQ=255)
                                                       gCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGTTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACa                                                         <  2:635710/150‑1 (MQ=255)
                                                        cTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAg                                                        >  1:964687/1‑150 (MQ=255)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCGTCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                     >  2:1263266/1‑136 (MQ=255)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                     <  1:1263266/136‑1 (MQ=255)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGaa                                                           >  2:255456/1‑146 (MQ=255)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGaa                                                           <  1:255456/146‑1 (MQ=255)
                                                          tttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGc                                                      <  1:1458664/150‑1 (MQ=255)
                                                          tttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGc                                                      <  1:1457671/150‑1 (MQ=255)
                                                              tAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgag                                                  >  1:1309648/1‑150 (MQ=255)
                                                              tAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgag                                                  >  2:1146843/1‑150 (MQ=255)
                                                              tAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgag                                                  >  1:1486184/1‑150 (MQ=255)
                                                                    tcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAg                                            >  1:1736492/1‑150 (MQ=255)
                                                                    tcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAg                                            >  1:1693789/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                          >  2:115380/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                          >  1:145939/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                          >  1:252931/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                          >  2:1827635/1‑150 (MQ=255)
                                                                         gTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCtt                                       <  2:1309648/150‑1 (MQ=255)
                                                                              tGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                            <  1:1717274/76‑1 (MQ=39)
                                                                              tGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                            >  2:1717274/1‑76 (MQ=39)
                                                                                cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                <  1:628663/150‑1 (MQ=255)
                                                                                 aGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGcc                               >  1:1000191/1‑150 (MQ=255)
                                                                                   aTCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTCGACCTGCCGc                             >  2:1080751/1‑150 (MQ=255)
                                                                                   aTCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGAAAATCACTCACCCTGTAGGATTTTGTTCTGCCGGGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGc                             <  1:233436/150‑1 (MQ=21)
                                                                                           gTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCa                      <  2:1486184/150‑1 (MQ=255)
                                                                                           gTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGCTTCACCCGCGAAGAGGCCGGAACAGGCCGCGAGTCTGCTCTTGACCTGCCGCTgcctca                      <  2:850743/150‑6 (MQ=37)
                                                                                            tCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCAc                     <  1:334472/150‑1 (MQ=255)
                                                                                             cATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACg                    >  1:681071/1‑150 (MQ=255)
                                                                                               tGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGc                  >  2:195010/1‑150 (MQ=255)
                                                                                                  gcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTc               <  2:1693789/150‑1 (MQ=255)
                                                                                                    aTCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGACCGCAATCACGCGCTAGGAGGGGGGTTATCCAGAAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGa             >  1:611106/1‑150 (MQ=25)
                                                                                                     tCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGaa            <  1:92862/150‑1 (MQ=255)
                                                                                                      ccTGTTTGCTGTGCGACAGCGCCGATTCTAATGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAAc           <  1:1434156/150‑1 (MQ=38)
                                                                                                      ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATTACCCTGTAGGAGCGGGTTCCCCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAAc           <  2:522706/150‑1 (MQ=255)
                                                                                                         gTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGa        <  1:1532480/150‑1 (MQ=255)
                                                                                                          tttGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGAc       >  2:610160/1‑150 (MQ=255)
                                                                                                          tttGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGAc       <  1:115380/150‑1 (MQ=255)
                                                                                                            tGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGCGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACaa     >  1:1301219/1‑150 (MQ=255)
                                                                                                               tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcc  >  2:1925110/1‑150 (MQ=255)
                                                                                                                 |                                                                                                                                                   
TTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545571‑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: