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.A77.8% 93.3 / 22.6 36coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (6/2);  new base A (11/17);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_vch_tpiA1,544,0220CA100.0% 115.4 / NA 36L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/19);  total (17/19)

GAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑‑CAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_vch_tpiA/1543881‑1544170
                                                                                                                                                |||                                                                                                                                                     
gagaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑‑Aaaaaaga                                                                                                                                                >  1:771548/1‑146 (MQ=255)
                    tAAATAGCATCCAGAAACAGAAGGCAAGAAGTAGACAACTGATTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAATCGCAAAAAAAGATAACGACTATTTTTCTCAtt                                                                                                                            <  2:492801/149‑1 (MQ=255)
                        tAGCATCCAGACAGAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑‑AAAAACGATAACGACTATTTTTTTAAtttttc                                                                                                                        >  2:485778/1‑147 (MQ=255)
                        tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGc                                                                                                                        >  2:668060/1‑149 (MQ=255)
                         aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:481702/149‑1 (MQ=255)
                            aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATCTCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAACATAAGGACAATTTTTCTCATTTTGCAgg                                                                                                                     >  1:255447/1‑148 (MQ=255)
                              ccAGACACAGGAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTACGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  1:224174/149‑1 (MQ=255)
                                    acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGGGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                            >  1:333702/1‑149 (MQ=255)
                                       gAAGGCATGAAGTAGACAATTGACTGATATCCAAAGATTGCTTCGAGCAATATGTTTTCCTTTATTTTACCAGATGGCCGGATTTTTTGGTTAAAAAAAAAAAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  2:576923/149‑1 (MQ=255)
                                              agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATGGTTCCTCACAt                                                                                                  >  1:768501/1‑149 (MQ=255)
                                                  gTACACAACTGACTGATAACCAGAGATTTATTCGAGCAATATGTTTTCCTTTATTGACCCAGAAGGCCGGATTTTCTGGTTTAAAAAAAAAAAAAAAGATAACGACTATTATTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                              <  2:739685/149‑1 (MQ=255)
                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCTCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          <  1:732953/149‑1 (MQ=255)
                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCTA‑‑AAAAAGATAAGGATTATTTTTCTCATTTTATAGTGCCATATTATTCTTCACATGCAAACAt                                                                                          >  2:149966/1‑149 (MQ=255)
                                                             actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACACATTgtaggt                                                                                   >  2:594644/1‑144 (MQ=255)
                                                               tgaTATCCAGAGATTTCTTCGAGCAATATGTGTTCCTTTATTTTCCCAAAAGGCCGCAGTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTCACCGTGCCGTATCGTCAGTCACATGCACACATGTGTAAtaa                                                                                 >  2:586303/1‑147 (MQ=255)
                                                                    tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                            >  2:623464/1‑149 (MQ=255)
                                                                      cAGAGATATCTTCGATCAATAAGTTTTCCTTTATTATCCCACAAGGCCGGATTTTCTTGTTTACAGACGAA‑‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          <  2:596147/149‑1 (MQ=255)
                                                                            tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    >  1:488301/1‑149 (MQ=255)
                                                                               cttcGAGCATTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                 <  1:217907/149‑1 (MQ=255)
                                                                                  cGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg                                                              <  1:668060/149‑1 (MQ=255)
                                                                                   gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  1:463328/149‑1 (MQ=255)
                                                                                   gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTACTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  1:62639/149‑1 (MQ=255)
                                                                                       aaTATGTTTTCCTAAAAATAACAAGAAGGCCGAATTTTCTGGTTTCGAATCGAG‑‑AACAAAGATAACGAATATATTAAACATCTTCAAGAATTATATTGTTCCTCTCATGCATAGATAAAAAAAAAAAAAAAGACAAAACACGAGGTAAg                                                         <  2:798592/149‑1 (MQ=255)
                                                                                           tGTTCTCCTTTTTTTTCCCAGTAGGCTGGATTTTCTTGTTTCCAGTGGCGAAAAAAGATGACGCCTCTTTTGTTCATTTTGCAGTGCCATCTTGTTCCTCTCATTAACACATTGGTAATGAAAAAAAGACAGAACACGAGGTAAGgcgc                                                     <  2:609828/149‑1 (MQ=255)
                                                                                               ttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAAGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 >  1:787929/1‑149 (MQ=255)
                                                                                               ttCCTTTATTTTCCAAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 <  2:538651/149‑1 (MQ=255)
                                                                                                 ccgttATTTTCCCAGAATCTCGGATTTTCTTGTTTCCAGGCGCGAAAAAAGATAATGACTATTTTTCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAATGACAAAACACGAGGTAAGGCGCAATAGc                                               <  2:815067/146‑1 (MQ=255)
                                                                                                  cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                              >  1:437616/1‑149 (MQ=255)
                                                                                                    ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  1:391882/149‑1 (MQ=255)
                                                                                                    ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  1:268863/149‑1 (MQ=255)
                                                                                                                 aGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:571048/149‑1 (MQ=255)
                                                                                                                      ggATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg                          >  1:757763/1‑149 (MQ=255)
                                                                                                                          tttCTGGTTTCCAGTCGCG‑‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt                      >  2:25834/1‑149 (MQ=255)
                                                                                                                            tCTGGTTTCCGGTCGAGAAAAAATATAACGACTAGTGTCACAATTTTTCGGTGCCATATTTTTCCTAGCATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGTAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  2:392358/149‑1 (MQ=255)
                                                                                                                                      aGTCGCG‑‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAGTTGGGTGAGGTGCt          >  2:484057/1‑149 (MQ=255)
                                                                                                                                      aGTCGCG‑‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  1:420052/1‑149 (MQ=255)
                                                                                                                                           tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        >  1:572324/8‑149 (MQ=255)
                                                                                                                                                 aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:189859/148‑1 (MQ=255)
                                                                                                                                                 aaaaaGATAACGACGATTTTTGTCATTTTGCAGTGCCATATTGTTCGTCAGAGGCACACAGTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:477008/148‑1 (MQ=255)
                                                                                                                                                |||                                                                                                                                                     
GAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑‑CAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_vch_tpiA/1543881‑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: