Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,002,963 G→T 5.2% H138N (CAT→AAT)  rfbD ← dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,002,9630GT 5.2% 847.5 / 24.7 248H138N (CAT→AAT) rfbDdTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (116/119);  new base T (3/10);  total (119/129)
Fisher's exact test for biased strand distribution p-value = 8.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GGGTGCGCCAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTGCTAATA  >  NC_002947/2002814‑2003112
                                                                                                                                                     |                                                                                                                                                     
gggTGCGCCAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATg                                                                                                                                                       <  2:844601/150‑1 (MQ=255)
 ggTGCGCCAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGc                                                                                                                                                      >  1:2362941/1‑150 (MQ=255)
 ggTGCGCCAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGc                                                                                                                                                      >  2:2254357/1‑150 (MQ=255)
       ccAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCt                                                                                                                                                >  2:169366/1‑150 (MQ=255)
        cAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGCATGCCATTCTCGccag                                                                                                                                               <  1:615164/150‑3 (MQ=255)
        cAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGCATGCCATTCTCGccag                                                                                                                                               <  1:4398498/150‑3 (MQ=255)
        cAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCtt                                                                                                                                               >  2:3059255/1‑150 (MQ=255)
        cAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCAGGCGAGCATGCCATTCTCGccag                                                                                                                                               <  1:2760700/150‑3 (MQ=255)
         aaaTTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTc                                                                                                                                              >  2:452343/1‑150 (MQ=255)
         aaaTTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCAGGCGAGCATGCCATTCTCGccagc                                                                                                                                              <  1:5156223/150‑4 (MQ=255)
          aaTTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTcc                                                                                                                                             <  1:5043960/150‑1 (MQ=255)
             tGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGc                                                                                                                                          >  1:2256958/1‑150 (MQ=255)
               gTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCtt                                                                                                                                        <  2:3002310/150‑1 (MQ=255)
                tCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTg                                                                                                                                       <  2:4955190/150‑1 (MQ=255)
                tCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTg                                                                                                                                       <  1:2547000/150‑1 (MQ=255)
                tCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTg                                                                                                                                       <  2:2573520/150‑1 (MQ=255)
                 cTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTgg                                                                                                                                      >  1:2869964/1‑150 (MQ=255)
                  tGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGt                                                                                                                                     >  1:2107273/1‑150 (MQ=255)
                   gCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGAGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTc                                                                                                                                    <  2:2845225/150‑1 (MQ=255)
                     tACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCg                                                                                                                                  >  2:323223/1‑150 (MQ=255)
                      aCCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCgg                                                                                                                                 >  2:2020504/1‑150 (MQ=255)
                         aCGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCg                                                                                                                              <  2:551542/150‑1 (MQ=255)
                              cAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGAc                                                                                                                         <  2:3901937/150‑1 (MQ=255)
                              cAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGAc                                                                                                                         <  2:2481702/150‑1 (MQ=255)
                               aaCGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGAcc                                                                                                                        >  1:1601375/1‑150 (MQ=255)
                               aaCGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGAcc                                                                                                                        >  2:1556889/1‑150 (MQ=255)
                               aaCGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGAcc                                                                                                                        >  1:1926164/1‑150 (MQ=255)
                                    ctcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAg                                                                                                                   >  1:3465495/1‑150 (MQ=255)
                                    ctcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGATCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAg                                                                                                                   >  1:997963/1‑150 (MQ=255)
                                     tcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgg                                                                                                                  >  2:3831917/1‑150 (MQ=255)
                                     tcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgg                                                                                                                  <  2:3982819/150‑1 (MQ=255)
                                     tcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgg                                                                                                                  >  2:4322590/1‑150 (MQ=255)
                                     tcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgg                                                                                                                  >  2:2216540/1‑150 (MQ=255)
                                     tcACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgg                                                                                                                  >  2:1333917/1‑150 (MQ=255)
                                      cACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAggg                                                                                                                 <  2:3650079/150‑1 (MQ=255)
                                      cACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAggg                                                                                                                 <  1:4503065/150‑1 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTg                                                                                                                                       <  1:4250030/127‑1 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgggg                                                                                                                >  2:314644/1‑150 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgggg                                                                                                                >  1:4082641/1‑150 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgggg                                                                                                                <  1:4800900/150‑1 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAgggg                                                                                                                >  1:5264253/1‑150 (MQ=255)
                                       aCGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCGGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGACTGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGACCgggg                                                                                                                >  2:1128833/1‑150 (MQ=255)
                                           tCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCg                                                                                                            >  2:518461/1‑150 (MQ=255)
                                             cgccgcCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGt                                                                                                          <  1:3216012/150‑1 (MQ=255)
                                               ccgccAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAg                                                                                                        >  1:4590820/1‑150 (MQ=255)
                                                    aGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCa                                                                                                   >  2:2069455/1‑150 (MQ=255)
                                                    aGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGGGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCa                                                                                                   <  1:1128833/150‑1 (MQ=255)
                                                        gcagcaTGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGttt                                                                                               >  1:2033499/1‑150 (MQ=255)
                                                           gcaTGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTtcct                                                                                            <  2:2344703/150‑1 (MQ=255)
                                                            caTGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTtcctc                                                                                           >  1:5454116/1‑150 (MQ=255)
                                                            caTGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTtcctc                                                                                           <  2:3228617/150‑1 (MQ=255)
                                                            caTGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTtcctc                                                                                           >  2:331614/1‑150 (MQ=255)
                                                                   ttGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAgcgc                                                                                    >  1:5607127/1‑150 (MQ=255)
                                                                            ttGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGcc                                                                                                                               <  2:3453363/98‑1 (MQ=255)
                                                                            ttGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGActgccg                                                                           >  2:4626376/1‑150 (MQ=255)
                                                                                 cccGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGActgccgctgcc                                                                      <  2:710567/150‑1 (MQ=255)
                                                                                  ccGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGActgccgctgcac                                                                     >  1:1777262/1‑148 (MQ=255)
                                                                                   cgcgcgCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGActgccgctgccgg                                                                    >  2:504437/1‑149 (MQ=255)
                                                                                   cgcgcgCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGt                                                                    >  1:414539/1‑150 (MQ=255)
                                                                                       cgcAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCaaa                                                                <  1:3154145/150‑1 (MQ=255)
                                                                                       cgcAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCaaa                                                                <  2:1740873/150‑1 (MQ=255)
                                                                                       cgcAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACGGCAGCGGCAGTCagc                                                                <  2:4776033/150‑3 (MQ=255)
                                                                                         cagcaTACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATa                                                              <  1:907801/150‑1 (MQ=255)
                                                                                           gcaTACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACa                                                            <  1:2927597/150‑1 (MQ=255)
                                                                                                 acCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCg                                                      >  2:1057566/1‑150 (MQ=255)
                                                                                                 acCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCg                                                      >  1:2754568/1‑150 (MQ=255)
                                                                                                 acCAAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCg                                                      >  2:1666740/1‑150 (MQ=255)
                                                                                                   ccAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTtcct                                                                                            >  2:1170404/1‑110 (MQ=255)
                                                                                                   ccAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGt                                                    <  1:718115/150‑1 (MQ=255)
                                                                                                    cAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTg                                                   >  1:3032930/1‑150 (MQ=255)
                                                                                                    cAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTg                                                   <  1:3969443/150‑1 (MQ=255)
                                                                                                     aGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTgg                                                  >  2:4426864/1‑150 (MQ=255)
                                                                                                     aGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTgg                                                  <  2:3545985/150‑1 (MQ=255)
                                                                                                          gTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAg                                             <  1:1081569/150‑1 (MQ=255)
                                                                                                          gTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAg                                             <  2:5030410/150‑1 (MQ=255)
                                                                                                           tGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGt                                            <  2:1438680/150‑1 (MQ=255)
                                                                                                           tGCGCAGGACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGt                                            <  2:2478131/150‑1 (MQ=255)
                                                                                                             cgcAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGa                                          <  1:2383356/150‑1 (MQ=255)
                                                                                                                   aCCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCTTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAAc                                    <  2:418397/150‑1 (MQ=255)
                                                                                                                    ccACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAAcc                                   >  1:2987562/1‑150 (MQ=255)
                                                                                                                     cACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCa                                  <  1:169366/150‑1 (MQ=255)
                                                                                                                     cACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCa                                  <  1:1398216/150‑1 (MQ=255)
                                                                                                                     cACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCa                                  <  1:1192043/150‑1 (MQ=255)
                                                                                                                          ccTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGgcgc                             <  1:1316336/150‑1 (MQ=255)
                                                                                                                            ttGGCGCCGCTGGCGAGAATGGCGTGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGccc                           >  2:929004/1‑150 (MQ=255)
                                                                                                                            ttGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGccc                           >  1:2231127/1‑150 (MQ=255)
                                                                                                                               gcgcCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCaaa                        <  1:377395/150‑1 (MQ=255)
                                                                                                                                   cgcTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAgctg                    >  2:5163753/1‑150 (MQ=255)
                                                                                                                                   cgcTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAgctg                    >  1:1361208/1‑150 (MQ=255)
                                                                                                                                    gcTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAgctgc                   <  1:3955166/150‑1 (MQ=255)
                                                                                                                                      aggCGAGCATGCCATTCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGctgctg                 >  2:4821950/2‑150 (MQ=255)
                                                                                                                                        gCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGtt               >  1:3134837/1‑150 (MQ=255)
                                                                                                                                            gaATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCAc           <  1:3167782/150‑1 (MQ=255)
                                                                                                                                            gaATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCAc           <  1:3580683/150‑1 (MQ=255)
                                                                                                                                              aTGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGt         >  1:5215914/1‑150 (MQ=255)
                                                                                                                                              aTGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGt         <  2:2373185/150‑1 (MQ=255)
                                                                                                                                               tgccATTCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTg        >  1:4722800/4‑150 (MQ=255)
                                                                                                                                               tGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTg        <  1:4322590/150‑1 (MQ=255)
                                                                                                                                                   aTGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTGCTaa    >  2:550093/1‑150 (MQ=255)
                                                                                                                                                     gCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTGCTAATa  <  1:3385786/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTGCGCCAAATTGGTCTGCTACCACGCTCAACGTCTCACGCTCCGCCGCCAGGCGCAGCATGGTCTTGGCAAAATTGTGCCCGCGCGCAGCATACACCCAGCTGGTGCGCAGTACCACGGCCTTGGCGCCGCTGGCGAGAATGGCATGCTCGCCTTCCAGCTTGGTCCGGCCGTAGACCGAAAGGGGGCCGGTAGGCGCAGTTTCCTCCCAGCGCTGACTGCCGCTGCCGTCAAATACATAATCGGTGGAATAGTGAATCAACCAGGCGCCCAAAGCTGCTGTTTCACGTGCTAATA  >  NC_002947/2002814‑2003112

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