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% 139.2 / 4.8 42coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base A (18/21);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 147.5 / NA 42L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/21);  total (21/21)

AACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543887‑1544170
                                                                                                                                        |||                                                                                                                                                    
aaCGGTAACTACGCGAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGCCTGATATCCAGAGATGTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGAc                                                                                                                                         >  4:223766/1‑150 (MQ=255)
          gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTtctc                                                                                                                               >  1:509868/1‑150 (MQ=255)
          gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTAAAGAAGAGAAAAAAGATAACGACTATTTTtctc                                                                                                                               <  2:587233/150‑1 (MQ=255)
                         cAGACACAGAATGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTACCAGAAGGCCGGATTTTCTGGTTAAAAAAAAAAAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                <  4:210005/150‑1 (MQ=255)
                                aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  1:579599/150‑1 (MQ=255)
                                aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTATTc                                                                                                         >  3:97584/1‑150 (MQ=255)
                                  aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                       >  1:480102/1‑150 (MQ=255)
                                   aGGCAAGAAGTAGACAACTGACTGATATCCAGAAATTGCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATGGTTCGCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      <  3:182592/150‑1 (MQ=255)
                                      cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCCTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                   >  1:366719/1‑150 (MQ=255)
                                         gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTGTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                                <  2:329164/150‑1 (MQ=255)
                                             tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                            >  4:59592/1‑150 (MQ=255)
                                                    ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  1:280624/150‑1 (MQ=255)
                                                    ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  4:206935/150‑1 (MQ=255)
                                                          gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa                                                                               <  2:392425/150‑1 (MQ=255)
                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAt                                                                                                                                      >  4:81115/1‑90 (MQ=255)
                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAt                                                                                                                                      <  3:81115/90‑1 (MQ=255)
                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  2:349272/1‑150 (MQ=255)
                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  3:212896/1‑150 (MQ=255)
                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  3:161778/1‑150 (MQ=255)
                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         <  1:604306/150‑1 (MQ=255)
                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGACGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         <  4:131914/150‑1 (MQ=255)
                                                                cAGAGATTTCCTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  3:161765/1‑150 (MQ=255)
                                                                     aTTTCTTTGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGACTGATTTACAGTTCACCAGTCGAGAAAAAACATACCGACAATTAATCTCAACTTGCAGTGCCATATTGTTCCTCACCTACACAAAAAAAAAAAAAAAAAAAGACaaaa                                                                    <  4:148173/150‑1 (MQ=255)
                                                                     aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    <  3:59592/150‑1 (MQ=255)
                                                                     aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGACGAATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAATAAAAATAAAAAAAAGACaaaa                                                                    <  2:497593/150‑1 (MQ=255)
                                                                         cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                         <  1:431699/93‑1 (MQ=255)
                                                                         cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                         >  2:431699/1‑93 (MQ=255)
                                                                            cGAGCAATATGTTCTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  4:326580/150‑1 (MQ=255)
                                                                                  atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               <  4:338730/94‑1 (MQ=255)
                                                                                  atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  3:338730/1‑94 (MQ=255)
                                                                                       ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATCAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                                  >  4:212131/1‑150 (MQ=255)
                                                                                            cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  1:115314/150‑1 (MQ=255)
                                                                                            cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  3:41682/1‑150 (MQ=255)
                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  1:37303/1‑150 (MQ=255)
                                                                                               tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAACACAGATAACGACTATTTTTCTCAATTTGCAGTACCATATTGTTCCTCACATGCACCCATTGGTAATGAAACAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                          >  2:285109/1‑150 (MQ=255)
                                                                                                 ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                        <  3:329084/150‑1 (MQ=255)
                                                                                                  tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       >  4:102639/1‑150 (MQ=255)
                                                                                                          aaGGCCGGATTTTCTGGTTTACAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACAAAAAAAAAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  4:246726/150‑1 (MQ=255)
                                                                                                               cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                         <  3:202451/55‑1 (MQ=255)
                                                                                                               cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                         >  4:202451/1‑55 (MQ=255)
                                                                                                                      tCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   >  1:339348/1‑150 (MQ=255)
                                                                                                                           tttCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg              >  2:445800/1‑150 (MQ=255)
                                                                                                                                   tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  3:187360/143‑1 (MQ=255)
                                                                                                                                     gcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg     >  3:146227/6‑150 (MQ=255)
                                                                                                                                     gcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg     <  4:97584/145‑1 (MQ=255)
                                                                                                                                        aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:397006/148‑1 (MQ=255)
                                                                                                                                        aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  2:398876/3‑150 (MQ=255)
                                                                                                                                        |||                                                                                                                                                    
AACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543887‑1544170

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