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.A92.9% 38.6 / ‑2.4 14coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (6/7);  total (7/7)
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.39e-01
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 41.2 / NA 14L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (7/7);  total (7/7)

TAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTA  >  NC_000913_3_pae_tpiA/1543905‑1544140
                                                                                                                      |||                                                                                                                     
tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAtt                                                                                                       <  1:210220/136‑1 (MQ=255)
        aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCAt                                                                                               >  1:192100/1‑136 (MQ=255)
              aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                         <  1:185519/136‑1 (MQ=255)
                                      ctgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATTca                                                                 >  2:170782/1‑136 (MQ=255)
                                                 agaTTTTTTCGACAAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGGTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCCTTTTGCAGTACCTGCGTGTTCCTCACATGCCCACATTGGTaa                                                      >  2:179634/1‑136 (MQ=255)
                                                   aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCAGATTGTTCCTCACATGCACACATTGGTACTg                                                    >  2:193402/1‑136 (MQ=255)
                                                     ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                  >  2:7720/1‑136 (MQ=255)
                                                           gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                            <  1:7720/136‑1 (MQ=255)
                                                                                    cgagagGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                   <  1:183714/134‑1 (MQ=255)
                                                                                         aGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                       <  1:199196/79‑1 (MQ=255)
                                                                                         aGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                       >  2:199196/1‑79 (MQ=255)
                                                                                                   ttCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtc                                                                                                   >  1:77385/1‑41 (MQ=255)
                                                                                                   ttCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtc                                                                                                   <  2:77385/41‑1 (MQ=255)
                                                                                                     ttGGTTTTCAGTCGCGAAAAAAGATAAGGACTATTTTTCTCTTTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTa  <  1:170782/135‑1 (MQ=255)
                                                                                                                      |||                                                                                                                     
TAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTA  >  NC_000913_3_pae_tpiA/1543905‑1544140

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