Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,629,456 A→C 8.6% *322Y (TAA→TAC mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,629,4560AC 8.6% 202.2 / 8.0 67*322Y (TAA→TACmucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base C (4/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTCCA  >  NC_002947/1629311‑1629604
                                                                                                                                                 |                                                                                                                                                    
aGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATgg                                                                                                                                                   >  7:259519/1‑149 (MQ=255)
   tAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCATTACTGgcgc                                                                                                                                                <  6:195995/149‑1 (MQ=255)
       cccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCtt                                                                                                                                            <  8:212335/149‑1 (MQ=255)
        ccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCttt                                                                                                                                           >  6:2434/1‑149 (MQ=255)
           aCCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATTCCCAGGCGCGCCAGTAATGGCGCGCTTTTGc                                                                                                                                        <  4:225850/149‑1 (MQ=255)
           aCCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGAAGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGc                                                                                                                                        >  5:169680/1‑149 (MQ=255)
              tCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTg                                                                                                                                     >  1:13770/1‑149 (MQ=255)
              tCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACCCGACGCCCTGTCGATGCGCACCCCGGAGGCCCAGCCGCGCCAGTAATGGCGCGCTTTTGCTTg                                                                                                                                     >  2:209776/1‑149 (MQ=255)
                ggCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                                          >  6:3758/1‑138 (MQ=255)
                ggCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                                          <  5:3758/142‑5 (MQ=255)
                 gCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGcc                                                                                                                                   >  3:50130/1‑148 (MQ=255)
                   aGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAg                                                                                                                                <  4:159089/149‑1 (MQ=255)
                         gcgTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAg                                                                                                                           <  1:247389/148‑1 (MQ=255)
                               ccTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                                    <  7:29601/149‑1 (MQ=255)
                               ccTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                                    <  4:50130/149‑1 (MQ=255)
                                cTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTgg                                                                                                                   >  7:247240/1‑149 (MQ=255)
                                cTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTgg                                                                                                                   <  5:55960/149‑1 (MQ=255)
                                   aaTACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGc                                                                                                                 <  3:279891/148‑1 (MQ=255)
                                    aTACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCa                                                                                                                >  7:78991/1‑148 (MQ=255)
                                        acCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGc                                                                                                           >  1:271601/1‑149 (MQ=255)
                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                         >  8:21747/1‑149 (MQ=255)
                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                         >  3:13082/1‑149 (MQ=255)
                                           cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCt                                                                                                         >  2:300703/1‑148 (MQ=255)
                                               ggCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaat                                                                                                    >  2:297209/1‑149 (MQ=255)
                                                gCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTACTGGCGCGCTTTTGCGTGCCGAGCCAGGGCTGGTGGGCAGACGCCtgaaat                                                                                                    >  8:255341/1‑148 (MQ=255)
                                                gCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGc                                                                                                                 >  8:226448/1‑135 (MQ=255)
                                                      tGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTa                                                                                             <  7:64939/149‑1 (MQ=255)
                                                        aTGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTaaa                                                                                           <  8:247240/149‑1 (MQ=255)
                                                        aTGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTaaa                                                                                           <  8:25724/149‑1 (MQ=255)
                                                        aTGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTaa                                                                                            <  7:21747/148‑1 (MQ=255)
                                                         tGGGGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGTCGCGCTTTTGCGTGTCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAg                                                                                          <  1:209776/149‑1 (MQ=255)
                                                          ggTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAg                                                                                          >  5:305785/1‑148 (MQ=255)
                                                          ggTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGACCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAg                                                                                          >  5:286713/1‑148 (MQ=255)
                                                             gACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCAtt                                                                                      >  2:115944/1‑149 (MQ=255)
                                                               ccGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTg                                                                                     >  8:271597/1‑148 (MQ=255)
                                                                cGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGt                                                                                    >  3:303123/1‑148 (MQ=255)
                                                                    gTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAttt                                                                               <  5:137153/149‑1 (MQ=255)
                                                                        gCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCaa                                                                           <  2:288740/149‑1 (MQ=255)
                                                                        gCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCaa                                                                           <  5:2434/149‑1 (MQ=255)
                                                                        gCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCaa                                                                           <  4:303123/149‑1 (MQ=255)
                                                                          gAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATc                                                                         >  6:153304/1‑149 (MQ=255)
                                                                                 cGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTg                                                                              <  6:30333/137‑1 (MQ=255)
                                                                                 cGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTg                                                                              >  5:30333/1‑137 (MQ=255)
                                                                                        tACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAAtttt                                                           <  2:271601/149‑1 (MQ=255)
                                                                                        tACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATGGTCATTTGCAATCCATCGGCAAAtttt                                                           >  6:96543/1‑149 (MQ=255)
                                                                                         aCGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTc                                                          >  8:99070/1‑149 (MQ=255)
                                                                                          cGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTc                                                          <  6:305785/148‑1 (MQ=255)
                                                                                              aGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAg                                                     <  7:99070/149‑1 (MQ=255)
                                                                                                  cgcgTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCa                                                 <  5:153304/149‑1 (MQ=255)
                                                                                                         cccTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTg                                          <  8:78991/149‑1 (MQ=255)
                                                                                                          ccTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTg                                          >  8:172604/1‑148 (MQ=255)
                                                                                                           cTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGaa                                        <  4:13082/149‑1 (MQ=255)
                                                                                                              tCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAaga                                      <  6:169680/148‑1 (MQ=255)
                                                                                                              tCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAag                                       <  1:286437/147‑1 (MQ=255)
                                                                                                                  tGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTg                                 <  1:300703/149‑1 (MQ=255)
                                                                                                                        cccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCttt                           >  6:116926/1‑149 (MQ=255)
                                                                                                                        cccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCttt                           >  6:116924/1‑149 (MQ=255)
                                                                                                                                  ccAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgt                 >  4:87462/1‑149 (MQ=255)
                                                                                                                                   cAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTATCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg                >  8:76097/1‑149 (MQ=255)
                                                                                                                                    aagcgcgcCATTACTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg                >  7:154547/3‑148 (MQ=255)
                                                                                                                                    aagcgcgcCATTACTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtgtg                >  7:226008/3‑148 (MQ=255)
                                                                                                                                        gcgcCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAggg           <  1:115944/149‑1 (MQ=255)
                                                                                                                                        gcgcCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAggg           <  7:255341/149‑1 (MQ=255)
                                                                                                                                        gcgcCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAggg           <  6:286713/149‑1 (MQ=255)
                                                                                                                                            cAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGa                                         <  3:120181/115‑1 (MQ=255)
                                                                                                                                            cAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGa                                         >  4:120181/1‑115 (MQ=255)
                                                                                                                                             aGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCtt      >  6:13363/1‑149 (MQ=255)
                                                                                                                                                 aTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTcc   <  1:297209/148‑1 (MQ=255)
                                                                                                                                                 aTGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTCCa  <  7:271597/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTGTGCAGGGGTCTTTCCA  >  NC_002947/1629311‑1629604

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: