Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,7480TG29.2% 31.7 / 14.6 24intergenic (‑172/+20)htpX/prcputative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (8/9);  new base G (3/4);  total (11/13)
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.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATC  >  NC_000913/1912601‑1912890
                                                                                                                                                   |                                                                                                                                              
cAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAtt                                                                                                                                               >  3:117492/1‑149 (MQ=255)
      gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTtgt                                                                                                                                         <  4:117492/149‑1 (MQ=255)
          ttttcGATACCCACATCGTATGGGTTACGCGTCTATTTTCAAGTCTGGATAGTGCGTAAGTGACCAGATAGACACAACTTTACATTTTGTAGCTTCTGATTTAGGGCATCTTGTCGCTGTTATATAATCAGGCTCAAGAAATTGTGCCt                                                                                                                                     <  6:128295/145‑1 (MQ=255)
          tatGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCt                                                                                                                                     >  1:161530/1‑149 (MQ=255)
            tGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                   <  2:161530/149‑1 (MQ=255)
             gaatctcccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtgat                                                                                                                                  <  7:17752/143‑1 (MQ=255)
                    ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtgat                                                                                                                                  >  8:17752/1‑142 (MQ=255)
                    ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                   >  2:152212/1‑141 (MQ=255)
                    ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAAGAAATTGTGCCtga                                                                                                                                   <  1:152212/141‑1 (MQ=255)
                                        gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgt                                                                                                          <  3:129871/146‑1 (MQ=255)
                                        gACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgt                                                                                                          >  4:129871/1‑146 (MQ=255)
                                                                         ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                    <  8:93847/103‑1 (MQ=255)
                                                                         ccAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                    >  7:93847/1‑103 (MQ=255)
                                                                          cAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                                                          >  3:60170/1‑144 (MQ=255)
                                                                          cAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGAt                                                                     <  4:60170/149‑1 (MQ=255)
                                                                                                                aCGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                          <  7:51742/138‑1 (MQ=255)
                                                                                                                aCGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                          >  8:51742/1‑138 (MQ=255)
                                                                                                                  ggCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCgg                             >  8:4074/1‑149 (MQ=255)
                                                                                                                    cATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGa                                                                      <  5:116477/106‑1 (MQ=255)
                                                                                                                    cATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGa                                                                      >  6:116477/1‑106 (MQ=255)
                                                                                                                              cTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc                                            <  2:31429/122‑1 (MQ=255)
                                                                                                                                   aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATcc                               >  5:22087/1‑130 (MQ=255)
                                                                                                                                   aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATcc                               <  6:22087/130‑1 (MQ=255)
                                                                                                                                             gCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcatc  <  7:4074/149‑1 (MQ=255)
                                                                                                                                             gCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATACACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatcatc  >  4:7868/1‑149 (MQ=255)
                                                                                                                                                   tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAg        >  4:10597/1‑137 (MQ=255)
                                                                                                                                                   tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAg        <  3:10597/137‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATC  >  NC_000913/1912601‑1912890

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