Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,683 T→C 13.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,6830TC13.3% 299.9 / 31.1 120intergenic (‑775/+102)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (53/51);  new base C (8/8);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 123 total aligned reads displayed.

TTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACCCC  >  NC_002947/5545571‑5545832
                                                                                                                |                                                                                                                                                     
   ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                               >  1:1027229/1‑150 (MQ=33)
   ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                               >  2:1290202/1‑150 (MQ=33)
        cATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCTAAAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGaa                                                                                                          >  1:1450425/1‑150 (MQ=255)
         aTAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAAc                                                                                                         <  1:2632265/150‑1 (MQ=38)
             ccaaGTATTCTGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCaa                                                                                                     <  1:2821327/148‑1 (MQ=255)
             cAAAGTATTCGGCGGCTGACCTACCCGGATCGGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCaa                                                                                                     >  2:1893993/1‑150 (MQ=255)
              aaaGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAAt                                                                                                    <  1:2143704/150‑1 (MQ=255)
              aaaGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAAt                                                                                                    >  1:2943023/1‑150 (MQ=255)
                 gTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAc                                                                                                 <  2:1027229/150‑1 (MQ=255)
                       ggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTgg                                                                                                                             <  1:1221705‑M2/36‑1 (MQ=255)
                       ggcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTa                                                                                           <  2:2503277/150‑1 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           <  2:2074665‑M2/38‑1 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           >  1:2074665‑M2/80‑117 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCtacgg                                                                                                               <  2:1858300‑M2/50‑3 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCtacgg                                                                                                               >  1:1858300‑M2/80‑127 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTacggggat                                                                                                           <  2:1582282‑M2/54‑3 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                                                                                                                >  2:2818628‑M2/80‑128 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                                                                                                                <  1:2818628‑M2/49‑1 (MQ=255)
                        gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACGGGgat                                                                                                          >  1:1582282‑M2/80‑132 (MQ=255)
                          ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCa                                                                                                                                            <  2:2043326‑M2/21‑1 (MQ=255)
                          ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCa                                                                                                                                            >  1:2043326‑M2/78‑98 (MQ=255)
                           gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAtc                                                                                                                       <  1:2552080‑M2/42‑1 (MQ=255)
                           gttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAtc                                                                                                                       >  2:2552080‑M2/77‑118 (MQ=255)
                           gcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAg                                                                                       >  1:1107021/1‑150 (MQ=255)
                            ttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAt                                                                                                                        <  2:2780661‑M2/41‑1 (MQ=255)
                            ttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAt                                                                                                                        >  1:2780661‑M2/76‑116 (MQ=255)
                            ttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGcgcg                                                                                                                                                    <  1:1489535‑M2/13‑1 (MQ=255)
                              gACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              <  1:1758764/124‑1 (MQ=35)
                              gACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              >  2:1758764/1‑124 (MQ=35)
                              gACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCgg                                                                                    <  2:2026050/150‑1 (MQ=255)
                               aCCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCggg                                                                                   <  1:1831682/150‑1 (MQ=255)
                                ccTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTTATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTNATCTCCTACAGGNGAACGCAATCACCCTGTAGGAGCGGGt                                                                                  <  1:3041049/150‑1 (MQ=255)
                                ccTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGt                                                                                  >  2:1915441/1‑150 (MQ=255)
                                 cTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCt                                                                                                                   <  1:2248070/116‑1 (MQ=11)
                                 cTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCt                                                                                                                   >  2:2248070/1‑116 (MQ=11)
                                 cTACCCGGATCAGTCGGCGAATGCTTTNCTAGATCTCGATNTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGtt                                                                                 >  1:2659684/1‑150 (MQ=255)
                                   aCCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCa                                                                               >  1:1175691/1‑150 (MQ=255)
                                    cccGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCAc                                                                              >  2:3005666/1‑150 (MQ=255)
                                     ccGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCAcc                                                                             >  1:1508799/1‑150 (MQ=255)
                                          tCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              <  2:1138549/128‑1 (MQ=255)
                                          tCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              >  1:1138549/1‑128 (MQ=255)
                                             gTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                     <  2:1474336/150‑1 (MQ=255)
                                              tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTc                                                                                >  1:471781/1‑138 (MQ=255)
                                              tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTc                                                                                <  2:471781/138‑1 (MQ=255)
                                                ggCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGGCGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                  <  1:1937572/150‑1 (MQ=255)
                                                 gCGAATGCTTTTCTAGATCTAGATGTCAGTGAAGATCTGACGTAATGCACATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGcc                                                                 >  2:350200/1‑150 (MQ=255)
                                                  cGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCg                                                                <  2:1175691/150‑1 (MQ=255)
                                                   gAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCTACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCgg                                                               <  2:1450425/150‑1 (MQ=255)
                                                    aaTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              <  1:2641540/102‑1 (MQ=35)
                                                    aaTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              >  2:2641540/1‑102 (MQ=35)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                                                                                                                >  1:175953/1‑95 (MQ=34)
                                                         ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                                                                                                                <  2:175953/95‑1 (MQ=34)
                                                         gtttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGTGTTCACCCGCGAAGAGGCCGGAACAgg                                                         >  2:2107471/2‑150 (MQ=255)
                                                                 aTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCTGGAACAGGCTga                                                     <  2:755856/146‑1 (MQ=255)
                                                                 aTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgaga                                                   >  2:729335/1‑148 (MQ=255)
                                                                 aTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTgaga                                                   <  1:729335/148‑1 (MQ=255)
                                                                 aTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTga                                                     >  1:755856/1‑146 (MQ=255)
                                                                   ctcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCa                                               >  1:1902888/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAAt                                                                                                    >  2:1210008/1‑94 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAAt                                                                                                    <  1:1210008/94‑1 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                            >  2:1043835/1‑150 (MQ=255)
                                                                      gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCATCAGCGCCTCCCGGGAGGGCTGGGATATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAgcg                                            >  1:57948/1‑150 (MQ=37)
                                                                         gTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCtt                                         <  2:1670856/150‑1 (MQ=255)
                                                                           cAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGa                                       <  2:2659684/150‑1 (MQ=255)
                                                                           cAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGa                                       <  1:3005666/150‑1 (MQ=255)
                                                                            aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAc                                                            <  2:683884/128‑1 (MQ=255)
                                                                            aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAc                                                            >  1:683884/1‑128 (MQ=255)
                                                                            aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAc                                                            >  1:683169/1‑128 (MQ=255)
                                                                             gTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTg                                                                                             <  1:703660/94‑1 (MQ=255)
                                                                             gTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTg                                                                                             >  2:703660/1‑94 (MQ=255)
                                                                             gTGCAGATCTGACCTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAcc                                     >  1:1352846/1‑150 (MQ=255)
                                                                                cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                  <  2:1902888/150‑1 (MQ=255)
                                                                                cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                  >  1:1534075/1‑150 (MQ=255)
                                                                                 aGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGcc                                 <  2:2896896/150‑1 (MQ=255)
                                                                                  gATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCg                                >  1:1664421/1‑150 (MQ=255)
                                                                                    tCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGa                                       >  2:1782552/1‑141 (MQ=255)
                                                                                    tCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCt                              >  2:2685945/1‑150 (MQ=255)
                                                                                    tCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCt                              >  2:1717578/1‑150 (MQ=255)
                                                                                      tGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGa                            >  2:2175754/1‑150 (MQ=255)
                                                                                        aCGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACg                          <  1:2175754/150‑1 (MQ=255)
                                                                                        aCGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACg                          >  2:1076144/1‑150 (MQ=255)
                                                                                            cATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGcc                                 >  2:17438/1‑139 (MQ=255)
                                                                                            cATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGcc                                 <  1:17438/139‑1 (MQ=255)
                                                                                             aTGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACgg                     >  2:22215/1‑150 (MQ=255)
                                                                                             aTGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCCTGACCTGCCGCTGACGCACgg                     >  1:263343/1‑150 (MQ=255)
                                                                                               gcgcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCt                   >  1:2100789/1‑150 (MQ=255)
                                                                                                 gcATCCTGTTTTCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCaa                                                                                                     <  1:1507135/66‑1 (MQ=255)
                                                                                                 gcATCCTGTTTTCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCaa                                                                                                     >  2:1507135/1‑66 (MQ=255)
                                                                                                 gcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTc                 <  1:1043835/150‑1 (MQ=255)
                                                                                                 gcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTc                 <  1:2685945/150‑1 (MQ=255)
                                                                                                    tCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCcgc                                                                          <  1:2268450/90‑1 (MQ=255)
                                                                                                    tCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGaa              <  2:57948/150‑1 (MQ=255)
                                                                                                     ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              >  2:2249138/1‑69 (MQ=255)
                                                                                                     ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              <  1:2249138/69‑1 (MQ=255)
                                                                                                     ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAAc             <  2:1508799/150‑1 (MQ=255)
                                                                                                     ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAAc             <  2:1352846/150‑1 (MQ=255)
                                                                                                      cTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa            <  1:1815795/150‑1 (MQ=255)
                                                                                                      cTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa            >  2:2086426/1‑150 (MQ=255)
                                                                                                      cTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa            >  2:1137746/1‑150 (MQ=255)
                                                                                                                tgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGTTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcccc  >  2:2567431/1‑150 (MQ=255)
                                                                                                                |                                                                                                                                                     
TTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACCCC  >  NC_002947/5545571‑5545832

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