Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,490,3400CT64.7% 7.8 / 10.0 17G164S (GGT→AGT) pabAaminodeoxychorismate synthase subunit 2
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (5/6);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACA  >  NC_000913/3490134‑3490551
                                                                                                                                                                                                              |                                                                                                                                                                                                                   
aCTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCggg                                                                                                                                                                         <  1:267870/251‑1 (MQ=255)
 cTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATAAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAAt                                                                                                                                                                                          >  1:171591/1‑233 (MQ=255)
                                          gCCATCTCATGATCACCCTTTTACGCATAAACAAATGTGAAATTATAACCACAAAAAATGCATAAAAAATCACTAAAGGGCAATCAGAAGGCAGCGATGCAGGAAATTACCCAGCAGTGGAGTTCCTTGTTCGCTAAGAATCCTTTCTGGATGCAACGGCACACTGTCCAGCTCCCACTGGCGAGGGCGAATCCCCATAATCTCCCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAg                                                                                                                                <  1:186082/250‑1 (MQ=255)
                                                 cATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAACAAATCACTAAATGGCAATCAGAAATAAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGATAAGAATAATTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAACGCACGCTGGGAATGAGTCAGGTTc                                                                                                                        >  1:234508/1‑251 (MQ=255)
                                                                  gCATACACACATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGt                                                                                                       >  1:269787/1‑251 (MQ=255)
                                                                                     tataACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTcc                                                                                                                                                                                                                 <  1:74968/126‑1 (MQ=255)
                                                                                     tataACCACAAAATATGCATAAAAAANCACTAANNGGNAATCAGAAATCAGCGATGCNGGNNNTNNNCNCGCAGNTGNNGTCCTTGTTNGCTAAGAATACTTTCTGGATGGAACTGCACACCTTcc                                                                                                                                                                                                                 >  2:74968/1‑126 (MQ=255)
                                                                                                                              cAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATTATCTCGCGGGTGTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGGAAGGAGTCAGTTTCCACCACCAGCGAATGGTAGCTTGTCACGGTAAGTTGAGTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTa                                            <  1:287425/250‑1 (MQ=255)
                                                                                                                                         gATGCAGGAAATTAGCCAGCAGTGGAGGGCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACtggcgatggcga                                                                                                                                                                                            <  1:49360/95‑1 (MQ=255)
                                                                                                                                             cAGGAAATTAGCCAGCAGTTGATGTACTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAAtc                                                                                                                                                                                >  1:250547/1‑103 (MQ=255)
                                                                                                                                                   aTTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATTCTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGTTAGCGTGTCACGGTAAGTGGATTTTCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAAGCGTCGATGTTTTGCCGTGcn                      <  1:96232/251‑2 (MQ=255)
                                                                                                                                                    ttAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATAGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATAc                                                               >  1:132534/1‑209 (MQ=255)
                                                                                                                                                              nttGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGc                                      >  1:193121/2‑224 (MQ=255)
                                                                                                                                                              gTTGATGTCCTTGTTCGCTAAGAATACTTTATGGATGGAACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCTGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTccacca                                                                                                                   >  1:272191/1‑147 (MQ=255)
                                                                                                                                                                                          ttCTGGATGGAACTGCACACTTTCCAGATCCCACtggcgatggcg                                                                                                                                                                                             <  1:128367/45‑1 (MQ=255)
                                                                                                                                                                                              ggatggaACTGCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGGCAAAGCACGCTTGTAATGAGTCAGGTTCCACCCCCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGc                                                                          <  1:161537/156‑1 (MQ=255)
                                                                                                                                                                                                        gCACACTTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGCATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGcaca  >  1:187828/1‑218 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                   
ACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGATCACCCTGTTACGCATAAACAAATGTGAAATTATAACCACAAAATATGCATAAAAAATCACTAAATGGCAATCAGAAATCAGCGATGCAGGAAATTAGCCAGCAGTTGATGTCCTTGTTCGCTAAGAATACTTTCTGGATGGAACTGCACACCTTCCAGATCCCACTGGCGATGGCGAATCCCCATAATCTCTCGGGTTTCGCTCCAGGCCGTCACGTCAAAGCACGCTGGTAATGAGTCAGGTTCCACCACCAGCGAATGGTAGCGTGTCACGGTAAGTGGATTTGCCAGCCCCCGAAATACGCCCTCACCGTTATGTGTAATCGGCGAGGTTTTGCCGTGCATGACCTTTGCGGCGCGCACA  >  NC_000913/3490134‑3490551

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