Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_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_vch_tpiA1,544,0211.A82.3% 80.1 / 12.8 34coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base A (16/12);  total (19/15)
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_vch_tpiA1,544,0220CA100.0% 92.4 / NA 34L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/15);  total (19/15)

TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_vch_tpiA/1543877‑1544169
                                                                                                                                                  |||                                                                                                                                                    
tATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑a‑aa                                                                                                                                                   >  1:190046/1‑149 (MQ=255)
     agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAATGCAAAAAGTAAACAACTAACAGATAACCACAGATTTCTTCGAGCAATCTGTTTTCCTTTATTTTCCCAGAAGGCCGGAGTTTCTGGTTTCCAGTCGCG‑A‑a‑aa‑aagat                                                                                                                                             >  1:156714/1‑145 (MQ=255)
       aataacggtaacggCGTTAAATAGCATCCAGACACAGAAGGCAATAAGTAGACAACTGACTGATTTCCAGAGATTTCTTCTAGCAATATTTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTTGTTTCCAGTCGCG‑AAA‑AA‑AGATaa                                                                                                                                            <  2:103869/147‑1 (MQ=255)
       aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATaa                                                                                                                                            <  2:311372/149‑1 (MQ=255)
          aaCGGGAAGCGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGa                                                                                                                                         <  2:170477/149‑1 (MQ=255)
                   tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTAGAGCAATATGTGTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCTCGAA‑A‑AA‑AGGTAACGACTATTTTtc                                                                                                                                >  2:452706/1‑149 (MQ=255)
                              gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCACAAAATTCTTCGAGCAATATGCTTTTCTCTATTTTCCCAGCATCCCGGATTTTCTGCTTTCCAATCCCGAA‑A‑AA‑AGAAACCATCTTTTTTTCTGAtttctaat                                                                                                                     >  2:164421/1‑144 (MQ=255)
                                aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTAATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTTTTAtttttttctt                                                                                                                   >  2:483299/1‑143 (MQ=255)
                                  ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                 >  1:532660/1‑149 (MQ=255)
                                    aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGGGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                               >  2:505612/1‑149 (MQ=255)
                                    aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                               >  2:381056/1‑149 (MQ=255)
                                          aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                         >  1:326404/1‑149 (MQ=255)
                                              ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCTGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                     <  1:202094/149‑1 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGCGTAATATGTTTTCCTTTATTTTCCCAGAAGTCCGTATTTTGTGGTTTCCAGTCTCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                          <  1:104403/149‑1 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTGGCGAA‑A‑AA‑AGCTAAAAACTATATTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                          >  2:275970/1‑149 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGGGCAATATGTTTTCCTTAATTTTCCCAGAAGGGCGGATTTTCTGGTTTCCAGTCGCG‑A‑A‑AACAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                       >  1:436717/1‑149 (MQ=255)
                                                             actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                      <  2:326404/149‑1 (MQ=255)
                                                                gatttatatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGAGAA‑A‑AA‑AGATACCGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                   <  2:90518/144‑1 (MQ=255)
                                                                   tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGAGACTAATTCATACTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                                >  1:518955/1‑149 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGACGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                          <  2:436717/149‑1 (MQ=255)
                                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAA‑AA‑AGATAATGACTATTTTTCTCATTTTGCAGTGTCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                        <  1:275970/149‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGCTTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCAACAAATGTTAATGAAAAAAAATCaaacaccg                                                               >  1:261072/1‑144 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGCTTTTCTGGTTTCCAGTAGAGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                               >  1:77819/1‑149 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCAAGTCGCG‑AAA‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                               >  1:49599/1‑149 (MQ=255)
                                                                                         gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                                          >  2:171900/1‑149 (MQ=255)
                                                                                          cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                         <  1:219997/149‑1 (MQ=255)
                                                                                            atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAgg                                                       >  1:184107/1‑149 (MQ=255)
                                                                                                  tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGAAAACGACAATTTTTCTCAATTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                                 >  1:62519/1‑149 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                             >  2:154486/1‑149 (MQ=255)
                                                                                                             ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                      <  1:303668/149‑1 (MQ=255)
                                                                                                                              tgtCTGGTTTCCATTCGAGAA‑A‑AA‑AAATAAAGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACTAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt                     <  2:354101/147‑1 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                  <  2:532660/149‑1 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCA‑A‑ATAA‑AGATAACGACTATTTTTTGCATTTTGAAGTGCAATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                  <  1:252366/149‑1 (MQ=255)
                                                                                                                                    gTTTCCAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGAGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGa               >  1:368024/1‑149 (MQ=255)
                                                                                                                                        ccAGTCGCGAA‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACATGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTg           <  1:171900/149‑1 (MQ=255)
                                                                                                                                             tcgcgaaa‑aa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa       >  1:502938/8‑149 (MQ=255)
                                                                                                                                                  aaa‑aa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACTAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa  <  1:72562/147‑1 (MQ=255)
                                                                                                                                                  |||                                                                                                                                                    
TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑A‑AA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAA  >  NC_000913_3_vch_tpiA/1543877‑1544169

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