Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,536,304 C→T G105D (GGC→GAC)  narV ← nitrate reductase Z subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,536,3040CT91.4% 96.9 / NA 35G105D (GGC→GAC) narVnitrate reductase Z subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/17);  minor base G (0/3);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

AGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGT  >  NC_000913/1536058‑1536532
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    
aGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAAt                                                                                                                                                                                                                                  >  1:208386/1‑251 (MQ=255)
                            ccctctcTGACGATCCACCACGTAAATCTACAATTCTTTGCGCCCAACCCACCAGCTTCATCATTCCGCTACCATCAGTATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGTCACTGGATCAGCAAAATGCTCATGATGAGGATACCCGGCGTAGTGGAAGGGGCACGCACGCGCTGATTTATCAGCCTGCGCCACAGCATTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATAccgccg                                                                                                                                                                                                      <  2:169098/248‑1 (MQ=255)
                                 tctGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCAc                                                                                                                                                                                                 >  1:19168/1‑251 (MQ=255)
                                  ctGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCAcc                                                                                                                                                                                                >  1:449454/1‑251 (MQ=255)
                                             ccacGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAg                                                                                                                                                                                     <  1:335420/251‑1 (MQ=255)
                                              cacgataaaggacatgctcttGCGCCCACCACACCAGATCCATCCTTACGCTACCATAAGGATACTGTGCTGAAAACGGTATGGTACCGAGCCCCAACAGACACGGGAGCAGCAAATGGCTCATGATGATGATACCCGGCGTAGTGGAAGTGGCACGCACACGCTGATTTGTCAGCCTGCGCCACAGCATACCTGCGCCGGCAATAAACGTCAACACGCCGCAATTACCGCCGAGCACCATCGACATGAgt                                                                                                                                                                                    <  2:436519/232‑1 (MQ=255)
                                               acGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAgtt                                                                                                                                                                                   <  1:169601/251‑1 (MQ=255)
                                                 ggAAAGTGACAATGCTTTGCGCCCAACCCCCCAGCTTCATCATTGCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGCCACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGTTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACTCCGCAAATACCGCCGAGCACCCTCGCCATGAgttgt                                                                                                                                                                                 <  2:412739/251‑1 (MQ=255)
                                                              gCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACt                                                                                                                                                                     <  1:426284/250‑1 (MQ=255)
                                                                     gcCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGCGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTAAATCAACGTCAACCCACCGCTAATACCGCCGAGCACCACAGCCATGAGTTGTTTCGCTGCCACGGGCaaaa                                                                                                                                                              >  2:214693/1‑250 (MQ=255)
                                                                     ccccAACCCACCAGCTTCATCATTTCGCTACCATCAGGATTCTGTGCTGAAAACGTTTTGGTACTGAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCCCGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACCCCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCaaaaa                                                                                                                                                             <  2:449454/250‑1 (MQ=255)
                                                                         aaccaacctacTTCATCATTTCCCTACCATTAGGTTACTGATCTGAAAGCGGTATGGTACTTGGCCCCAACCGACACTGGATCAGCAAAATGCTCATGATGATGATAACCGTCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCTTCAATCAACGTCAACGCCCCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGATGCCACTGGCAAAAACCAc                                                                                                                                                         <  2:417945/241‑1 (MQ=255)
                                                                                 agcatcagcattccgctcccaataGGATAATGTGACGAAAACGGAATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGGCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGGTGTTTCGCTGCCACTGGCAAAAACCACGCGTACAt                                                                                                                                                 <  2:78556/228‑1 (MQ=255)
                                                                                     tcatcaTTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATcc                                                                                                                                               <  1:88761/249‑1 (MQ=255)
                                                                                     tcatcaTTGCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACTCCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCa                                                                                                                                              <  1:42272/250‑1 (MQ=255)
                                                                                                cTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTa                                                                                                                                   >  2:412648/1‑250 (MQ=255)
                                                                                                  aCCAGCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAAc                                                                                                                                 <  1:463614/250‑1 (MQ=255)
                                                                                                         ggataatgagcagagaaaggTATGGTACTTAGCCCCAAAAGAAACTGGCTCAGCGAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTACTGCGCCGTCAATCAACGTCCACACACCGCAAATACCGCCGAGCGCCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACAACGCGAACATCCAGGGCGGCGATAACATGCCGa                                                                                                                          >  2:81974/19‑250 (MQ=255)
                                                                                                                              tagggtacatgCCCCAAACAGCCAATGGATCAGCAAAATGCCCATGATGATGATATCCGGCGTAGTGGAAGGGGCACGCACGCGCTGATTTGTCACCCTGCTCCACAGCAGTCCGGCGCCGGCAATCAACGTCAACCCACCGCAAATACCGCCTAGCACCATCGCCATGGGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCTGAACAAGTGCCCGCAGAAAATc                                                                                                      <  2:152312/240‑1 (MQ=255)
                                                                                                                              aggttcgttcgCCCCAACAGACCCTGGTACTGCAAAAGGCTTAGGAGTATGATATCCGGCGTAGTGTACGTGGCACGCAGTCGCTGAGTTGCCAGCCAGCGCCCCCGCAGTCCGGCGCCGTCAATCAACGTCAACACACCGCAAATCCCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCCCGCGTACATCCAGTGCGGCGTTAACATGCCGATCATGTGCCCGAAGAATATccc                                                                                                    <  2:163665/241‑1 (MQ=255)
                                                                                                                               tGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCCTGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCACCACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTAGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACACGTGCCCGAAGAAAATcccc                                                                                                   >  2:199097/1‑251 (MQ=255)
                                                                                                                                    cgtAGCCCCAACAGAACCTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGGGGAAGTGGCACGCACGCGCTGACTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCTCAAATAACGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAATCCACGCGTACATCCAGTGCGTCGTTAACATGCCTAACAAGTGCCCGAAGAAAAGCCCCAAAAt                                                                                              <  2:105548/249‑1 (MQ=255)
                                                                                                                                                        ggATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTc                                                                           <  1:89076/251‑1 (MQ=255)
                                                                                                                                                             aGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGAcc                                                                       >  1:2260/1‑250 (MQ=255)
                                                                                                                                                                                tgatATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGttt                                                    >  1:501681/1‑250 (MQ=255)
                                                                                                                                                                                 gatATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATAcc                                                                                                                                                                                                          >  1:323999/1‑99 (MQ=255)
                                                                                                                                                                                 gatATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGGCAATCAACGTCAACACACCGCAAATAcc                                                                                                                                                                                                          <  2:323999/99‑1 (MQ=255)
                                                                                                                                                                                       cGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTc                                                                           >  1:86759/1‑220 (MQ=255)
                                                                                                                                                                                                       gcacgcacGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGgc                            <  1:199097/251‑1 (MQ=255)
                                                                                                                                                                                                         acgcacgcGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGgcgc                          >  1:169282/1‑251 (MQ=255)
                                                                                                                                                                                                               gcgcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTTAGGAGGCGCGCCAgg                    >  1:240637/1‑251 (MQ=255)
                                                                                                                                                                                                                cgcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGgt                   <  1:538705/251‑1 (MQ=255)
                                                                                                                                                                                                                    gATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCTCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTa                <  1:412648/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcCACACCAGTCCTGCGACGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGAACATATCACCACCCCGCGTTCATCGAGCATTTGACTTGAGGAGGCGCGCCAGGGGTATTGACCGTAGTCGt  >  2:239532/1‑249 (MQ=255)
                                                                                                                                                                                                                                             tCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCAttt                                         >  1:533739/1‑200 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    
AGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTCGT  >  NC_000913/1536058‑1536532

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