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→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/19);  total (17/19)

GAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑‑C‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_vch_tpiA/1543881‑1544170
                                                                                                                                                |||                                                                                                                                                     
gagaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑‑A‑a‑aaaaga                                                                                                                                                >  1:771548/1‑146 (MQ=255)
                    tAAATAGCATCCAGAAACAGAAGGCAAGAAGTAGACAACTGATTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAATCGC‑AA‑A‑A‑AAA‑GATAACGACTATTTTTCTCAtt                                                                                                                            <  2:492801/149‑1 (MQ=255)
                        tAGCATCCAGACAGAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑‑A‑A‑AAACGATAACGACTATTTTTTTAAtttttc                                                                                                                        >  2:485778/1‑147 (MQ=255)
                        tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGc                                                                                                                        >  2:668060/1‑149 (MQ=255)
                         aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:481702/149‑1 (MQ=255)
                            aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATCTCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAAA‑A‑AAA‑CATAAGGACAATTTTTCTCATTTTGCAgg                                                                                                                     >  1:255447/1‑148 (MQ=255)
                              ccAGACACAGGAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTACGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  1:224174/149‑1 (MQ=255)
                                    acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGGGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                            >  1:333702/1‑149 (MQ=255)
                                       gAAGGCATGAAGTAGACAATTGACTGATATCCAAAGATTGCTTCGAGCAATATGTTTTCCTTTATTTTACCAGATGGCCGGATTTTTTGGTTAAAAAAAAA‑AA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  2:576923/149‑1 (MQ=255)
                                              agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATGGTTCCTCACAt                                                                                                  >  1:768501/1‑149 (MQ=255)
                                                  gTACACAACTGACTGATAACCAGAGATTTATTCGAGCAATATGTTTTCCTTTATTGACCCAGAAGGCCGGATTTTCTGGTTTAAAAAAAA‑AA‑A‑A‑AAA‑GATAACGACTATTATTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                              <  2:739685/149‑1 (MQ=255)
                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCTC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          <  1:732953/149‑1 (MQ=255)
                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCTA‑‑A‑A‑AAA‑GATAAGGATTATTTTTCTCATTTTATAGTGCCATATTATTCTTCACATGCAAACAt                                                                                          >  2:149966/1‑149 (MQ=255)
                                                             actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACACATTgtaggt                                                                                   >  2:594644/1‑144 (MQ=255)
                                                               tgaTATCCAGAGATTTCTTCGAGCAATATGTGTTCCTTTATTTTCCCAAAAGGCCGCAGTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTCACCGTGCCGTATCGTCAGTCACATGCACACATGTGTAAtaa                                                                                 >  2:586303/1‑147 (MQ=255)
                                                                    tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                            >  2:623464/1‑149 (MQ=255)
                                                                      cAGAGATATCTTCGATCAATAAGTTTTCCTTTATTATCCCACAAGGCCGGATTTTCTTGTTTACAGACGA‑A‑‑AAA‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          <  2:596147/149‑1 (MQ=255)
                                                                            tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    >  1:488301/1‑149 (MQ=255)
                                                                               cttcGAGCATTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                 <  1:217907/149‑1 (MQ=255)
                                                                                  cGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg                                                              <  1:668060/149‑1 (MQ=255)
                                                                                   gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  1:463328/149‑1 (MQ=255)
                                                                                   gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTACTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  1:62639/149‑1 (MQ=255)
                                                                                       aaTATGTTTTCCTAAAAATAACAAGAAGGCCGAATTTTCTGGTTTCGAATCGA‑G‑‑A‑ACAAA‑GATAACGAATATATTAAACATCTTCAAGAATTATATTGTTCCTCTCATGCATAGATAAAAAAAAAAAAAAAGACAAAACACGAGGTAAg                                                         <  2:798592/149‑1 (MQ=255)
                                                                                           tGTTCTCCTTTTTTTTCCCAGTAGGCTGGATTTTCTTGTTTCCAGTGGC‑GA‑A‑A‑AAA‑GATGACGCCTCTTTTGTTCATTTTGCAGTGCCATCTTGTTCCTCTCATTAACACATTGGTAATGAAAAAAAGACAGAACACGAGGTAAGgcgc                                                     <  2:609828/149‑1 (MQ=255)
                                                                                               ttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAAGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 >  1:787929/1‑149 (MQ=255)
                                                                                               ttCCTTTATTTTCCAAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACTACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 <  2:538651/149‑1 (MQ=255)
                                                                                                 ccgttATTTTCCCAGAATCTCGGATTTTCTTGTTTCCAGGCGC‑GA‑A‑A‑AAA‑GATAATGACTATTTTTCTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAATGACAAAACACGAGGTAAGGCGCAATAGc                                               <  2:815067/146‑1 (MQ=255)
                                                                                                  cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                              >  1:437616/1‑149 (MQ=255)
                                                                                                    ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  1:391882/149‑1 (MQ=255)
                                                                                                    ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  1:268863/149‑1 (MQ=255)
                                                                                                                 aGGCCGGATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:571048/149‑1 (MQ=255)
                                                                                                                      ggATTTTCTGGTTTCCAGTCGC‑GA‑A‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg                          >  1:757763/1‑149 (MQ=255)
                                                                                                                          tttCTGGTTTCCAGTCGC‑G‑‑AAA‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt                      >  2:25834/1‑149 (MQ=255)
                                                                                                                            tCTGGTTTCCGGTCGA‑GA‑A‑A‑AAA‑TATAACGACTAGTGTCACAATTTTTCGGTGCCATATTTTTCCTAGCATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGTAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  2:392358/149‑1 (MQ=255)
                                                                                                                                      aGTCGC‑G‑‑AAA‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAGTTGGGTGAGGTGCt          >  2:484057/1‑149 (MQ=255)
                                                                                                                                      aGTCGC‑G‑‑AAA‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  1:420052/1‑149 (MQ=255)
                                                                                                                                           tcgcgaaa‑aaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        >  1:572324/8‑149 (MQ=255)
                                                                                                                                                 aa‑aaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:189859/148‑1 (MQ=255)
                                                                                                                                                 aa‑aaa‑GATAACGACGATTTTTGTCATTTTGCAGTGCCATATTGTTCGTCAGAGGCACACAGTGGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:477008/148‑1 (MQ=255)
                                                                                                                                                |||                                                                                                                                                     
GAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑‑C‑A‑AAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  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: ‑