Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,812,7920CA46.7% 56.3 / 17.8 46T21N (ACC→AAC) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/19);  new base A (21/0);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 3.06e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACC  >  NC_000913/3812642‑3812929
                                                                                                                                                      |                                                                                                                                          
cGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAc                                                                                                                                            <  2:879439/151‑1 (MQ=255)
   aaaGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAcccc                                                                                                                                         >  1:156537/1‑151 (MQ=255)
   aaaGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAaccc                                                                                                                                         >  2:878573/1‑151 (MQ=255)
     aGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTg                                                                                                                                       <  2:117017/151‑1 (MQ=255)
        gAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTTCCTTTAAAAACCCTGaa                                                                                                                                     >  1:197929/1‑150 (MQ=255)
        gAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGaa                                                                                                                                     <  1:984954/150‑1 (MQ=255)
               aGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTATAAAAACCCTGAACTGGTGcg                                                                                                                             >  1:893365/1‑151 (MQ=255)
               aGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGGACTAGGGcg                                                                                                                             >  1:289993/1‑151 (MQ=255)
                 gTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGcgtc                                                                                                                           <  2:840882/151‑1 (MQ=255)
                 gTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGGGcg                                                                                                                             >  1:933977/1‑149 (MQ=255)
                   gAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGcgtcg                                                                                                                          <  1:261530/150‑1 (MQ=255)
                   gAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAAATGGTGcg                                                                                                                             >  1:241709/1‑147 (MQ=255)
                               acaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTATAAAAACCCCGAACTGGGGCGGCGTTTgcgc                                                                                                                   >  2:963020/1‑145 (MQ=255)
                                caaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATcgcg                                                                                                            <  1:325701/151‑1 (MQ=255)
                                      tATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAAATGGGGCGTCGTTTGCGCGATCGCgggg                                                                                                         >  1:890124/1‑148 (MQ=255)
                                      tATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGGACTGGTGCGGCGTTTTCGCGGTCGCgggg                                                                                                         >  2:486853/1‑148 (MQ=255)
                                       aTTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGcc                                                                                                       >  1:823036/1‑149 (MQ=255)
                                         tgtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGGCCATAAAAGCCCTCAAATGGGGCGTCGTTTTGGCGATCGCgggg                                                                                                         >  1:800687/1‑145 (MQ=255)
                                          gtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACg                                                                                                   <  1:941619/150‑1 (MQ=255)
                                           tgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTAATGCTGCCGACAAAAACCCCGAACTGGGGCGGCGTTTGCGCGATGGCCGGGCCCACgg                                                                                                  >  2:540635/1‑149 (MQ=255)
                                             aTAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCCATAAAAACCCCGAAACGGTGCGGCGTTTTCGCGATCGCGGGGCCCACGTcc                                                                                                >  1:66372/1‑150 (MQ=255)
                                                       cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAGCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCAt                                                                                      <  1:540307/150‑1 (MQ=255)
                                                       cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGcc                                                                                        >  1:618918/1‑148 (MQ=255)
                                                         aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGGTTCCCCTAAAAACCCTGAACTGGTGCGGCGGTTGGGCGATCCCGGGGCAGAGGAACGCGGAGccgag                                                                                     >  2:465824/1‑146 (MQ=255)
                                                            tCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCCTTAAAAAACCTGCACTGGTGCGTCCTTTGCGCGATAGCGGGGGCGGCGGCCGCGGAAGCATGACCgc                                                                                >  2:813979/1‑150 (MQ=255)
                                                            tCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCg                                                                                 <  2:471375/150‑1 (MQ=255)
                                                             cTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTTAACTGGGGCGGCGGTTGCGGGGTCGCGGGGGCGACGTCCGCGGAGCCATGAcc                                                                                  >  2:834493/1‑148 (MQ=255)
                                                               ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGaa                                                                               <  1:818033/149‑1 (MQ=255)
                                                                 tgtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGGGCGTCGTTTGCGCGATCGCGGGGCCGACGGCCGCGTAGCCATGACCGAAGCgg                                                                           >  1:492858/1‑151 (MQ=255)
                                                                  gtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGCATTGTCGCCTCTAAAAACCCCGAAATTGTGCGTCGTTTGCGCGGACGCGGGGCCGGCGTCCGCGGCGCCACGACCGAAACgg                                                                           >  2:921751/1‑150 (MQ=255)
                                                                    cACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCa                                                                         <  1:937767/150‑1 (MQ=255)
                                                                      cccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGCTGCCTATAAAAACCCCGAACTGGGGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGCGGCCATGACCGAAGCgg                                                                           >  2:664987/1‑146 (MQ=255)
                                                                        cTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGGCCGCGTAGGCATGACCGAAGCGGCaaaa                                                                      >  1:635361/1‑149 (MQ=255)
                                                                         tACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGGGCGGCGGTTGCGCGGTCGCGGGGGCGACGTCCGCGTAGCCCTGACCGAAGCGGCAAAAGcc                                                                   >  2:700952/1‑151 (MQ=255)
                                                                                         atGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTa                                                   <  2:771824/151‑1 (MQ=255)
                                                                                          tGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGCTGCCCATAAAAACCCCGAACTGGGGCGGCGGTTTCGCGATCGCGGGGCCGACGGCCCCGGAGCCATGACCGAAGCGGCAAAAGCCCTTTTCAcccccccttg                                                  >  2:33036/1‑145 (MQ=255)
                                                                                             gCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCtt                                               <  2:343782/151‑1 (MQ=255)
                                                                                                tGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTTAACTGGTGCGTCGTTTGGGCGGTCGCGGGGGCGGCGTCCGCGGAGCCATGACCGGAGGGGCAAAAGCCTTTATTACCCCCCCTTGCtt                                               >  2:195879/1‑148 (MQ=255)
                                                                                                           aaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGATCTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTg                                <  1:842461/151‑1 (MQ=255)
                                                                                                           aaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTg                                <  2:502036/151‑1 (MQ=255)
                                                                                                              aTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGGTTTTTGGt                              >  1:97598/1‑151 (MQ=255)
                                                                                                              aTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGGTTTTgggt                             >  1:197680/1‑149 (MQ=255)
                                                                                                                      cGGCGTTAGCGGCGGTATTGCTGCCTATAAACCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCCCTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGt                      <  1:10328/150‑1 (MQ=255)
                                                                                                                      cGGCGTTAGCGGCGGGATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGGCGACGTCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCCCCTAGCTTGCAGGC‑ggggtcttggtattcc                        >  2:353986/1‑134 (MQ=255)
                                                                                                                        gCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTc                   <  2:776207/151‑1 (MQ=255)
                                                                                                                          gTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTcc                  <  2:66369/150‑1 (MQ=255)
                                                                                                                           ttAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTAt                           <  2:962520/140‑1 (MQ=255)
                                                                                                                                 gcggTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTct          <  2:654852/151‑1 (MQ=255)
                                                                                                                                         gctgcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAcc  >  1:755776/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACC  >  NC_000913/3812642‑3812929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑