Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,067 G→A 6.5% intergenic (‑111/+17) PP_0386 ← / ← rpoD sensory box protein/RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,0670GA 6.5% 177.0 / 2.5 61intergenic (‑111/+17)PP_0386/rpoDsensory box protein/RNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base A (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

ATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGG  >  NC_002947/470919‑471210
                                                                                                                                                    |                                                                                                                                               
aTCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTg                                                                                                                                                 <  6:55596/149‑1 (MQ=255)
      aTCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgg                                                                                                                                           <  7:227740/149‑1 (MQ=255)
         aGCAACAGAATGGCCGGCTTCAGGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCgggg                                                                                                                                         <  4:34936/148‑1 (MQ=255)
          tcAACAGAATGGCAGGCTTCAAGGACATGTCAGCCTCGCACTCAGGTTGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCCGCCTGGGCCGGGGtt                                                                                                                                       <  8:285708/148‑1 (MQ=255)
            aaCAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGTAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGtt                                                                                                                                       >  7:285708/1‑147 (MQ=255)
               aGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGt                                                                                                                                   >  6:86911/1‑148 (MQ=255)
                   tGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGtcatca                                                                                                                              <  5:46064/149‑1 (MQ=255)
                               aGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGagg                                                                                                                    >  3:134145/1‑147 (MQ=255)
                                       gCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAgcgc                                                                                                           <  3:40594/148‑1 (MQ=255)
                                       gCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAgcg                                                                                                            >  2:61736/1‑147 (MQ=255)
                                        cAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAg                                                                                                         <  1:49007/149‑1 (MQ=255)
                                             tCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCt                                                                                                    >  6:291230/1‑149 (MQ=255)
                                                 aCTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCt                                                                                                >  2:226192/1‑149 (MQ=255)
                                                  cTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTg                                                                                               <  2:22041/149‑1 (MQ=255)
                                                     aGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCACGTGTGCAACGGAAGGCGCGGTTTATGTCACCGGCCACTAAAGAAAAACCCCGGACTGGGCCGGGCTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                              >  8:179651/1‑147 (MQ=255)
                                                        gggtCTGCAGTCGCGGGGAGACGGGCATTGCCGCAGGGGGGCAACGGCAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACGCGTCGAGGAAGGAGCGCAGGTGCGCGCTGCGAGTc                                                                                         <  5:22320/148‑1 (MQ=255)
                                                          gtCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                       >  3:37634/1‑149 (MQ=255)
                                                              gAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaagg                                                                                                                >  6:64988/1‑120 (MQ=255)
                                                              gAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGaggaagg                                                                                                                <  5:64988/120‑1 (MQ=255)
                                                                  tCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGgcg                                                                                 >  3:256273/1‑147 (MQ=255)
                                                                    gAGTGTAGACGGGCATTGCCGCAGGTGTGAAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGc                                                                             >  5:190106/1‑149 (MQ=255)
                                                                         tAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTgc                                                                         <  4:134145/148‑1 (MQ=255)
                                                                          aGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTgg                                                                                   >  8:77714/1‑137 (MQ=255)
                                                                          aGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTgg                                                                                   <  7:77714/137‑1 (MQ=255)
                                                                                 cATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCct                                                                >  8:42836/1‑149 (MQ=255)
                                                                                  aTTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCctt                                                               >  7:72142/1‑149 (MQ=255)
                                                                                            gtgtGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATc                                                     >  3:211530/1‑149 (MQ=255)
                                                                                               tgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCctt                                                               <  7:288345/136‑1 (MQ=255)
                                                                                               tgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCctt                                                               >  8:288345/1‑136 (MQ=255)
                                                                                               tgCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGa                                                  <  1:61736/149‑1 (MQ=255)
                                                                                                gCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                       <  5:217913/111‑1 (MQ=255)
                                                                                                gCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCgg                                                                                       >  6:217913/1‑111 (MQ=255)
                                                                                                      gAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATcc                                           <  4:256273/149‑1 (MQ=255)
                                                                                                          gCGGTCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTc                                       <  7:42836/149‑1 (MQ=255)
                                                                                                          gCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGCGGAAGGAGAGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGAGCAACGCCTTCGCTTCGATCTGACGGATCCGCtg                                       >  3:228599/1‑148 (MQ=255)
                                                                                                           cGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                              <  1:143407/93‑1 (MQ=255)
                                                                                                           cGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGc                                                                                              >  2:143407/1‑93 (MQ=255)
                                                                                                           cGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCa                                      <  5:86911/149‑1 (MQ=255)
                                                                                                            gggCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCAc                                     >  1:198413/1‑149 (MQ=255)
                                                                                                              gCTTTATGTCACGCGCCACTAAAGAAAAACCCCGGCCTGGGCCGCGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATACGCTCAcgc                                   >  6:238861/1‑149 (MQ=255)
                                                                                                                tttATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGt                                 >  5:207834/1‑149 (MQ=255)
                                                                                                                    tGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGTCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTAc                              <  6:190106/148‑1 (MQ=255)
                                                                                                                        aCGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCa                          >  6:292088/1‑148 (MQ=255)
                                                                                                                         cGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                        <  8:72142/149‑1 (MQ=255)
                                                                                                                          gggCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAAc                       >  8:201412/1‑149 (MQ=255)
                                                                                                                                 tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                >  1:133380/1‑149 (MQ=255)
                                                                                                                                 tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                >  6:176424/1‑149 (MQ=255)
                                                                                                                                 tAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                >  7:92759/1‑149 (MQ=255)
                                                                                                                                  aaaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGACGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc               >  1:59476/1‑149 (MQ=255)
                                                                                                                                   aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGcc               <  5:292088/148‑1 (MQ=255)
                                                                                                                                   aaGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGc                <  5:176424/147‑1 (MQ=255)
                                                                                                                                    aGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCa             >  4:204997/1‑149 (MQ=255)
                                                                                                                                       aaagCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                        <  5:279199/135‑1 (MQ=255)
                                                                                                                                       aaagCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCaaa                        >  6:279199/1‑135 (MQ=255)
                                                                                                                                       aaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgct                                                           >  5:263520/1‑100 (MQ=255)
                                                                                                                                       aaaaCCCCGGCCCAGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCcttcgct                                                           <  6:263520/100‑1 (MQ=255)
                                                                                                                                         aaCCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACctct        >  2:132586/1‑149 (MQ=255)
                                                                                                                                           ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTAc                              >  2:96832/1‑125 (MQ=255)
                                                                                                                                           ccccGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTAc                              <  1:96832/125‑1 (MQ=255)
                                                                                                                                              cGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAg   <  6:207834/149‑1 (MQ=255)
                                                                                                                                               ggCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGATGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg  <  1:226192/149‑1 (MQ=255)
                                                                                                                                               ggCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAgg  <  7:179651/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATCAGGATCAGCAACAGAATGGCCGGCTTCAAGGACATGGCAGCCTCGCACTCAGGTGGTCTGAAGTCGAGTGTAGACGGGCATTGCCGCAGGTGTGCAACGGAAGGCGGGCTTTATGTCACGGGCCACTAAAGAAAAACCCCGGCCTGGGCCGGGGTTTTGTCATCACTCGTCGAGGAAGGAGCGCAGGTGCTCGCTGCGAGTCGGGTGGCGCAGCTTGCGCAACGCCTTCGCTTCGATCTGACGGATCCGCTCACGCGTTACGTCAAACTGCTTGCCCACCTCTTCGAGG  >  NC_002947/470919‑471210

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: ‑