Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1450CG22.2% 119.1 / 14.3 45intergenic (+128/‑70)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (24/11);  new base G (0/10);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAT  >  NC_000913/3781065‑3781225
                                                                                |                                                                                  
cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                                    <  2:821325/81‑1 (MQ=255)
 gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                  >  2:821390/1‑81 (MQ=255)
  ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                 <  2:821326/81‑1 (MQ=255)
              tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                    <  1:821328/81‑1 (MQ=255)
              tgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                    <  2:821327/81‑1 (MQ=255)
                taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                   <  1:821329/80‑1 (MQ=255)
                       tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                           >  1:821409/1‑81 (MQ=255)
                         aCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTg                                                         >  2:821388/1‑81 (MQ=255)
                         aCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTg                                                         >  1:821397/1‑81 (MQ=255)
                            ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                       >  1:821405/1‑80 (MQ=255)
                             gaaaaCAACCCCCCCCTTTTTTTTCCCCCCCCCCCCCTCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTc                                                      <  1:821338/78‑1 (MQ=255)
                              gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                    >  2:821407/1‑81 (MQ=255)
                               aaaCGACACCCCACTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTc                                                      <  2:821330/78‑1 (MQ=255)
                                  cGACACCCACCTTTTTTTTCCCCCCCCCCCCCTCGGGGGAGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAAt                                                <  2:821337/81‑1 (MQ=255)
                                    acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGTGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                              <  2:821331/81‑1 (MQ=255)
                                    acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                              >  1:821391/1‑81 (MQ=255)
                                      acCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                            >  2:821392/1‑81 (MQ=255)
                                           aCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGCCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                       <  2:821332/81‑1 (MQ=255)
                                           aCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                       >  2:821393/1‑81 (MQ=255)
                                             tctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                     >  1:821410/1‑81 (MQ=255)
                                                  ttttCTCCCTCGCCCCCCGGGGGGGGGGCCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                 <  2:821335/80‑1 (MQ=255)
                                                  ttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                  <  1:821334/80‑1 (MQ=255)
                                                     tccccccCGCCCCGCGGGGGGGGGGGAGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                              <  1:821355/77‑1 (MQ=255)
                                                           tCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                        >  2:821395/1‑80 (MQ=255)
                                                            cGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                        >  2:821401/1‑79 (MQ=255)
                                                            cGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                        >  2:821412/1‑79 (MQ=255)
                                                             gCCCCTCGGGGGAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                      <  2:821336/80‑1 (MQ=255)
                                                             gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                     >  2:821399/1‑81 (MQ=255)
                                                             gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                     >  2:821413/1‑81 (MQ=255)
                                                             gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                     >  1:821402/1‑81 (MQ=255)
                                                             cccccccCGGGGGGGGGGCCGGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                      <  2:821340/80‑1 (MQ=255)
                                                              ccccccGGGGGAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                    <  1:821343/81‑1 (MQ=255)
                                                                cccccGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                  <  2:821342/78‑1 (MQ=255)
                                                                ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                  >  2:821403/1‑81 (MQ=255)
                                                                 cccggggaAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                 <  2:821341/77‑1 (MQ=255)
                                                                 cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                 >  2:821411/1‑81 (MQ=255)
                                                                  tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                 >  1:821396/1‑80 (MQ=255)
                                                                   ccGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                 >  1:821400/1‑79 (MQ=255)
                                                                    cGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg              >  1:821406/1‑81 (MQ=255)
                                                                      gggAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa            >  2:821421/1‑81 (MQ=255)
                                                                        gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg          >  2:821394/1‑81 (MQ=255)
                                                                          gggggCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt        <  2:821352/79‑1 (MQ=255)
                                                                          gagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc         >  2:821398/1‑80 (MQ=255)
                                                                           ggggCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa       <  2:821346/80‑1 (MQ=255)
                                                                           agGCCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa       <  2:821345/81‑1 (MQ=255)
                                                                            gggCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa      <  2:821347/81‑1 (MQ=255)
                                                                                cGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt  <  2:821348/81‑1 (MQ=255)
                                                                                |                                                                                  
CGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAT  >  NC_000913/3781065‑3781225

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: