Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 1,544,021:1 +A 78.1% coding (40/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.A78.1% 74.4 / 18.7 32coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base A (11/13);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_pae_tpiA/1543874‑1544158
                                                                                                                                                     |                                                                                                                                           
ttgtATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATGTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTc‑gcga                                                                                                                                             <  2:72740/146‑2 (MQ=255)
ttgtATTGAGATTAACGGTAACGGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATAGCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTc‑gcga                                                                                                                                             <  2:72833/146‑2 (MQ=255)
     ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAAGCAAGAAGGAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATTTTTTCCTTTATTTTCCCACAAGGCCGGATTTTCTGGTTTCCAGTC‑CCG‑A‑aaa‑aa                                                                                                                                     >  2:557105/1‑148 (MQ=255)
                      tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTtc                                                                                                                     >  1:506092/1‑149 (MQ=255)
                         gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑A‑AAACAGATAACGACTATTTTTCTCa                                                                                                                  >  1:666407/1‑149 (MQ=255)
                               tAGCATCAAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGCTGACTTCGAGCAAGAAGTTTTCCTTTATTTTCCCATAAGGCCGGATTTTCTGGTTTCCATTG‑GCGTA‑AAA‑AGATAACGACTATTTTTCTCATCTTGc                                                                                                            >  2:708838/1‑149 (MQ=255)
                                 gtgtCCAGACGCAGATGTCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTACTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAg                                                                                                          <  2:281470/146‑1 (MQ=255)
                                    tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                       >  2:603081/1‑149 (MQ=255)
                                      cAGACACAGAGGGGAAGAAGGAGACAACGGACGGATCTCCAGAGATTTCAGCGAGCAATATGTTTTCCTTTATTTTACCAGAAGGCCGGATTTTCTGGATTCCAATC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                     <  2:282663/149‑1 (MQ=255)
                                       aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCATTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTACAGTC‑GCGAA‑AAA‑AGATAACCACTATTCTTCTATTTTTGCAGTGCCat                                                                                                    >  2:426551/1‑149 (MQ=255)
                                            caGAAGGCAAGAAGGAGACAACTGACGGATAACCAGAGATTTCTTCGAGCAATATGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                               <  1:113996/149‑1 (MQ=255)
                                                    aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                       <  1:133968/149‑1 (MQ=255)
                                                                      tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGCTTTTCTGGTTTCCAGTC‑TCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGACATATTGTTCCTCACATGCACACTATGGTAATGa                                                                     >  2:162538/1‑149 (MQ=255)
                                                                             cAGAGATTTCTTCGAGCAATATGTTTTCCTTTCTTTTCCCGGAAGGCAGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGTATATTTTTCTCATTTTGCAGTGCACTATTGTTACTCTCATGCACACATATGTGATAaaaataat                                                              >  2:551858/1‑148 (MQ=255)
                                                                             aagagATTTCTCCAAGCAACATGCCTTCCCTTCTTTACCCAGAAGACCCCATTATCTAGTTACCACAC‑AAGAA‑AAA‑AGAGAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACCTGCACACATTGGTAATGAAAAAAAg                                                              <  2:528636/148‑1 (MQ=255)
                                                                              agagCTTTCCTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCATATGCACACATTGGTAATGAAAAAAAGa                                                             <  1:87248/149‑1 (MQ=255)
                                                                               gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                            >  1:135646/1‑149 (MQ=255)
                                                                                        tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                   <  1:494769/149‑1 (MQ=255)
                                                                                          gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                 >  1:129884/1‑149 (MQ=255)
                                                                                          gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                 >  1:130389/1‑149 (MQ=255)
                                                                                             cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                              >  2:662880/1‑149 (MQ=255)
                                                                                                tatGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                           <  1:169927/149‑1 (MQ=255)
                                                                                                     tttCCTTTATTTTCTCGGAGGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCACATGCACACAGTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAAt                                      <  2:512371/149‑1 (MQ=255)
                                                                                                          tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGccg                                 >  1:55246/1‑148 (MQ=255)
                                                                                                             aTTGTCCGAGAAGGACGGTTTTTCTGGTTTCGAGTC‑GGGAA‑AAA‑AGATAACGACGATTGTTCTGATTTTGCAGTGCCATATGGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                              <  2:116081/149‑1 (MQ=255)
                                                                                                                          gCCGGGTTTTCTGGTTTCCAGTCGGTA‑A‑AAA‑AGATAACGACTATTTTTCTCATTTTGTAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg                 <  2:61135/149‑1 (MQ=255)
                                                                                                                           taggttttcCTGTTTTCCAGTCTACA‑A‑AAA‑AGACAACGACTATGTTTCGAAATTTGCAGTGCCATTTTGTTCCTCACATGCACACATTGGTAAGGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                <  2:291915/144‑1 (MQ=255)
                                                                                                                             ggATTTTCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaaaaaag              >  1:145301/1‑144 (MQ=255)
                                                                                                                             ggATTTTCCGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTCGAATTAAGGagg              >  1:91114/1‑147 (MQ=255)
                                                                                                                                  ttCTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg         <  1:588047/149‑1 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg       <  1:475788/149‑1 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAAGAATTggg       >  2:293888/1‑149 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTA‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg       <  2:135646/149‑1 (MQ=255)
                                                                                                                                         ttCCAGTC‑GCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGATATTAGAATTAAGGATGAATTGGGTGAAg  <  1:603081/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_pae_tpiA/1543874‑1544158

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: ‑