Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 818,406 A→G 8.0% A90A (GCT→GCC)  PP_0704 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947818,4060AG 8.0% 172.2 / 4.9 62A90A (GCT→GCC) PP_0704ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (3/2);  total (33/29)
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.71e-01

GATGCTCAAACCGTCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAGGTATCCTGGG  >  NC_002947/818262‑818554
                                                                                                                                                |                                                                                                                                                     
gATGCTCAAACCGTCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGc                                                                                                                                                     >  4:177313/1‑147 (MQ=255)
    cTCAAACCGTCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑cc                                                                                                                                              >  8:119650/1‑149 (MQ=255)
             tCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTgg                                                                                                                                       >  2:341484/1‑147 (MQ=255)
             tCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTa                                                                                                                                     <  6:354342/149‑1 (MQ=255)
                aGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAAGCAATCCAGGAAGACAACCTTATCCACATCGCTTGGCAATTCTTCTTAGGGCACGTGAGCCAGG‑CCCACCTGGTAGGc                                                                                                                                  >  2:315425/1‑149 (MQ=255)
                 gcggggCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGTCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGc                                                                                                                                  <  4:229275/145‑1 (MQ=255)
                 gCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGcc                                                                                                                                 <  1:315425/149‑1 (MQ=255)
                       ccAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGcc                                                                                                                                                    >  5:78098/1‑125 (MQ=255)
                       ccAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGcc                                                                                                                                                    <  6:78098/125‑1 (MQ=255)
                                    ttGGCGGACGCGCGGTTTCAACTGCTGCAACGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTCCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCgcggcggc                                                                                                              >  5:4745/1‑149 (MQ=255)
                                        cGGAAGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCATCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGaa                                                                                                          >  2:96596/1‑149 (MQ=255)
                                            cgcgcgGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTgg                                                                                                      <  7:233868/149‑1 (MQ=255)
                                                cgGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGa                                                                                                  <  7:119650/149‑1 (MQ=255)
                                                   tttCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAgg                                                                                                                                                 >  1:94131/1‑100 (MQ=255)
                                                   tttCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAgg                                                                                                                                                 <  2:94131/100‑1 (MQ=255)
                                                       aaCTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGTGGCGGAACTGGTCGATCATCaa                                                                                           <  1:341484/149‑1 (MQ=255)
                                                          tgctgcAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCgcggcggcg                                                                                                             >  7:239849/1‑128 (MQ=255)
                                                          tgctgcAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCgcggcggcg                                                                                                             <  8:239849/128‑1 (MQ=255)
                                                             tgcAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACg                                                                                     >  6:59628/1‑149 (MQ=255)
                                                               cAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGTATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGGATCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCa                                                                                       <  1:96596/146‑1 (MQ=255)
                                                               cAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCaa                                                                                           >  7:19927/1‑141 (MQ=255)
                                                               cAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCaa                                                                                           <  8:19927/141‑1 (MQ=255)
                                                                       ggTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTCCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCGCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGGCACGGGCGAtgctg                                                                           <  6:4745/149‑1 (MQ=255)
                                                                        gTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGAtgctgc                                                                          <  1:133461/149‑1 (MQ=255)
                                                                         tCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATgctgct                                                                         <  3:177313/149‑1 (MQ=255)
                                                                           gAGGCGTTACAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGGCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATgctgct                                                                         >  1:163303/1‑147 (MQ=255)
                                                                            aGGCGCTCCAGGCTATCGACGATGACACGCGTATCCTCTTCGCTGGGCACTTCGGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGGAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAg                                                                      <  7:178740/149‑1 (MQ=255)
                                                                              gCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGGGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGa                                                                       <  5:48043/146‑1 (MQ=255)
                                                                              gCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGGGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGa                                                                       >  6:48043/1‑146 (MQ=255)
                                                                                gTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAg                                                                  >  7:72206/1‑149 (MQ=255)
                                                                                gTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGATCCAGCGCGCGTCGTCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAg                                                                  <  2:73027/149‑1 (MQ=255)
                                                                                    cAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGgcgc                                                              >  8:82467/1‑149 (MQ=255)
                                                                                      ggCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTc                                                                                                                            <  1:222215/85‑1 (MQ=255)
                                                                                      ggCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTc                                                                                                                            >  2:222215/1‑85 (MQ=255)
                                                                                       gCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGgcgcgc                                                            <  8:272088/148‑1 (MQ=255)
                                                                                       gCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGgcgcgc                                                            >  7:272088/1‑148 (MQ=255)
                                                                                         tATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGt                                                         >  7:200567/1‑149 (MQ=255)
                                                                                         tATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGt                                                         >  3:296566/1‑149 (MQ=255)
                                                                                             gAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCt                                                      <  4:296566/148‑1 (MQ=255)
                                                                                              aGGATGACAAGCTTATCCTCTTCGCTGGGCGCTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGat                                                                                                 >  5:354238/1‑104 (MQ=255)
                                                                                              aGGATGACAAGCTTATCCTCTTCGCTGGGCGCTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGat                                                                                                 <  6:354238/104‑1 (MQ=255)
                                                                                              aGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCttt                                                    >  7:173374/1‑149 (MQ=255)
                                                                                              aGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCttt                                                    <  6:66006/149‑1 (MQ=255)
                                                                                                gATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCttt                                                    >  6:137881/1‑147 (MQ=255)
                                                                                                 aTGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAg                                                 <  7:82467/149‑1 (MQ=255)
                                                                                                     cAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAgcggc                                             <  8:173374/149‑1 (MQ=255)
                                                                                                      aaGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGAtgc                                                                             <  7:149492/116‑1 (MQ=255)
                                                                                                      aaGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGAtgc                                                                             >  8:149492/1‑116 (MQ=255)
                                                                                                        gCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGa                                           >  8:154959/1‑148 (MQ=255)
                                                                                                              ccTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTgg                                    >  6:180749/1‑149 (MQ=255)
                                                                                                               ctctTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTggc                                   >  3:97962/1‑149 (MQ=255)
                                                                                                                tctTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTggcg                                  >  7:310847/1‑149 (MQ=255)
                                                                                                                              ttCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCaacaac                     >  1:181780/1‑148 (MQ=255)
                                                                                                                              ttCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAAcc                    >  2:94811/1‑149 (MQ=255)
                                                                                                                                cTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCaacaac                     <  2:181780/146‑1 (MQ=255)
                                                                                                                                cTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAAcc                    <  8:72206/147‑1 (MQ=255)
                                                                                                                                      ctaccAGGAGGGCCTGG‑CTCACCTGGTAGGGCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAg            >  1:30740/5‑149 (MQ=255)
                                                                                                                                      ctaccAGGAGGGCCTGG‑CTCACCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGcaa                         <  6:300671/132‑1 (MQ=255)
                                                                                                                                      ctaccAGGAGGGCCTGG‑CTCACCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGcaa                         >  5:300671/5‑136 (MQ=255)
                                                                                                                                       gggCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAg            <  5:180749/148‑1 (MQ=255)
                                                                                                                                           aGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAGGTATc       >  4:111785/1‑149 (MQ=255)
                                                                                                                                             gTGAGCCAGG‑CCCTACTGGTAGGCCCGCTCCAGCGCCAGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGAGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCTCGAAggaaatcg     >  1:184225/1‑143 (MQ=255)
                                                                                                                                                aGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAGGTATCCTggg  >  3:286316/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GATGCTCAAACCGTCTAGCTGGGCCAGCAGGAACGCTTGGCGGACGCGCGGTTTCAACTGCTGCAAGGCGCGGTCGAGGCGTTCCAGGCTATCGAGGATGACAAGCTTATCCTCTTCGCTGGGCACTTCTGCTTCGGGCAGGTGAGCCAGG‑CCCTCCTGGTAGGCCCGCTCCAGCGCGCGGCGGCGGAACTGGTCGATCATCAAGCCACGGGCGATGCTGCTGAGGTAGGCGCGCGGTTCTTTGAGCGGCGACATCTGGCGGGCGCGCAACAACCGCACGAAGGTATCCTGGG  >  NC_002947/818262‑818554

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