Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,777,5630AC40.4% 63.4 / 22.5 52intergenic (+255/+37)PP_5476/PP_1585hypothetical protein/antidote protein
Reads supporting (aligned to +/- strand):  ref base A (5/26);  new base C (21/0);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 6.85e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGA  >  NC_002947/1777434‑1777694
                                                                                                                                 |                                                                                                                                   
cTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCg                                                                                                                               <  2:849256/136‑1 (MQ=255)
     cGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgcg                                                                                                                          <  1:524304/136‑1 (MQ=255)
        cccAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCCTgggggcc                                                                                                                       >  1:1283861/1‑131 (MQ=255)
            ggCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGt                                                                                                                   <  1:821815/136‑1 (MQ=255)
                 tGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAtt                                                                                                              <  1:712214/136‑1 (MQ=255)
                 tGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAtt                                                                                                              <  1:689161/136‑1 (MQ=255)
                     gCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATc                                                                                                          <  1:733077/136‑1 (MQ=255)
                       cGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                        <  2:86144/136‑1 (MQ=255)
                          agagCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGCGCCCCGGGTTATTGTTAAggtt                                                                                                     >  2:528837/1‑134 (MQ=255)
                            agCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCCTGGCGCCCCTGTTTTTTGATCAGGAtgt                                                                                                   >  2:736124/1‑136 (MQ=255)
                                   aGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGAGGCGCCCCGGGTTATTTATAAAGTTTTGTTCaaa                                                                                            >  2:1160076/1‑135 (MQ=255)
                                   aGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGGGGCGCCCCAGTTTATTGATAAGGATGTGTTCaac                                                                                            >  2:331056/1‑135 (MQ=255)
                                   aGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCCCCACCCCCGAGGCGCCCCGGGTTTTTTATAAGGAAGTGTTAAAg                                                                                            >  1:431722/1‑136 (MQ=255)
                                     gTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCg                                                                                          <  2:927568/136‑1 (MQ=255)
                                      tctttGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGa                                                                                         <  1:939371/133‑1 (MQ=255)
                                           gtgtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGa                                                                                    <  1:772929/136‑1 (MQ=255)
                                             gtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCGGGGCGCCCCTGTTTATTGGAGAGGATGTTGTAAAGCGAccccacg                                                                                  >  1:1214742/1‑131 (MQ=255)
                                             gtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgcg                                                                                                                          >  2:922520/1‑96 (MQ=255)
                                             gtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgcg                                                                                                                          <  1:922520/96‑1 (MQ=255)
                                                 aTTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGGGCCCCGGGTTTTTTTTGAGGTTGTTTTGAATGGACCAGCCggggg                                                                              >  1:260656/1‑132 (MQ=255)
                                                 aTTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGCGCCCCAGGTTATTGATAAGGTTGTGGTGAATCGCCCAGACGCTgg                                                                              >  2:647930/1‑136 (MQ=255)
                                                   tGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGGGGGGCCCCCGTTTTTTTATCAGGTTGTGTTGAAGCGACCAGACGCggggc                                                                            >  1:729244/1‑134 (MQ=255)
                                                    ggggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGGGCCCCCCCCCGGGGGGCCCGGGTTTTTGATGATGTTGTGTTAAATCGGCCCGCCGCGGggcc                                                                           >  1:132043/1‑136 (MQ=255)
                                                    ggggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCCTGGGGCCCCGGTTTTTTGATCAGGATGTGTTAAAGCGACCCGGAGCGGgacc                                                                           >  2:1373541/1‑136 (MQ=255)
                                                    ggggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCCCCCCGCCCCGGGGGCCCCGGGTTTTTTTTGAGGATGTTTTGAATCGCCCAGCCGCGGggcc                                                                           >  1:43183/1‑136 (MQ=255)
                                                      ggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGCGCCCCGGGTTATTGATGAGGATGTGTTCAAGCGCCCAGACGCTGgaccgc                                                                         >  2:1046226/1‑135 (MQ=255)
                                                       gTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCGGGGGGCCCCGGGGTTTTGTTGAGGATGTGTTTTAGCGACCACCCGCTGggccggc                                                                        >  1:308432/1‑134 (MQ=255)
                                                       gTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgac                                                                        <  1:992062/136‑1 (MQ=255)
                                                         tGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGAccc                                                                      <  1:1107518/136‑1 (MQ=255)
                                                         tGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCCGGGCGCCCCTGTTTATTGATCAGGCTGTGTTGAAGCGGCCAGGCGCTGGGCCGAccc                                                                      >  2:857497/1‑136 (MQ=255)
                                                           cAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCGGGGGGCCCCGGGTTTTTGTTGAGGATGTGTTGAAGCGACCCGACGCGGGGCCGCCCCgc                                                                    >  1:22738/1‑135 (MQ=255)
                                                                  ccccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCCCCCCGCCCGGGGGGCCCCCGGTTATTGATCAAGATGTGTTAAAACGACCAGACGCGGGACCCCCCCcatgatata                                                             >  2:17712/1‑133 (MQ=255)
                                                                  ccccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCCCCCCCCCCCGGGCGCCCCGGGTTTTTGAACAGGATGTGTTAAAGCGACCCCACGCAGGGCCGCCCCCATGCtctc                                                             >  2:646515/1‑136 (MQ=255)
                                                                   cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCCGGGGGCCCCAGTTTTTTTATAAAGATGTGTTTAATCGACCAGACCCAGGGCCGACCCcatgataaca                                                            >  1:730985/1‑132 (MQ=255)
                                                                     cAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCCTGGCGCCCCGGGTTTTTTATAAGGATGTGTTGAAGCGACCAGACGCGGGACCGACCCCATGATCTAAta                                                          >  2:1192011/1‑135 (MQ=255)
                                                                         ctctcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCg                                                      <  2:94604/136‑1 (MQ=255)
                                                                         ctctcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCg                                                      <  1:1192011/136‑1 (MQ=255)
                                                                         ctctcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCg                                                      <  1:451128/136‑1 (MQ=255)
                                                                              cGTTTTTCGATGTCAGTCGCGTGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTg                                                 <  2:275684/136‑1 (MQ=255)
                                                                               gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCGTGGCGCCCCGGTTTTTTTATAAGGATGTGTTAAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGtttg                                                >  1:410542/1‑134 (MQ=255)
                                                                               gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCCCCCGGGGCGCCCCGGGTTTTTTATCAGGATGTGTTAAAGCGCCCAGCCGCTGGGCCGACCCCATGATCTGATCGCCGGGTTgg                                                >  2:240969/1‑135 (MQ=255)
                                                                                    tCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTc                                           <  2:470496/136‑1 (MQ=255)
                                                                                     cGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCGCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCg                                          <  1:1064134/136‑1 (MQ=255)
                                                                                                      ttATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                         <  2:1219916/136‑1 (MQ=255)
                                                                                                      ttATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                         <  2:458404/136‑1 (MQ=255)
                                                                                                                 cGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCa              <  1:1293419/136‑1 (MQ=255)
                                                                                                                      cTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTg         <  2:1341264/136‑1 (MQ=255)
                                                                                                                       tCGCGGCGCCCCGCCCGGGGCGCCCCGGTTTATTGTTCAGGGTGTGTTGAAGCGGCCAGGCGCTGGACCGCCCCGATGATCTCATCGCCGGGGTGCTGTTCGGCGGTGGGCAGATGGAAGGCCGGCTGCAGGTTGa        >  2:1293419/1‑136 (MQ=255)
                                                                                                                        cgcgGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAg       <  2:498639/136‑1 (MQ=255)
                                                                                                                        cgcgGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAg       <  1:838835/136‑1 (MQ=255)
                                                                                                                             cgcCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGCTCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGa  >  1:1262509/1‑136 (MQ=255)
                                                                                                                             cgcCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGa  <  2:813777/136‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
CTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGA  >  NC_002947/1777434‑1777694

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