Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,299,069 C→T 6.4% intergenic (+15/+25) pdxH → / ← PP_1130 pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,299,0690CT 6.4% 831.1 / 31.1 263intergenic (+15/+25)pdxH/PP_1130pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (137/109);  new base T (9/8);  total (146/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 286 total aligned reads displayed.

CTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCTG  >  NC_002947/1298921‑1299218
                                                                                                                                                    |                                                                                                                                                     
cTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCa                                                                                                                                                      >  1:2198515/1‑150 (MQ=255)
 tGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCGGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCaa                                                                                                                                                     >  2:253006/1‑150 (MQ=255)
   aGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                   <  1:194236/150‑1 (MQ=255)
   aGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                   <  1:6668117/150‑1 (MQ=255)
      aTTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGc                                                                                                                                                <  1:4193125/150‑1 (MQ=255)
         gggggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGcccc                                                                                                                                             <  2:4555738/150‑1 (MQ=255)
          ggggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCa                                                                                                                                            <  1:2748417/150‑1 (MQ=255)
          ggggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCa                                                                                                                                            >  1:5040788/1‑150 (MQ=255)
            ggggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgt                                                                                                                                          <  2:3603950/150‑2 (MQ=255)
               gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAt                                                                                                                                       <  2:7094243/150‑1 (MQ=255)
               gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAt                                                                                                                                       >  2:6084085/1‑150 (MQ=255)
               gATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAt                                                                                                                                       <  1:6939420/150‑1 (MQ=255)
                  aCCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGa                                                                                                                                    >  1:3280445/1‑150 (MQ=255)
                      tttGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                                <  1:6333329/150‑1 (MQ=255)
                      tttGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                                <  2:5159216/150‑1 (MQ=255)
                              ccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcccc                                                                                                                        >  2:1354446/1‑150 (MQ=255)
                              ccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcccc                                                                                                                        >  2:5495302/1‑150 (MQ=255)
                                 aGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGTTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGa                                                                                                                     >  1:5946871/1‑150 (MQ=255)
                                        cGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTg                                                                                                              <  1:5797239/150‑1 (MQ=255)
                                        cGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTg                                                                                                              >  1:4476029/1‑150 (MQ=255)
                                         gAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                             >  1:860087/1‑150 (MQ=255)
                                         gAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgc                                                                                                             <  2:3470285/150‑1 (MQ=255)
                                           gTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTCCCCCGGATAACCGgccg                                                                                                           >  1:3339388/1‑150 (MQ=255)
                                               tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTc                                                                                                       <  1:6176748/150‑1 (MQ=255)
                                                 gCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTcgc                                                                                                     <  2:4811990/150‑1 (MQ=255)
                                                  cAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATCACCTGCCGCCTcgcg                                                                                                    >  2:5934695/1‑150 (MQ=255)
                                                     gggcgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTc                                                                                                 <  2:1246028/150‑1 (MQ=255)
                                                        cgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCaa                                                                                              <  1:1485924/150‑1 (MQ=255)
                                                        cgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCaa                                                                                              <  1:6438530/150‑1 (MQ=255)
                                                            cAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                          <  1:2045217/150‑1 (MQ=255)
                                                             aaGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccacc                                                                                         <  2:4071499/150‑1 (MQ=255)
                                                              aGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAccc                                                                                        >  1:1691467/1‑150 (MQ=255)
                                                                  gtcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGcc                                                                                    >  2:6417942/1‑150 (MQ=255)
                                                                   tcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCa                                                                                   <  2:6597953/150‑1 (MQ=255)
                                                                   tcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGAGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCa                                                                                   >  2:5885038/1‑150 (MQ=255)
                                                                     tGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGt                                                                                 >  1:3108134/1‑150 (MQ=255)
                                                                      gCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                                <  1:3848383/150‑1 (MQ=255)
                                                                       cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTg                                                                               >  2:1920987/1‑150 (MQ=255)
                                                                       cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTg                                                                               >  1:977473/1‑150 (MQ=255)
                                                                          gATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCtt                                                                            >  2:4229209/1‑150 (MQ=255)
                                                                            tCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTAc                                                                          >  2:6215606/1‑150 (MQ=255)
                                                                             cGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACg                                                                         <  1:3595134/150‑1 (MQ=255)
                                                                              gTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGc                                                                        <  2:1534806/150‑1 (MQ=255)
                                                                                ttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGccc                                                                      >  1:5373083/1‑150 (MQ=255)
                                                                                  gAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCa                                                                    >  2:7307406/1‑150 (MQ=255)
                                                                                   aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCt                                                                                                                             >  1:2897205/1‑92 (MQ=255)
                                                                                    aCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAcc                                                                  >  2:1158787/1‑150 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCtt                                                                >  2:6935336/1‑150 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCtt                                                                >  2:2688192/1‑150 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCtt                                                                >  1:5568547/1‑150 (MQ=255)
                                                                                      tATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCtg                                                                >  2:6528041/1‑149 (MQ=255)
                                                                                        tCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtct                                                              <  2:6356923/150‑1 (MQ=255)
                                                                                         cGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtctc                                                             <  2:1629760/150‑1 (MQ=255)
                                                                                             tggtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGc                                                         >  2:5683492/1‑150 (MQ=255)
                                                                                               gtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAg                                                                                                                                           <  2:2039694/66‑1 (MQ=38)
                                                                                               gtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCcg                                                       <  1:688053/150‑1 (MQ=255)
                                                                                                tggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCcgc                                                      <  1:6749892/150‑1 (MQ=255)
                                                                                                  gACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGc                                                                              <  1:3405014/124‑1 (MQ=255)
                                                                                                  gACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGc                                                                              >  2:3405014/1‑124 (MQ=255)
                                                                                                     gggCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtct                                                              <  1:1216859/137‑1 (MQ=255)
                                                                                                     ccccAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTg                                                 >  1:1667650/4‑150 (MQ=255)
                                                                                                       gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGcccc                                                                                                                                             <  2:318173/56‑1 (MQ=37)
                                                                                                       gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgcg                                               >  1:469031/1‑150 (MQ=255)
                                                                                                       gCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgcg                                               >  1:1382572/1‑150 (MQ=255)
                                                                                                           tGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCaa                                           >  1:4902238/1‑150 (MQ=255)
                                                                                                            ggCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAc                                          <  1:4463155/150‑1 (MQ=255)
                                                                                                             gCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAcc                                         >  2:4702737/1‑150 (MQ=255)
                                                                                                             gCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAcc                                         >  1:483309/1‑150 (MQ=255)
                                                                                                               tGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCtt                                       >  1:3059233/1‑150 (MQ=255)
                                                                                                                gCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCttt                                      >  2:2990746/1‑150 (MQ=255)
                                                                                                                gCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCttt                                      <  2:2349659/150‑1 (MQ=255)
                                                                                                                  gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTcc                                    >  1:7093855/1‑150 (MQ=255)
                                                                                                                  gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTcc                                    >  2:6773988/1‑150 (MQ=255)
                                                                                                                       cGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCAt                               <  1:6546702/150‑1 (MQ=255)
                                                                                                                       cGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCAt                               <  2:3954631/150‑1 (MQ=255)
                                                                                                                           tGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                           >  2:4274850/1‑150 (MQ=255)
                                                                                                                                cccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgcg                      <  1:1566902/150‑1 (MQ=255)
                                                                                                                                  cTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                           <  1:2295231/143‑1 (MQ=255)
                                                                                                                                   tcaatggcTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGttt                   >  2:1734659/8‑150 (MQ=255)
                                                                                                                                   tGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGttt                   >  1:3800459/1‑150 (MQ=255)
                                                                                                                                     aGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGttttt                 <  1:6599754/150‑1 (MQ=255)
                                                                                                                                      gCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTc                <  2:927591/150‑1 (MQ=255)
                                                                                                                                      gCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTc                <  1:7359596/150‑1 (MQ=255)
                                                                                                                                       cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa               <  2:2175292/150‑1 (MQ=255)
                                                                                                                                           tGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGAt           >  1:752911/1‑150 (MQ=255)
                                                                                                                                           tGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGAt           >  2:3118166/1‑150 (MQ=255)
                                                                                                                                           tGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGAt           >  1:5487181/1‑150 (MQ=255)
                                                                                                                                            ggggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATc          <  2:53216/150‑1 (MQ=255)
                                                                                                                                             gggCCGCTTGGCGGCCCCAGTCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCg         >  1:4098798/1‑150 (MQ=255)
                                                                                                                                             gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCg         >  2:748314/1‑150 (MQ=255)
                                                                                                                                               gccgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGc       >  2:5996914/1‑150 (MQ=255)
                                                                                                                                               gccgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGc       >  1:5329564/1‑150 (MQ=255)
                                                                                                                                                  gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                                >  1:4003144/1‑74 (MQ=255)
                                                                                                                                                  gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGcccc    >  1:1139004/1‑150 (MQ=255)
                                                                                                                                                  gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTGTCAGGATCGGcccc    >  1:4375727/1‑150 (MQ=255)
                                                                                                                                                   ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   <  2:4447163/150‑1 (MQ=255)
                                                                                                                                                   ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   >  1:3689058/1‑150 (MQ=255)
                                                                                                                                                   ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   >  1:1089678/1‑150 (MQ=255)
                                                                                                                                                   ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt   >  2:2855748/1‑150 (MQ=255)
                                                                                                                                                   ccAAGCGGCCCCAGCCATTGATCAAGCAACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGGCAACCATACCAGCAGGAACGCGCGTATTTCAGGATCGGCCCCt   >  1:606657/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCTG  >  NC_002947/1298921‑1299218

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: