Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,421,2030TC53.5% 23.3 / 65.3 43H8H (CAT→CAC) yhdYamino‑acid transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (14/6);  new base C (0/23);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 4.95e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-08
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGC  >  NC_007779/3421103‑3421302
                                                                                                    |                                                                                                   
ttATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCAt                                                                                                     >  3:1396582/1‑101 (MQ=255)
 tATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATc                                                                                                    <  3:1451053/101‑1 (MQ=255)
 tATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATc                                                                                                    <  3:1533721/101‑1 (MQ=255)
   tCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATccc                                                                                                  >  1:687426/1‑101 (MQ=255)
        cTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCcgcg                                                                                             >  1:381365/1‑101 (MQ=255)
           aCTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGccc                                                                                          <  2:1439/101‑1 (MQ=255)
             tATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCgg                                                                                        <  2:112983/101‑1 (MQ=255)
                   gctgctGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGcc                                                                                  >  3:15318/1‑101 (MQ=255)
                            gaATATCTATACCCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCaa                                                                         <  3:702449/101‑1 (MQ=255)
                                  cTATACCCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCcgcg                                                                   <  4:157661/101‑1 (MQ=255)
                                   tataCCCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCAGAACTCAAGCcgcgc                                                                  <  4:118965/101‑1 (MQ=255)
                                      acccgccgCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGAt                                                               <  3:564451/99‑1 (MQ=255)
                                         cgccgcATGGCGATCGTTGAACGCTAGGGAGTCATAATGACAAAAGTTTTGCTGTCTCCCCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGgt                                                            <  3:559502/101‑1 (MQ=255)
                                         cgccgcATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGgt                                                            <  1:1372634/101‑1 (MQ=255)
                                          gccgcATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGgtg                                                           >  1:1278911/1‑101 (MQ=255)
                                          gccgcATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGgtg                                                           <  4:123997/101‑1 (MQ=255)
                                           ccgcATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGgtgt                                                          >  3:333680/1‑101 (MQ=255)
                                              cATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTggg                                                       >  3:186296/1‑101 (MQ=255)
                                               accgcgATCGTTGAACGCTAAGGAGTCATGATGAAAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGt                                                      <  3:880556/99‑1 (MQ=255)
                                               accgcgATCGTTAAACGCTAAGGGGTATGGAAGACAAAAGTTTTGCTGTCCCTCCCCCCGCGCCGGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGt                                                      <  1:1433252/99‑1 (MQ=255)
                                               aTCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCTCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGt                                                      >  3:1273842/1‑101 (MQ=255)
                                               aTCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGt                                                      >  1:1457679/1‑101 (MQ=255)
                                                 cgcgATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGc                                                    <  3:1474110/101‑1 (MQ=255)
                                                 cgcgATCGTTGAACGCGAAGGAGTCATGATGACAAAAGTTTTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGc                                                    <  3:694291/101‑1 (MQ=255)
                                                    gATCGTTAAACGGTAAGGGGTCAGGATGACAAAAGTTTGGCTGCCCCATCCCCCGCGCCCGGGGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGaa                                                 <  3:933944/101‑1 (MQ=255)
                                                     aTCGTGAAACGCTAAGAAGTGATGATGAAAAAAGTATTGCTGCCTCCCCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGaaa                                                <  3:1147423/101‑1 (MQ=255)
                                                      tCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCCCCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGCGGGTGCGaaaa                                               <  1:1382683/101‑1 (MQ=255)
                                                      tCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGaaaa                                               >  3:911267/1‑101 (MQ=255)
                                                      tCGTTGAACGCTAAGGAGTCATGATGAAAAAAGTTTTGCTGTCTCACCCCCCGCGCCNGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGaaan                                               <  3:1367679/101‑2 (MQ=255)
                                                          taaaCGCTAAGGAGTCATGATGAAAAAAGTTTTGCTGCCTCACCCCCCGCGCCGGGGGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATc                                           <  3:1466659/99‑1 (MQ=255)
                                                          tGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATc                                           <  1:448100/101‑1 (MQ=255)
                                                                       gTCATGATGACAAAAGTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTg                              <  3:681525/101‑1 (MQ=255)
                                                                        tCATGATGCCAAACGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                             >  4:66608/1‑101 (MQ=255)
                                                                        tCATGATAACAAAATTATTGCTGTCCCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                             <  3:1472557/101‑1 (MQ=255)
                                                                        tCAGGATGAGAAATGTTTTCTTGTCTCGCCCCCCGGCCCCGGCGAGCCATATCTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                             <  3:384561/101‑1 (MQ=255)
                                                                        tCAGGATGACAAAGGTATGGCTGTCTCACCCCCCGCCCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTgg                             <  1:1042617/101‑1 (MQ=255)
                                                                                      gttttGCTGCCTCCCCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgct               <  4:65299/98‑1 (MQ=255)
                                                                                      gTATTGCTGTCTCACCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCATCGGTATCAATAGCctgct               <  4:25000/101‑1 (MQ=255)
                                                                                          gtcTGTCTCACCCCCCGCCCCCGGGGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACt           <  1:925730/99‑1 (MQ=255)
                                                                                           gCTGCCGCACCCCCCGCGCCGGGGGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTa          <  3:1096195/101‑1 (MQ=255)
                                                                                                ctcATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGc     >  1:1214251/1‑101 (MQ=255)
                                                                                                 tcATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCt    >  3:338434/1‑101 (MQ=255)
                                                                                                   aTCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGc  >  1:683954/1‑101 (MQ=255)
                                                                                                   aTCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGGGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGc  >  1:957976/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
TTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGC  >  NC_007779/3421103‑3421302

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