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.A75.0% 102.7 / 33.3 40coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (5/5);  new base A (15/15);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01
*NC_000913_3_pae_tpiA1,544,0220CA97.5% 126.5 / NA 40L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (20/19);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543898‑1544170
                                                                                                                             |||                                                                                                                                                    
cgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtctc                                                                                                                               >  1:437049/1‑149 (MQ=255)
        aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCAtt                                                                                                                            <  2:336321/144‑1 (MQ=255)
        aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:527569/149‑1 (MQ=255)
        aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCAtt                                                                                                                            >  1:336321/1‑144 (MQ=255)
            tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   <  2:437049/149‑1 (MQ=255)
                     aGAAGGCAAGAAGTAGAAAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                          <  1:281035/149‑1 (MQ=255)
                            aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTCTTCATcaca                                                                                                   >  2:483240/1‑149 (MQ=255)
                               aagTAGACAACTTACATATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTCTCTACTCTCCAATCCCG‑CAAAA‑ACACAACACCCACATCTCCCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                                <  2:769569/149‑1 (MQ=255)
                                 gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATCTGGTTCCTCACATGcac                                                                                              >  1:574920/1‑149 (MQ=255)
                                 gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaa                                                                                              >  2:504079/1‑148 (MQ=255)
                                     aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGATGATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  2:558867/138‑1 (MQ=255)
                                     aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     >  1:558867/1‑138 (MQ=255)
                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        >  2:59619/1‑149 (MQ=255)
                                               gaTATACAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACCTTAACACATTGGTAATGaa                                                                                >  2:268436/1‑149 (MQ=255)
                                                  atCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGACGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCAAATGCACAGATGGGTCATGaaaaa                                                                             >  1:297047/1‑149 (MQ=255)
                                                    ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGGTTTCCAGTAGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTAGTAATGaaaaaaa                                                                           <  2:596704/149‑1 (MQ=255)
                                                         gaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCTGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      <  1:592274/149‑1 (MQ=255)
                                                          aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                                                     >  2:664712/1‑149 (MQ=255)
                                                            ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aaga                                                                                                                                               >  2:405323/1‑69 (MQ=255)
                                                            ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aaga                                                                                                                                               <  1:405323/72‑4 (MQ=255)
                                                             tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAaca                                                                  >  1:270220/1‑149 (MQ=255)
                                                              cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          >  1:723474/1‑140 (MQ=255)
                                                              cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          <  2:723474/140‑1 (MQ=255)
                                                                tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               >  1:715913/1‑149 (MQ=255)
                                                                tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  2:645399/149‑1 (MQ=255)
                                                                   aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                            >  2:90394/1‑149 (MQ=255)
                                                                             tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                                  >  1:757146/1‑149 (MQ=255)
                                                                                  tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGGCAGGCGCCAtgtgca                                             >  2:741537/1‑144 (MQ=255)
                                                                                   ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  1:90394/149‑1 (MQ=255)
                                                                                         tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttctta                                      >  2:756358/1‑149 (MQ=255)
                                                                                             ggAAGGGCGGATTTTGTGGTTGCCTGTCGCGAA‑AAA‑AGATAACGTCTATTTTTCTCATTTTGCAGTGCCATCTTGTTCTTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                  <  2:641617/148‑1 (MQ=255)
                                                                                               gaGGTCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTa                                <  2:95174/148‑1 (MQ=255)
                                                                                                aGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  1:386086/149‑1 (MQ=255)
                                                                                                         tttCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:540456/140‑1 (MQ=255)
                                                                                                         tttCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               >  1:540456/1‑140 (MQ=255)
                                                                                                           tCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:664712/149‑1 (MQ=255)
                                                                                                           tCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  2:757146/149‑1 (MQ=255)
                                                                                                            cgggTTTCCTGTGGCG‑A‑AAATAGATGACGTCTATTTTTGTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGATTTAAGGATGAATTggg                   <  2:316753/147‑1 (MQ=255)
                                                                                                            cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  1:518790/149‑1 (MQ=255)
                                                                                                               gTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa                >  1:439660/1‑149 (MQ=255)
                                                                                                                        tcgcgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        <  1:309372/142‑1 (MQ=255)
                                                                                                                        tcgcgaaaaa‑aGATAACGACAATTTGTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        <  2:176779/142‑1 (MQ=255)
                                                                                                                         cgcgaaaaa‑aGATAACGATTATGTTTCTCATTTTGAAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  2:574920/143‑1 (MQ=255)
                                                                                                                           cgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       >  2:636955/5‑115 (MQ=255)
                                                                                                                           cgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       <  1:636955/111‑1 (MQ=255)
                                                                                                                           cgaaaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGAGAAGTGCTGAtgg     >  2:467911/5‑147 (MQ=255)
                                                                                                                             aaaaa‑aGATAACGAATACTTTTCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGAGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa   <  1:268436/147‑1 (MQ=255)
                                                                                                                              aaaa‑aTATAACGACCAATTTTCTCATTTTACAGTGCCAAATTGTTCCTCACATGCACACATTGGTAAAGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:696117/148‑1 (MQ=255)
                                                                                                                              aaaa‑aGATAACGACTATTTCTATAATTATGCTGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAGACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  2:673853/2‑149 (MQ=255)
                                                                                                                             |||                                                                                                                                                    
CGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543898‑1544170

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