Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1500TG23.9% 94.0 / 11.4 46intergenic (+133/‑65)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/11);  new base G (0/11);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 5.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAC  >  NC_000913/3781079‑3781230
                                                                        |                                                                                 
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)
                                                                cggggggaTGCGGGGAAACGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa          <  2:821351/77‑1 (MQ=255)
                                                                  cGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt       <  2:821348/81‑1 (MQ=255)
                                                                    ggggggAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc      <  1:821349/81‑1 (MQ=255)
                                                                       gTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa   >  1:821408/1‑81 (MQ=255)
                                                                        tGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc  >  2:821418/1‑81 (MQ=255)
                                                                        |                                                                                 
TGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAC  >  NC_000913/3781079‑3781230

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