Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 1,544,022 C→A L13F (TTG→TTT narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA1,544,0220CA96.2% 70.9 / NA 26L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (13/12);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_vch_tpiA/1543897‑1544157
                                                                                                                              |                                                                                                                                          
gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCATAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtct                                                                                                                   <  1:73178/149‑1 (MQ=255)
  gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCAAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCATTCGCGAAAAAAGATAACGACTATTTTTCTCa                                                                                                                 >  2:249514/1‑149 (MQ=255)
        tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATAGCCAGACATTACTCCGAGCAATATGTGTTCCTTTATTTTCCCCGAAGGCCGGATTTTCGCGTTTCCAGNCGCGAAATAANANANAGACTATTTTTCTCATATTGc                                                                                                           >  1:485284/1‑149 (MQ=255)
                        caGGCAANAAGTAGACAGCTGCTTAATGTCCAGAGATTTCTTCGATCAATATGTATACCTTTATTTTCCCTGTAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                           <  2:485284/148‑1 (MQ=255)
                         aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                          >  1:83605/1‑149 (MQ=255)
                              agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                     <  2:83605/149‑1 (MQ=255)
                                       cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTGCAGTGCCATATTGTTc                                                                                            >  1:431296/1‑132 (MQ=255)
                                       cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTGCAGTGCCATATTGTTc                                                                                            <  2:431296/132‑1 (MQ=255)
                                            gactgaTATCCAGAGATTTGTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                       <  1:426770/149‑1 (MQ=255)
                                                gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCACATGCACACATTGGTAATGaa                                                                   <  1:7744/149‑1 (MQ=255)
                                                         agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTACCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                          >  1:296628/1‑149 (MQ=255)
                                                               cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                    <  1:490457/149‑1 (MQ=255)
                                                                ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                   >  2:457629/1‑149 (MQ=255)
                                                                 tctaGCAATATGTTTTGCTTTGTTTTCCCATAATGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTTCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                  <  2:146451/146‑1 (MQ=255)
                                                                         tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGGGTAAGgc                                          >  1:51214/1‑149 (MQ=255)
                                                                                  cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcgccacaacc                                 >  2:551773/1‑142 (MQ=255)
                                                                                   tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                >  1:326806/1‑149 (MQ=255)
                                                                                   tttATTTTCCCAGAAGGCCGGATTTTCTGGATTCCAGACGCGACACAAGATAACAACTCTTTCACTCCTCTTGCAGCGCCAAATTCTTCCTCACATAAAAACCTAGGTCATGAACAAAAGACAACACCCGACGTATGGCGGAATCGccc                                >  2:13806/1‑148 (MQ=255)
                                                                                   ttaagtttCACAGAAGGGCGGATTTTTTGGTTTTCAGGCGCGAAATAAGATAACGAGTATTTTTCTAATTTTGCACCGCCATATTGTTCCTGACATGTACATATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                <  2:28876/144‑1 (MQ=255)
                                                                                             cAGACGGCAGGATGTTCTGTTTTTCAGCTACGAAAAAAGATAACGATGATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACCCGCGGTAAGGCGCAATAGCCAGTTATTAGaa                      <  2:218353/149‑1 (MQ=255)
                                                                                             cAGAATGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                      <  2:390182/149‑1 (MQ=255)
                                                                                              aGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                     >  1:555775/1‑149 (MQ=255)
                                                                                                     cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTACCCTAATGTTCCTCACATGCACACATTGGTAATGAAAAAAAGATAAAACACGGGGTAAGGGGCAATAGCCAGTTATTAGAATTAAGGAt              >  1:179445/1‑149 (MQ=255)
                                                                                                     cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAAAGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGAAATAGCCAGTTATTAGAATTAAGGAt              >  2:560872/1‑149 (MQ=255)
                                                                                                       gatgtggtGGTTTCCGGCCGCGAATAAAGATAACGACTTTTTTTCGTGTTGTGCATTGCGATATTGTTCCGCAAAGACACACATTTTTAATGAAAAAAAGACAAAACACGAGTTAAGGCTCAATAGCCAGTTATTAGAATTAAGGACGa            <  2:212765/142‑1 (MQ=255)
                                                                                                                 tttCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCATCACCTACACTCATTGGTAATGAAAAACAGACAAAACACGAGATAAGGCGCAATAGACAGTTATTAAAATAAAGAATCAATTGGATGaa  >  2:572656/1‑149 (MQ=255)
                                                                                                                              |                                                                                                                                          
GCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAA  >  NC_000913_3_vch_tpiA/1543897‑1544157

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: