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

TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  NC_000913_3_pae_tpiA/1543874‑1544158
                                                                                                                                                     |                                                                                                                                           
ttgtATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATGTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcga                                                                                                                                             <  2:72740/146‑2 (MQ=255)
ttgtATTGAGATTAACGGTAACGGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATAGCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcga                                                                                                                                             <  2:72833/146‑2 (MQ=255)
     ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAAGCAAGAAGGAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATTTTTTCCTTTATTTTCCCACAAGGCCGGATTTTCTGGTTTCCAGTCCCGAaaaaa                                                                                                                                     >  2:557105/1‑148 (MQ=255)
                      tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtc                                                                                                                     >  1:506092/1‑149 (MQ=255)
                         gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAACAGATAACGACTATTTTTCTCa                                                                                                                  >  1:666407/1‑149 (MQ=255)
                               tAGCATCAAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGCTGACTTCGAGCAAGAAGTTTTCCTTTATTTTCCCATAAGGCCGGATTTTCTGGTTTCCATTGGCGTAAAAAGATAACGACTATTTTTCTCATCTTGc                                                                                                            >  2:708838/1‑149 (MQ=255)
                                 gtgtCCAGACGCAGATGTCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTACTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAg                                                                                                          <  2:281470/146‑1 (MQ=255)
                                    tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                       >  2:603081/1‑149 (MQ=255)
                                      cAGACACAGAGGGGAAGAAGGAGACAACGGACGGATCTCCAGAGATTTCAGCGAGCAATATGTTTTCCTTTATTTTACCAGAAGGCCGGATTTTCTGGATTCCAATCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                     <  2:282663/149‑1 (MQ=255)
                                       aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCATTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTACAGTCGCGAAAAAAGATAACCACTATTCTTCTATTTTTGCAGTGCCat                                                                                                    >  2:426551/1‑149 (MQ=255)
                                            caGAAGGCAAGAAGGAGACAACTGACGGATAACCAGAGATTTCTTCGAGCAATATGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                               <  1:113996/149‑1 (MQ=255)
                                                    aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                       <  1:133968/149‑1 (MQ=255)
                                                                      tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGCTTTTCTGGTTTCCAGTCTCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGACATATTGTTCCTCACATGCACACTATGGTAATGa                                                                     >  2:162538/1‑149 (MQ=255)
                                                                             cAGAGATTTCTTCGAGCAATATGTTTTCCTTTCTTTTCCCGGAAGGCAGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGTATATTTTTCTCATTTTGCAGTGCACTATTGTTACTCTCATGCACACATATGTGATAaaaataat                                                              >  2:551858/1‑148 (MQ=255)
                                                                             aagagATTTCTCCAAGCAACATGCCTTCCCTTCTTTACCCAGAAGACCCCATTATCTAGTTACCACACAAGAAAAAAGAGAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACCTGCACACATTGGTAATGAAAAAAAg                                                              <  2:528636/148‑1 (MQ=255)
                                                                              agagCTTTCCTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCATATGCACACATTGGTAATGAAAAAAAGa                                                             <  1:87248/149‑1 (MQ=255)
                                                                               gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                            >  1:135646/1‑149 (MQ=255)
                                                                                        tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                   <  1:494769/149‑1 (MQ=255)
                                                                                          gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                 >  1:129884/1‑149 (MQ=255)
                                                                                          gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                 >  1:130389/1‑149 (MQ=255)
                                                                                             cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                              >  2:662880/1‑149 (MQ=255)
                                                                                                tatGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                           <  1:169927/149‑1 (MQ=255)
                                                                                                     tttCCTTTATTTTCTCGGAGGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCACATGCACACAGTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAAt                                      <  2:512371/149‑1 (MQ=255)
                                                                                                          tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGccg                                 >  1:55246/1‑148 (MQ=255)
                                                                                                             aTTGTCCGAGAAGGACGGTTTTTCTGGTTTCGAGTCGGGAAAAAAGATAACGACGATTGTTCTGATTTTGCAGTGCCATATGGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                              <  2:116081/149‑1 (MQ=255)
                                                                                                                          gCCGGGTTTTCTGGTTTCCAGTCGGTAAAAAAGATAACGACTATTTTTCTCATTTTGTAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg                 <  2:61135/149‑1 (MQ=255)
                                                                                                                           taggttttcCTGTTTTCCAGTCTACAAAAAAGACAACGACTATGTTTCGAAATTTGCAGTGCCATTTTGTTCCTCACATGCACACATTGGTAAGGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                <  2:291915/144‑1 (MQ=255)
                                                                                                                             ggATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaaaaaag              >  1:145301/1‑144 (MQ=255)
                                                                                                                             ggATTTTCCGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTCGAATTAAGGagg              >  1:91114/1‑147 (MQ=255)
                                                                                                                                  ttCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg         <  1:588047/149‑1 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg       <  1:475788/149‑1 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAAGAATTggg       >  2:293888/1‑149 (MQ=255)
                                                                                                                                    cTGGTTTCCAGTAGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg       <  2:135646/149‑1 (MQ=255)
                                                                                                                                         ttCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGATATTAGAATTAAGGATGAATTGGGTGAAg  <  1:603081/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
TTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAG  >  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: