Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 1,544,022 1 bp→AA coding (39/1389 nt) narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,544,0211.A79.0% 40.2 / 8.0 19coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (9/6);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01
*NC_000913_3_pae_tpiA1,544,0220CA94.7% 46.6 / NA 19L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (11/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATG  >  NC_000913_3_pae_tpiA/1543896‑1544146
                                                                                                                               |||                                                                                                                               
tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATACCGACTATTTTtc                                                                                                          >  2:76591/1‑149 (MQ=255)
         tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCAtttttc                                                                                                 >  2:236207/1‑147 (MQ=255)
              tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                            >  1:503531/1‑149 (MQ=255)
                cAGAAACAGAAGGTAAGAAGTAAACAACTGGCTGATATCCAGAGATTTCTTCGAGTAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAATAACGCATAAAGATAACGACTCTGTTTCTCATTTTGCAGTGCCa                                                                                          <  2:88812/149‑1 (MQ=255)
                               agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCGGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACAt                                                                           >  2:403124/1‑149 (MQ=255)
                                  agcaGACAACCGACTGATATCCAGCGATTTCTTCGACTAATATGTGTTCCTTTATTTTCTCAGCAGGCCGGATTTTCCGGTTTCCAGTCGCGAAAAGAGATAACCACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCTCCCGTGct                                                                        <  2:482100/146‑2 (MQ=255)
                                         aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTACTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                 <  1:434276/149‑1 (MQ=255)
                                            tgactgaTATCAAGATGGTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCCCCCATTGGTa                                                              >  2:324598/1‑149 (MQ=255)
                                                           gaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                              <  1:245115/149‑1 (MQ=255)
                                                                 gttgAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTTCAGTGCCATGTTGTTCCTCACAGGCACACATCGGTAATGCAAAAAACACAAAACACg                                        >  2:138321/4‑149 (MQ=255)
                                                                 ctacagaaaTAGGTTTTCCTTTATTTTCCCAGAAGGCCGAATTTTCTGGTTTGCAGTCGCGAAAAAGGATAACGACTATTTTTCTCATTTTGCAGTGCCATCTTGTTCCTCACATTCACACATTGTTATGAGAAAAAAGACTAAACACg                                        >  2:114638/8‑149 (MQ=255)
                                                                                   cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCAATATTGTTCCTCACATGCACACATTGGTAGTGAAAAAACAACAAAAAACGAGGTAAGGCGCAATGGcc                      >  2:360296/1‑149 (MQ=255)
                                                                                    tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                            <  2:442770/126‑1 (MQ=255)
                                                                                    tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                            >  1:442770/1‑126 (MQ=255)
                                                                                    tttATTTTCCCAAGAGCCTGGATTTTCGGGTTTCCAGTCGCGAAAAAAGATAACGACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                     <  1:68960/149‑1 (MQ=255)
                                                                                           tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta              <  1:180524/149‑1 (MQ=255)
                                                                                                  aGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGGCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAGGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa       <  1:82073/149‑1 (MQ=255)
                                                                                                    gCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg     >  1:330921/1‑149 (MQ=255)
                                                                                                       ggATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg  >  2:442548/1‑149 (MQ=255)
                                                                                                                               |||                                                                                                                               
TGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATG  >  NC_000913_3_pae_tpiA/1543896‑1544146

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: