Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,543,636 C→T G142D (GGT→GAT)  narU ← nitrate/nitrite transporter NarU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,6360CT97.0% 107.3 / NA 33G142D (GGT→GAT) narUnitrate/nitrite transporter NarU
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/17);  minor base G (0/1);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATT  >  NC_000913/1543405‑1543875
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
ggTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCAGTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCaa                                                                                                                                                                                                                              <  1:24929/251‑1 (MQ=255)
    cccAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGcg                                                                                                                                                                                                                          <  1:342255/251‑1 (MQ=255)
         aTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGcgataac                                                                                                                                                                                                                     <  1:376767/251‑1 (MQ=255)
                cGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAat                                                                                                                                                                                                              >  2:252699/1‑251 (MQ=255)
                cGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATGGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCAGCTGCAAAACCGCATAGCAAAGCGATAACGATAAat                                                                                                                                                                                                              <  2:247360/251‑1 (MQ=255)
                 ggCATTCGCCAGCGACATCACAGAACCGACGGCCTGCGGGACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGCTAAat                                                                                                                                                                                                              >  2:378223/1‑250 (MQ=255)
                        gCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCaaa                                                                                                                                                                                                      <  1:270342/251‑1 (MQ=255)
                             cGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGCCCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAg                                                                                                                                                                                                    <  1:228183/248‑1 (MQ=255)
                             cGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAg                                                                                                                                                                                                    >  2:228183/1‑248 (MQ=255)
                                   cACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc                                                                                                                                                                                            <  1:330414/250‑1 (MQ=255)
                                               ggCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGCACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACgg                                                                                                                                                                               >  2:32898/1‑251 (MQ=255)
                                                                     gaactaaaGGCAAATACAGGTAAAAAAAGGACCAGCGGTGCAACCAGCTGCATTTCACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTGGCACGGCAATTCCGAGCCAGACGCAAgg                                                                                                                                                         <  2:402277/246‑1 (MQ=255)
                                                                          aaaGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAgg                                                                                                                                                         >  1:402277/1‑246 (MQ=255)
                                                                                      aggaagaaaactgaaaaGCGGTGCAACCCGCTGCATTTCACTTACACATAAGTTTACTAATCCGCGATTAATCCCAACAACACTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCACGCTCGAAGCATAGTTTGCATCAGCAAAACCGCATCGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGACAGAAGCAAGGAATAACCAGGATTGc                                                                                                                                          >  2:192045/17‑249 (MQ=255)
                                                                                      aGGTACAAAAATGCCCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGaa                                                                                                                                                       <  2:49206/236‑1 (MQ=255)
                                                                                      aGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGaa                                                                                                                                                       >  1:49206/1‑236 (MQ=255)
                                                                                          aCAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACt                                                                                                                                    <  1:152869/251‑1 (MQ=255)
                                                                                            aaanaTGACCAGCGGTTCAACAAGCTACATTACACTGACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGATCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAATCAAAGTTTGCATCTGCAAAACCGCATATCAAAGCGATAACTATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACTCATGGAATAATCCGGATTGCAGTACTaa                                                                                                                                  <  2:368996/251‑1 (MQ=255)
                                                                                               aaTGACCAGCGGAGCAACCAGCTGAATTACACTTACACCTAAGTTTCCTAATCCGCCCTTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATCTCCCACAAGGAGTATTCGGATTTTGCACGGCAATTCAGAGCCAGACGCAAGGAATAATCAGGATTGCAGTAATaaaa                                                                                                                                >  2:34553/1‑250 (MQ=255)
                                                                                                                         ttacacttacacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAatat                                                                                                      >  1:1163/1‑250 (MQ=255)
                                                                                                                           ccacttacacCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAg                                                                                                    <  2:104382/249‑1 (MQ=255)
                                                                                                                                     cTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATaa                                                                                           >  1:141790/1‑249 (MQ=255)
                                                                                                                                                        aTTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCATGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCaata                                                                      <  1:374893/251‑1 (MQ=255)
                                                                                                                                                                 aaGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGcc                                                               >  1:114136/1‑249 (MQ=255)
                                                                                                                                                                           ccTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAggg                                                   <  1:4100/251‑1 (MQ=255)
                                                                                                                                                                           ccTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAggg                                                   <  1:252699/251‑1 (MQ=255)
                                                                                                                                                                                                 agaaACTGATATTGCCCATGCTCGAAGCAAAGTTTGGATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGt                             >  1:126972/1‑251 (MQ=255)
                                                                                                                                                                                                        gATATTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCgg                      <  1:121552/251‑1 (MQ=255)
                                                                                                                                                                                                         atatTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGgt                     >  2:35315/1‑251 (MQ=255)
                                                                                                                                                                                                           atTGCCCATGCTCGAAGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGgta                    <  1:341126/250‑1 (MQ=255)
                                                                                                                                                                                                                          aGCAAAGTTTGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGCAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCg     >  2:223735/1‑250 (MQ=255)
                                                                                                                                                                                                                            cAAAGTTTGCATCTGCAACACCGCACAGCAAGGCGATAACGACAGATATACCAAACGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGCATAATCAGGATTGCAGTACTACAAACCGTCCATCGGCGTGCACCGAATATAGGCACCATACAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGGAATGCGGTTAATACAAAGAGTTGATCGGTCGTAAAATTAAAACCGAtt  >  2:165835/1‑251 (MQ=255)
                                                                                                                                                                                                                            cAAAGTTGGCATCTGCAAAACCGCATAGCAAAGCGATAACGATAAATAGCCCAAAAGGAGTAGTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAAc                                                         <  2:382888/196‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
GGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTTTGCACGGCAATTCCGAGCCAGACGCAAGGAATAATCAGGATTGCAGTACTAAAAACCGTCCATCGGCGTCCACCGAATATAGGCACCATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATT  >  NC_000913/1543405‑1543875

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: