Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 593,437 A→C 6.4% E391A (GAG→GCG)  PP_0503 → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947593,4370AC 6.4% 1226.4 / 37.5 397E391A (GAG→GCG) PP_0503MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (173/194);  new base C (11/19);  total (184/213)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 399 total aligned reads displayed.

TTGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCAG  >  NC_002947/593288‑593585
                                                                                                                                                     |                                                                                                                                                    
 tGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAg                                                                                                                                                     >  1:4063751/1‑150 (MQ=255)
   cTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgc                                                                                                                                                   >  2:6624273/1‑150 (MQ=255)
    tATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgcc                                                                                                                                                  >  2:4742096/1‑150 (MQ=255)
       tCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGgctgcc                                                                                                                                               <  2:2939501/150‑1 (MQ=255)
         gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCt                                                                                                                                             >  2:7442908/1‑150 (MQ=255)
         gCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCt                                                                                                                                             >  1:5369114/1‑150 (MQ=255)
          cTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCTGCCGCCCtc                                                                                                                                            <  2:7305726/150‑2 (MQ=255)
          cTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTg                                                                                                                                            <  1:1644033/150‑1 (MQ=255)
           tGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCTCGGCCGCCCtcg                                                                                                                                           <  2:121572/150‑3 (MQ=255)
               cAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATg                                                                                                                                       >  2:390107/1‑150 (MQ=255)
               cAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATg                                                                                                                                       >  1:6554354/1‑150 (MQ=255)
                aGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGt                                                                                                                                      >  1:1381043/1‑150 (MQ=255)
                  cGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTc                                                                                                                                    <  2:4586761/150‑1 (MQ=255)
                    gTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAg                                                                                                                                  >  1:4285035/1‑150 (MQ=255)
                     tttGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAgg                                                                                                                                 >  2:1457595/1‑150 (MQ=255)
                              ggTGTTCGCGAGATGGGTGATTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCACTGTTCAGGTCAGAACtt                                                                                                                        >  2:4670146/1‑150 (MQ=255)
                               gtgtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTg                                                                                                                       <  1:2315062/150‑1 (MQ=255)
                                 gtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                                     <  1:7344422/150‑1 (MQ=255)
                                 gtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                                     <  2:6439245/150‑1 (MQ=255)
                                 gtTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTa                                                                                                                     <  1:5699870/150‑1 (MQ=255)
                                  ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                                    <  1:390107/150‑1 (MQ=255)
                                  ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                                    <  1:7237131/150‑1 (MQ=255)
                                  ttGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                                    >  2:6334159/1‑150 (MQ=255)
                                  gtGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTTTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGTCCGCCCTGCAATGTTCAGGTCAGAACTTGTAg                                                                                                                    <  1:2866377/149‑1 (MQ=255)
                                    ggCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGcc                                                                                                                  >  1:7111368/1‑150 (MQ=255)
                                      cGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                                <  2:275583/150‑1 (MQ=255)
                                      cGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                                <  1:4596044/150‑1 (MQ=255)
                                      cGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                                <  2:2608631/150‑1 (MQ=255)
                                      cGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCaa                                                                                                                <  1:2751757/150‑1 (MQ=255)
                                                gtTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCAt                                                                                                      >  2:5659061/1‑150 (MQ=255)
                                                 tttGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATg                                                                                                     <  2:3996042/150‑1 (MQ=255)
                                                 ttggCACACTGGGGGTGTGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGTGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCTGATCTTGTAGCCAAGACCGACCATg                                                                                                     <  1:6214029/147‑1 (MQ=255)
                                                    gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtat                                                                                                  <  2:2212285/150‑1 (MQ=255)
                                                    gCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGtat                                                                                                  >  2:5441812/1‑150 (MQ=255)
                                                       cacTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATAcc                                                                                               <  2:4158945/150‑1 (MQ=255)
                                                          tGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCAc                                                                                            <  2:712729/150‑1 (MQ=255)
                                                          tGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCAc                                                                                            <  2:3538726/150‑1 (MQ=255)
                                                          tGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCAc                                                                                            <  2:901623/150‑1 (MQ=255)
                                                            gggCTGGGTTGGGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACgg                                                                                          <  1:1236553/150‑1 (MQ=255)
                                                             ggCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACggg                                                                                         >  1:5821385/1‑150 (MQ=255)
                                                              gCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgt                                                                                        >  1:1875517/1‑150 (MQ=255)
                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAgaccgac                                                                                                         <  1:4922801/132‑1 (MQ=255)
                                                               cTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtc                                                                                       <  2:3344380/150‑1 (MQ=255)
                                                                  ggTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                                    >  2:1041776/1‑150 (MQ=255)
                                                                     tGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtc                                                                                 <  2:3401355/150‑1 (MQ=255)
                                                                        ggTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                              <  2:7066864/150‑1 (MQ=255)
                                                                        ggTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                              <  2:3164845/150‑1 (MQ=255)
                                                                        agtgtgtGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgac                                                                              <  1:3095140/149‑1 (MQ=255)
                                                                           gtgtgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtc                                                                           <  2:1305282/150‑1 (MQ=255)
                                                                              tgggctgggcATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgacgtcgacgtcgac                                                                        <  2:247443/150‑1 (MQ=255)
                                                                              tgGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgtcgtcgacgtcgac                                                                        >  1:3541649/1‑150 (MQ=255)
                                                                               gggCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAcc                                                                       >  2:2809685/1‑150 (MQ=255)
                                                                                 gCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCtt                                                                     >  2:2449056/1‑150 (MQ=255)
                                                                                  ctgggCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGCGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTg                                                                    <  1:2665022/147‑1 (MQ=255)
                                                                                  cTTGGCCTGGGCATCAGTGCCGCCGCACTGTGCTACTGTCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCTACGTCGACGTCGACCTTg                                                                    <  1:5431957/150‑1 (MQ=255)
                                                                                   ttGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCACCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTgg                                                                   <  2:435250/150‑1 (MQ=255)
                                                                                    tGGCATGGGCATCAGTGCCGCCGCACTGTGCTTCTGGCTTTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGt                                                                  <  1:1013477/150‑1 (MQ=255)
                                                                                    tGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTTTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGt                                                                  <  2:5365738/150‑1 (MQ=255)
                                                                                      gCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCt                                                                >  2:6131004/1‑150 (MQ=255)
                                                                                        aTGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTg                                                              >  1:6418254/1‑150 (MQ=255)
                                                                                        aTGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTg                                                              >  1:7379838/1‑150 (MQ=255)
                                                                                           tagcacAGTGCGGCGGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGCCGTCGACGGCGACCCGGGTCATGCTg                                                           >  1:4705380/6‑150 (MQ=255)
                                                                                           ggCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTg                                                           >  2:4845040/1‑150 (MQ=255)
                                                                                             cATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGCCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTa                                                         <  1:3382774/150‑1 (MQ=255)
                                                                                              aTCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAg                                                        >  2:6066860/1‑150 (MQ=255)
                                                                                               tCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGc                                                       >  1:6094045/1‑150 (MQ=255)
                                                                                                  gggccgccGCACTGTGCTACGGGCTCTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa                                                    <  1:796541/148‑1 (MQ=255)
                                                                                                  gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa                                                    >  1:5023738/1‑150 (MQ=255)
                                                                                                  gTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCaa                                                    <  2:1665961/150‑1 (MQ=255)
                                                                                                    gccgccGCACTGTTCTCCGGGCTATGGCTGCACCGCGCGCGGGCGGCCGAGCCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAgg                                                  <  1:3136009/150‑1 (MQ=255)
                                                                                                    gccgccGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAgg                                                  <  2:5011809/150‑1 (MQ=255)
                                                                                                    gccgccGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAgg                                                  <  1:282854/150‑1 (MQ=255)
                                                                                                       gccgcACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCGTGGTCTTGCTGTAGCCAAGGGCg                                               >  1:1812379/1‑150 (MQ=255)
                                                                                                         cgcCCTGTGCTACTGGCTATGGCTGCACCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGt                                             <  1:5550391/150‑1 (MQ=255)
                                                                                                              tgtgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGc                                        <  2:5916071/150‑1 (MQ=255)
                                                                                                              tgtgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGc                                        >  2:2323923/1‑150 (MQ=255)
                                                                                                              tgtgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGc                                        >  2:5538558/1‑150 (MQ=255)
                                                                                                                tgCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTTGCCAAGGGCGGGCGGGCCg                                      >  1:435036/1‑150 (MQ=255)
                                                                                                                  cTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCATGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCggt                                    >  1:454523/1‑148 (MQ=255)
                                                                                                                   tACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCGTGCTGTAGACAAGGGCGGTCGGGCCGGTg                                   >  1:3525273/1‑150 (MQ=255)
                                                                                                                     cTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGAGCGTGGTGATGGTGTAGGCAAGGGCGGTCGGGCCgagatc                                 >  1:1516189/1‑145 (MQ=255)
                                                                                                                           aTGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACgg                                                                                          >  2:6203560/1‑87 (MQ=255)
                                                                                                                             ggCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCt                         >  2:6854290/1‑150 (MQ=255)
                                                                                                                               cTGCAACGCGCGCGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                       <  2:3586775/150‑1 (MQ=255)
                                                                                                                               cTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTg                       <  2:1971143/150‑1 (MQ=255)
                                                                                                                                   aaCGCGCGTGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgtgt                   <  2:7217091/150‑1 (MQ=255)
                                                                                                                                   aaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGgtgt                   <  2:1571417/150‑1 (MQ=255)
                                                                                                                                   aaCGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGGCGTCGTCGACGGCGTCCGTGGTCTTGGTGTATCCAAGGGGGGGCGGGCCGGTGCCGCTGGCCTTGgtgt                   >  1:3295120/1‑150 (MQ=255)
                                                                                                                                      gcgcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcga                <  2:524328/150‑1 (MQ=255)
                                                                                                                                        gcgAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatg              <  2:6960825/150‑1 (MQ=255)
                                                                                                                                           aGGGCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCGTGGtgtcgatgtcg           >  2:1756526/1‑150 (MQ=255)
                                                                                                                                            gggCGGCCTCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCATGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGGTGACGCTGGCCTTGGtgtcgatgtaga          >  1:6589088/1‑150 (MQ=255)
                                                                                                                                            gggCGGCCGCGCCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCACGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGtgtcgatgtcga          <  1:2464460/150‑1 (MQ=255)
                                                                                                                                                 gcCGCGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATCCCCCCGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc     <  1:2036445/150‑1 (MQ=255)
                                                                                                                                                 gcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATCCCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCACGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc     <  2:4481938/150‑1 (MQ=255)
                                                                                                                                                 gcCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAc     <  1:3613744/150‑1 (MQ=255)
                                                                                                                                                  ccGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGGCGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCATGGtgtcgatgtcgtgccact    >  1:5779783/1‑143 (MQ=255)
                                                                                                                                                  ccGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTGGACGTCGACCTTGGTGTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCGTTGGTGTCGATGTCGTCGTTcc    >  2:1870171/1‑150 (MQ=255)
                                                                                                                                                  ccGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTAcc    >  2:4918947/1‑150 (MQ=255)
                                                                                                                                                    gAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGgtcgac                                                                                    >  1:5321221/1‑68 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTGCTATTCGCTGGGCAGCGGTTTGGGGGCGGTGTTGGCGAGCTGGGTGTTTGCACACTGGGGCTGGGTTGCGGTGTGTGGGCTTGGCATGGGCATCAGTGCCGCCGCACTGTGCTACTGGCTATGGCTGCAACGCGCGAGGGCGGCCGAGGCCGCCCTGCAATGTTCAGGTCAGAACTTGTAGCCAAGACCGACCATGTATACCCACGGGTCGACGTCGACGTCGACCTTGGTCTTGCTGTAGCCAAGGGCGGTCGGGCCGTTGACGCTGGCCTTGGTGTCGATGTCGACGTACCAG  >  NC_002947/593288‑593585

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