Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,544,0211.A73.9% 56.6 / 17.4 23coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base A (8/9);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAT  >  NC_000913_3_pae_tpiA/1543879‑1544149
                                                                                                                                               |                                                                                                                                 
ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGCCTGATATCCACAGATTTCTGCGAGCAATATGTTTTCCTTTATCTTCCGAGAAGGCCGGACTTCTTGGGTCCCATTCGCG‑A‑AAAAg                                                                                                                            >  2:499070/1‑149 (MQ=255)
                cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTAttttt                                                                                                             <  1:62684/149‑1 (MQ=255)
                 tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtc                                                                                                            >  2:166722/1‑149 (MQ=255)
                       aaaTAGCATCCAGACACAGAAGGCATGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTGTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAtgg                                                                                                      <  1:558505/149‑3 (MQ=255)
                               tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATACCGACTATTTTTATCATTTTGCactac                                                                                              >  1:514992/1‑145 (MQ=255)
                                 cAGACACAGAAGGCCAGAAATAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                            >  2:384544/1‑149 (MQ=255)
                                  aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcg                                                                                                                                    <  1:346220/109‑1 (MQ=255)
                                  aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcg                                                                                                                                    >  2:346220/1‑109 (MQ=255)
                                                  aagTAAACAACTGACTGATATCCAGATATTTCTTCGACCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                           <  2:358624/149‑1 (MQ=255)
                                                    gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                         >  2:389225/1‑149 (MQ=255)
                                                          aaCTGACTGATATCCAGATATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                   <  2:282463/149‑1 (MQ=255)
                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                   <  2:566276/149‑1 (MQ=255)
                                                               actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                              >  1:295110/1‑149 (MQ=255)
                                                                       ccAGAGATTTGTTCGAGAAATGTGTTTTCCTTTATTTTGCCAGAAGGCCGGAGTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTGTTTTTCTCATTTTGCAGTGCCGTATTGTTCCTCACATGTACACATTGGTAATGaaaaaaa                                                      <  2:547969/149‑1 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                  <  2:599745/135‑1 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                  >  1:599745/1‑135 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                    >  2:196664/1‑149 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                    <  1:564702/149‑1 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                           >  2:504779/1‑149 (MQ=255)
                                                                                                   cgtttATTTTCCCAGAAGGTTGGATTTTCTGGTTACCAGTCGCGAA‑AAAAGATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCCCACATTCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                          <  2:110474/147‑1 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATACCGACTAGTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                        >  1:584958/1‑149 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                        >  1:596021/1‑149 (MQ=255)
                                                                                                        aTTTTCCCAGAAGGCCGGATTTTCTGGTTTGTACTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACAAGAGGTAAGGCGCAATAGCCAGtt                     >  2:445936/1‑149 (MQ=255)
                                                                                                                          aTTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGaa   <  1:166722/149‑1 (MQ=255)
                                                                                                                           ttttCTTGTTTCCAGGCGCGAA‑AAAAGTTAACGACTATTTTTCTCATTTTGCAGTGTCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAt  <  2:565033/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                 
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAT  >  NC_000913_3_pae_tpiA/1543879‑1544149

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