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.A86.7% 79.7 / 7.3 30coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (15/11);  total (17/13)
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
*NC_000913_3_pae_tpiA1,544,0220CA96.7% 86.2 / NA 30L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (17/12);  minor base G (0/1);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAT  >  NC_000913_3_pae_tpiA/1543882‑1544165
                                                                                                                                              |||                                                                                                                                              
agaTTAACGGTAACTGCGCTGAATAGCATCCAGACACAGAAGGAAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑aaaaagat                                                                                                                                         <  2:291912/149‑5 (MQ=255)
  aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATaa                                                                                                                                        >  2:52390/1‑149 (MQ=255)
  aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑A‑aaaaagaaag                                                                                                                                       >  1:37167/1‑143 (MQ=255)
         gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGTTAAAGACTAt                                                                                                                                 >  1:34933/1‑149 (MQ=255)
             cTGCGCTAAACAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTGCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTAttttt                                                                                                                             <  2:20082/149‑1 (MQ=255)
                   tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCAtt                                                                                                                       >  2:496336/1‑149 (MQ=255)
                                    caGAAGGGAAAAATTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTACTCATTTTGCAGTGCCATATTGt                                                                                                      <  2:366461/149‑1 (MQ=255)
                                        aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                  <  1:496336/149‑1 (MQ=255)
                                          gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                >  2:494328/1‑149 (MQ=255)
                                           cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                               >  1:577549/1‑149 (MQ=255)
                                               aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                           >  2:332890/1‑149 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGCATTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACAGATTg                                                                                   >  1:151978/1‑149 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                   <  1:52390/149‑1 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                   >  1:131387/1‑149 (MQ=255)
                                                       aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCAGTATTGTTCCTCACATGCACACCTTg                                                                                   >  2:416984/1‑149 (MQ=255)
                                                          tgacttaTTTTCAGCGATTGTTTTGAGTAATAGGTTGTCCTTTAATCGCCCGGAAGGCCTGATTTTCAGTGTGTCAGGCG‑GGAG‑AAAAGATAACGACTATGTTTCTGATGTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                <  2:127457/149‑1 (MQ=255)
                                                            actgaGATCAAGAGATTTCATCGAGCAATATGATTTACTTAATTTTAACAGAAGGCCGGATTTTCTGGTAAAAAAAAA‑AAAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                              <  2:181780/149‑1 (MQ=255)
                                                                    ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                      >  1:330024/1‑149 (MQ=255)
                                                                          aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                                                >  2:283323/1‑149 (MQ=255)
                                                                              ctccGAGCAATATGTTTTCCTTTAGTTTCCCAGAAGGCCGGAGTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACGATGTTGCTCATTTTGCAGTGCCATAGTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                            <  1:265034/146‑1 (MQ=255)
                                                                              attcGAGCAATCTGTTTTCCTTTATGTTTTCAGAAGGCTGGATTGTCAAGTTGCCATTCGCCA‑A‑AAAAGATAACGACTATTTTTCTCTTTCTGCAGTGCTATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                            <  2:551878/148‑1 (MQ=255)
                                                                               ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCATCACCTGCACACATTGGTAATGAAAAAAAGACAAAacact                                                           >  2:346082/1‑148 (MQ=255)
                                                                                  gaccAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                        <  1:204996/146‑1 (MQ=255)
                                                                                      aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTTTTTTTTTTTTTTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                    >  2:549482/1‑149 (MQ=255)
                                                                                                                 ggCTGGATTTTCTGGTTTCTAGTTG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                         <  2:37167/149‑1 (MQ=255)
                                                                                                                   ccGGATTTTCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                       >  1:321139/1‑149 (MQ=255)
                                                                                                                           tCTGGTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTGCCTCACATGCACACATTGGTAATGACAAAAAGACAAAACACGAGGTAAGGCGCATAAGCCAGTTATTAGAATAAAGGATAAATTgg               >  2:631150/1‑149 (MQ=255)
                                                                                                                              ggTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                             >  2:136151/1‑68 (MQ=255)
                                                                                                                              ggTTTCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                             <  1:136151/68‑1 (MQ=255)
                                                                                                                                  tCCAGTCG‑CGAA‑AAAAGATAACGACTATTTTTCTAATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGt        <  2:131387/149‑1 (MQ=255)
                                                                                                                                          cgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt  <  1:332890/143‑1 (MQ=255)
                                                                                                                                              |||                                                                                                                                              
AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCG‑CG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAT  >  NC_000913_3_pae_tpiA/1543882‑1544165

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