Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 1,544,022 1 bp→AA 100% 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.A90.6% 94.5 / 4.7 32coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (14/15);  total (16/16)
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.51e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 100.5 / NA 32L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/16);  total (16/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_pae_tpiA/1543895‑1544169
                                                                                                                                |||                                                                                                                                                   
cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTAttttt                                                                                                                                 <  1:926023/149‑1 (MQ=255)
 tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTtc                                                                                                                                >  2:593495/1‑149 (MQ=255)
  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTtct                                                                                                                               <  1:304142/149‑1 (MQ=255)
   cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTtctc                                                                                                                              >  1:290305/1‑149 (MQ=255)
                        agaagaCAAGAAGTAGACAAGTGACTGATATCAAGAGATTTCTTCAAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTACAAAAAAAAAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                         <  2:660607/149‑1 (MQ=255)
                         gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTTAAGTGCCATATTGTTc                                                                                                        >  2:23894/1‑149 (MQ=255)
                             gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                    >  1:704964/1‑149 (MQ=255)
                                       gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAGTTTCCCAGAAGGCCGGATGTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTACCATCTTGTTCCTCACATGccaaca                                                                                          >  2:85300/1‑144 (MQ=255)
                                           actgactgATATCCAGAGATTTCTTNGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                      <  1:621923/149‑1 (MQ=255)
                                               actgaTATTCAGTGTGTTCTCAGTGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                  <  2:96391/149‑1 (MQ=255)
                                                    tatCCAGAGATTTCTGCGAGCAAGATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTTACATGCCCACATTGGTAATGaaaa                                                                             >  2:87158/1‑149 (MQ=255)
                                                      tcaagaaATTTCTTCGAGCAATATGTTTTCCTTTAGTTCCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAA‑GATAACGACTATTTTTCTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                           <  2:77197/146‑1 (MQ=255)
                                                       ccAGATATTTCTTCGATCAATATGTTTTCCTTTCTTTTCCCAGAAGGCCCGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTCTTTTTCTCATTTTGCAGTGACATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                          <  1:440324/149‑1 (MQ=255)
                                                        cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                         >  2:166559/1‑149 (MQ=255)
                                                         agagATCTCCTCGAACAATACGTTTTCCTGTATTTTCCCAGAAGGCCGGTTTTTCTGGTTAAAAATAAAAAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                        <  2:519061/149‑1 (MQ=255)
                                                          gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCAAGTAGCG‑A‑AAAACGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCATCACATGCACACATTGGTAATGAAAAAAAGAc                                                                       >  2:69844/1‑149 (MQ=255)
                                                           agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                      >  1:274010/1‑149 (MQ=255)
                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACAACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGCCAAAACACg                                                               >  2:247880/1‑149 (MQ=255)
                                                                   tccaGCAATATATTTTCCTTTATTTTCCCAAAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                              <  2:306519/146‑1 (MQ=255)
                                                                   gtgAGTAATATTTTTTCCTTTATTTTCTCATAAGGCGGGATGTTCTGGTTGCCAGGCGCGAA‑AAAA‑GATAACGACTGTTTTTCTCATGTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                              <  1:322232/147‑1 (MQ=255)
                                                                      aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                           >  1:24909/1‑149 (MQ=255)
                                                                        cAATTTGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                         <  2:587621/149‑1 (MQ=255)
                                                                           tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACAAGGAAAGgc                                                      >  2:382832/1‑149 (MQ=255)
                                                                               ttttTCTTTATTGTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                  <  2:582313/149‑1 (MQ=255)
                                                                                  tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAg                                               >  2:468734/1‑149 (MQ=255)
                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATGGCCa                                            >  1:765235/1‑149 (MQ=255)
                                                                                            tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                     <  1:846108/149‑1 (MQ=255)
                                                                                                       cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                          <  1:593495/149‑1 (MQ=255)
                                                                                                       cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCCCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                          >  1:828301/1‑149 (MQ=255)
                                                                                                             ttCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                    <  2:304541/149‑1 (MQ=255)
                                                                                                              tCTGGTTTCCAGTCGTGAA‑AAAA‑GATAACCACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                   <  2:116277/149‑1 (MQ=255)
                                                                                                                  gTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa               >  2:789501/1‑149 (MQ=255)
                                                                                                                            cgcgaaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt      <  1:192965/143‑1 (MQ=255)
                                                                                                                               gaaaaac‑gaaaacgaGAATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGACAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATAGa   >  2:810834/4‑149 (MQ=255)
                                                                                                                                aaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa  >  1:578118/3‑149 (MQ=255)
                                                                                                                                |||                                                                                                                                                   
CTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_pae_tpiA/1543895‑1544169

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