Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,527,6020GC28.6% 88.4 / 20.4 42intergenic (+450/‑300)yncH/rhsEIPR020099 family protein/pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base G (24/6);  new base C (0/12);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT  >  NC_000913/1527537‑1527669
                                                                 |                                                                   
ccACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCg                                                                   >  2:506172/1‑68 (MQ=255)
  aCTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGt                                                                 <  2:491627/68‑1 (MQ=255)
      ttAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg                                                             <  1:868916/68‑1 (MQ=255)
      taaaGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCg                                                             <  1:642946/66‑1 (MQ=255)
       tAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGt                                                            >  2:206595/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:266048/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  2:481641/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:44814/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:470024/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  2:717741/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:548186/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:574563/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  2:175027/1‑68 (MQ=255)
          gCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGtata                                                         >  1:724184/1‑68 (MQ=255)
            tAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATAtt                                                       >  1:206097/1‑68 (MQ=255)
             aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTc                                                      <  2:44814/68‑1 (MQ=255)
             aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTc                                                      <  1:742292/68‑1 (MQ=255)
             aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc                                                      <  2:590659/68‑1 (MQ=255)
             aGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGCACCGTATATTc                                                      <  2:649403/68‑1 (MQ=255)
             aGTGTTTAACGACGTTAAGATGGCAATGTGCAGACCATGTTATCAGGGTGGGCCGGTAGCGTATATTc                                                      <  2:903420/68‑1 (MQ=255)
             aGTGGTTAACGACGTTAAGATGGCAACGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTc                                                      <  2:526848/68‑1 (MQ=255)
               tgtTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCt                                                    >  2:269749/1‑68 (MQ=255)
                 tttAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGAAGCGTATATTCCTta                                                  <  1:426248/68‑1 (MQ=255)
                  ttAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTtat                                                 <  2:724182/68‑1 (MQ=255)
                         cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAAcc                                          <  1:206595/68‑1 (MQ=255)
                         cgTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGCCGGTAGCGTACATTCCTTATACAAAcc                                          <  2:548771/68‑1 (MQ=255)
                                         tgCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATg                          >  1:95958/1‑68 (MQ=255)
                                          gCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACACACCATCATACAATTCAATGa                         >  2:456564/1‑68 (MQ=255)
                                                    ttATCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa               <  2:574563/68‑1 (MQ=255)
                                                      aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa               >  2:298920/1‑66 (MQ=255)
                                                      aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa               <  1:298920/66‑1 (MQ=255)
                                                      aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAat             >  2:114898/1‑68 (MQ=255)
                                                      aTCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGATATTGAat             >  2:865513/1‑68 (MQ=255)
                                                       tCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTTTTGAatt            >  2:903536/1‑68 (MQ=255)
                                                       tCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGATATTGAatt            >  1:592001/1‑68 (MQ=255)
                                                       tCAGGGTGGGCCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAatt            <  1:480152/68‑1 (MQ=255)
                                                         aGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattat          <  1:376112/68‑1 (MQ=255)
                                                          gggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAattatt         >  1:887747/1‑68 (MQ=255)
                                                           ggTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAAttattt        >  1:294209/1‑68 (MQ=255)
                                                            gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa               >  1:929056/1‑60 (MQ=255)
                                                            gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGa               <  2:929053/60‑1 (MQ=255)
                                                            gTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAAttattta       >  1:690592/1‑68 (MQ=255)
                                                                 gCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTAttgtt  >  2:272398/1‑68 (MQ=255)
                                                                 |                                                                   
CCACTATTAAGCTAGTGTTTAACGACGTTAAGATGGCAATGTGCAGATCATGTTATCAGGGTGGGGCGGTAGCGTATATTCCTTATACAAACCATCATTCATTTCAATGATTGTTATTGAATTATTTATTGTT  >  NC_000913/1527537‑1527669

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