Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA1,544,0211.A69.0% 89.4 / 39.7 42coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (11/2);  new base A (16/13);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 8.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGC  >  NC_000913_3_vch_tpiA/1543879‑1544161
                                                                                                                                               |                                                                                                                                               
ttGAGATTAACGGTAACTGCGCTACATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTGTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaa                                                                                                                                         <  1:337399/149‑2 (MQ=255)
 tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaag                                                                                                                                        >  1:242931/1‑147 (MQ=255)
  gagaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTTTGGTTTCCAGTCTCGAaaaaaga                                                                                                                                       >  2:577956/1‑146 (MQ=255)
   agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGCCTGCTATCCAGAGATTTCTTCGAGCAATATGATTTCCTTTATTTTGGCTGTAGCCCGGCTTTTCTGGTTTCCAGTCGCGACAAAAGAt                                                                                                                                      >  2:505904/1‑149 (MQ=255)
    gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCTATTCGCGTaaaaagata                                                                                                                                     >  2:815038/1‑144 (MQ=255)
    gaTTAACGGTAACTGCGCTAAATAGCATCAAGACACAGAAGGCAGGAAGTAGACGACAGACAGATATCCAAAGATTTCTTCGAGCACTATGTTTTCCTTTATTGTCCAAGAAGGCCGGATTTTCTCGTTTCCAGTCGCGAAAAAggata                                                                                                                                    >  1:703587/1‑145 (MQ=255)
      ttAACGGTAACTGCGCTAATTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGTAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAAc                                                                                                                                   <  2:606348/149‑1 (MQ=255)
      ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATACCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAAc                                                                                                                                   <  1:546495/149‑1 (MQ=255)
       tAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACg                                                                                                                                  >  1:212362/1‑149 (MQ=255)
                  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTTCGGTTGTGTATAAAGATAACGACTAttttttt                                                                                                                       >  1:361043/1‑147 (MQ=255)
                  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtct                                                                                                                       >  1:362141/1‑149 (MQ=255)
                  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTAGTTTCCAGTCGCGAAAAAAGATAACGAATATTTTtct                                                                                                                       >  1:572998/1‑149 (MQ=255)
                           aGCATCCAGACACAGAAGGAAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                              <  1:226629/149‑1 (MQ=255)
                             cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCATTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                            <  1:121951/149‑1 (MQ=255)
                              aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAAAAGGCCGGATTTTCTGGTTTCCAATCGCGAAAAAAGATAACGACTCTTTTTCTCATTTTGCatag                                                                                                           >  2:21903/1‑146 (MQ=255)
                                      acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                   >  1:746485/1‑149 (MQ=255)
                                        aGAAGGCAAGAAGTAGACAACTGGCTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                 >  1:331872/1‑149 (MQ=255)
                                        aGAACGCAAGAAGTAGACAACTGACTTATCACCAGAGCTGTCTTCGATCAATATGTTTTCCTTTATTTTCCCAGACGGCCGCTTTTTCTGGTTTCCAGTCGCGAAAAAAGATCACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                 <  2:507421/149‑1 (MQ=255)
                                           tcgCAAAAAGTACACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGACGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                              <  2:534670/147‑1 (MQ=255)
                                                agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTATCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                         >  1:646326/1‑149 (MQ=255)
                                                      aGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTATCATTTTGCAGTGCCATATTGTTCCTCACATGcaaac                                                                                   >  2:515730/1‑149 (MQ=255)
                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                               <  1:566520/149‑1 (MQ=255)
                                                          aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                               >  1:284324/1‑149 (MQ=255)
                                                                 tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTGTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCCTATTGTTCCTCACATTCAAACATTGGTAATGa                                                                        >  1:729594/1‑149 (MQ=255)
                                                                 tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGTTAATGa                                                                        >  1:450815/1‑149 (MQ=255)
                                                                        cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                 >  1:462564/1‑149 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATCCACACATTGTTAATGAAAAAAAGa                                                                >  2:278984/1‑149 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGCCTATTTTTCTCCTTTTGCAGTGCCATATTGTTCGTCCCATGCACTCTTTGGTAATGAAAAAAAGACACAACACg                                                       >  2:444463/1‑149 (MQ=255)
                                                                                   tcGAGAAATATGTTTCCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                      <  1:680287/149‑1 (MQ=255)
                                                                                     gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATGGGTAATGACACAAAGACAAAACACGAgg                                                    >  1:68033/1‑149 (MQ=255)
                                                                                         aaTATGTTTTCCTTTATTCTCCCAGAAGGCGGGATTTTCTAGTTTCCAGTCGCGAAAAAAGATAACGACTCTTTTTCTCATTTTGCAGAGCCACATTGTTCCCCACATGCACCCATAGGTAATCAAAAAAAGACAAAACACGAGCTAAg                                                >  2:580146/1‑149 (MQ=255)
                                                                                           tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTGTGGTTTCCAGTCGCGAAAAATGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCGCACATGCACACATTGGTAATGAAAAAAAACCAAAACACGAGGTAAGgc                                              >  1:489221/1‑149 (MQ=255)
                                                                                                          tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                          <  2:323833/138‑1 (MQ=255)
                                                                                                          tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                          >  1:323833/1‑138 (MQ=255)
                                                                                                           ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                              <  1:624271/149‑1 (MQ=255)
                                                                                                               cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                             >  2:207875/1‑114 (MQ=255)
                                                                                                               cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                             <  1:207875/114‑1 (MQ=255)
                                                                                                                 gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTTTTTTTTTTTAGTGCCATATTGTTCCTTACATGTACACATTGGTAATGAAAAAAAAACAAAATAGGTTGTACGTCGCAATAAGAAGTTATTAGAAtt                        >  2:324366/1‑149 (MQ=255)
                                                                                                                   cggCCGGAGTTTCTCGATTCAAGACGAGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATGGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                      <  2:451608/148‑1 (MQ=255)
                                                                                                                    ggCGGGATTTTCGGGTTGCAAGTCGCGAAAAAAGATAACAACTATTTTTCTCATTTTGAAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                     <  2:203282/149‑1 (MQ=255)
                                                                                                                      ccGGATTTTCGGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                   <  1:50586/149‑1 (MQ=255)
                                                                                                                                       cAGTCGCGAAAAAAGATAACGACTATTTTTTTTATTTTTCTTTGCCATATTGTTCCTCACATGAACACCTTGGTAATGAAAAAAAGACAAAACACGACGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATAAATTGAGTGAAATGc  >  2:233588/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGC  >  NC_000913_3_vch_tpiA/1543879‑1544161

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