Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,724,0990TG35.3% 10.3 / 17.3 17intergenic (+17/‑29)pssA/yfiMphosphatidylserine synthase; CDP‑diacylglycerol‑serine O‑phosphatidyltransferase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base G (3/3);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGC  >  NC_000913_3_pae_pgi/2723954‑2724239
                                                                                                                                                 |                                                                                                                                            
gAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACgg                                                                                                                                          <  2:154495/150‑1 (MQ=255)
    aTACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCACATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                      >  1:200717/1‑150 (MQ=255)
      accaccATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGgtttg                                                                                                                                    >  2:22391/1‑150 (MQ=255)
      accaccATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCACATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttgt                                                                                                                                    <  2:200717/150‑3 (MQ=255)
                                    cAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCtt                                                                                                      >  4:80085/1‑150 (MQ=255)
                                         tATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGt                                                                                                 <  3:80085/150‑1 (MQ=255)
                                                               gTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                          >  2:165177/1‑103 (MQ=255)
                                                               gTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                          <  1:165177/103‑1 (MQ=255)
                                                                     aaaCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCTTACAGGACGGGgtt                                                                                                                                      <  1:107168/85‑1 (MQ=255)
                                                                     aaaCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgtt                                                                                                                                      >  2:107168/1‑85 (MQ=255)
                                                                                            tATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTg                                              >  3:119465/1‑150 (MQ=255)
                                                                                                   aTCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCggg                                       <  4:119465/150‑1 (MQ=255)
                                                                                                       aCCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAgg                                   >  2:74822/1‑150 (MQ=255)
                                                                                                                                    aaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCa                                     <  3:22651/118‑1 (MQ=255)
                                                                                                                                    aaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCa                                     >  4:22651/2‑119 (MQ=255)
                                                                                                                                     aaCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGa     <  1:74822/150‑1 (MQ=255)
                                                                                                                                        cccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGc  >  1:15995/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GAGCATACCACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGATGC  >  NC_000913_3_pae_pgi/2723954‑2724239

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