Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,140,3520GC21.9% 71.7 / 15.2 32R150G (CGC→GGC) hybEhydrogenase 2‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (10/15);  new base C (4/3);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTG  >  NC_000913/3140218‑3140490
                                                                                                                                      |                                                                                                                                          
ccAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGtgtg                                                                                                                              <  4:150238/149‑1 (MQ=255)
              cacGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGct                                                                                                                                  >  1:82314/1‑130 (MQ=255)
              cacGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGTCGCGCCCTGct                                                                                                                                  <  2:82314/131‑2 (MQ=255)
                cggcggTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGt                                                                                                              >  7:25447/1‑149 (MQ=255)
                        aCGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCa                                                                                                      <  4:32105/149‑1 (MQ=255)
                                   aTCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGtgtg                                                                                                                              <  3:10948/114‑1 (MQ=255)
                                   aTCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGtgtg                                                                                                                              >  4:10948/1‑114 (MQ=255)
                                      gTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAt                                                                                         <  8:52163/148‑1 (MQ=255)
                                         ctgctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGt                                                                                     <  8:25447/149‑1 (MQ=255)
                                            ctCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTc                                                                                                                <  8:64683/119‑1 (MQ=255)
                                               cGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                             <  8:165606/103‑1 (MQ=255)
                                               cGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                             >  7:165606/1‑103 (MQ=255)
                                                    gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt                                                                                                                               <  5:180218/96‑1 (MQ=255)
                                                    gCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCgtgt                                                                                                                               >  6:180218/1‑96 (MQ=255)
                                                                    aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTg                                                                                    >  7:135968/1‑123 (MQ=255)
                                                                    aaCGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTg                                                                                    <  8:135968/123‑1 (MQ=255)
                                                                                     gagaCAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCt                                         >  1:204265/1‑149 (MQ=255)
                                                                                                   aTTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                             <  7:146661/51‑1 (MQ=255)
                                                                                                   aTTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTgg                                                                                                                             >  8:146661/1‑51 (MQ=255)
                                                                                                          ccACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGTCGTTgg                                                             <  8:1757/108‑1 (MQ=255)
                                                                                                          ccACTCCTGCGACCAAACAGCAGGGCGCCACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCAATCGTGCGCAGTCATCCGTCAGGCGTTgg                                                             >  7:1757/1‑108 (MQ=255)
                                                                                                              tCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATc                                                 >  3:140263/1‑116 (MQ=255)
                                                                                                              tCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATc                                                 <  4:140263/116‑1 (MQ=255)
                                                                                                                        cacacgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggc                            >  8:1572/6‑127 (MQ=255)
                                                                                                                        cacacgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAgcggc                            <  7:1572/122‑1 (MQ=255)
                                                                                                                        cacacgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga           >  4:63021/6‑144 (MQ=255)
                                                                                                                        cacacgcagGGCGTCGCGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGga           <  3:63021/139‑1 (MQ=255)
                                                                                                                        aaaCAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAg             >  6:118346/1‑142 (MQ=255)
                                                                                                                        aaaCAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAg             <  5:118346/142‑1 (MQ=255)
                                                                                                                            agcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCctct                                                       <  3:58493/96‑1 (MQ=255)
                                                                                                                            agcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCctct                                                       >  4:58493/1‑96 (MQ=255)
                                                                                                                            agcagGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTg  <  2:204265/149‑1 (MQ=255)
                                                                                                                                     cgACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCTGTCAGGCGTTGGCCCTCTTCAATCGACATGCtgtg                                      <  2:41493/104‑1 (MQ=255)
                                                                                                                                     cgACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCTGTCAGGCGTTGGCCCTCTTCAATCGACATGCtgtg                                      >  1:41493/1‑104 (MQ=255)
                                                                                                                                      |                                                                                                                                          
CCAATTTCCAGCCACACGGCGGTGACGCGCTTAACATCGTGCTGCTCCGCCTGCCGTTGGATAATTTCAACGGCGCTCTGGCAAAGAGACAACTCATGCATTTTCGCCACTCCTGCGACCAAACAGCAGGGCGCGACGCCCTGCGTGTGGTACATCCGGATTCGTGACTGGCAGCGAAAGGATCATTCGTGCGCAGTCATCCGTCAGGCGTTGGCCCTCTTCAATCGACATGCTGTGCGAAAGCGGCGACATCAGCGAACAGGAGAGATATTG  >  NC_000913/3140218‑3140490

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