Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,787 G→A 10.8% intergenic (+261/+279) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,7870GA10.8% 106.3 / 2.7 39intergenic (+261/+279)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base G (8/26);  new base A (2/2);  total (10/29)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGC  >  NC_000913/1286638‑1286927
                                                                                                                                                     |                                                                                                                                            
cggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                                                                                                                                              <  2:1124197/149‑1 (MQ=255)
   cccGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                                           >  1:690377/1‑150 (MQ=255)
     cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                                         >  1:302712/1‑150 (MQ=255)
      ggTAAACGTGGAGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCa                                                                                                                                        <  1:408097/150‑1 (MQ=255)
          aaCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTc                                                                                                                                    <  1:250438/150‑1 (MQ=255)
           aCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCg                                                                                                                                   <  2:846379/150‑1 (MQ=255)
           aCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCg                                                                                                                                   >  2:1349510/1‑150 (MQ=255)
           aCGTGGTGCCATAGGTCAAAATTTGCTCGTAGTAACAGTTTTTGTTTTTTTAAGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                   <  2:893507/150‑1 (MQ=17)
            cGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGa                                                                                                                                  <  2:806879/150‑1 (MQ=255)
              gggagCCCTCGGGCAAAATTCGCGCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAc                                                                                                                                <  1:136487/146‑1 (MQ=18)
               ggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAcc                                                                                                                               <  1:557913/150‑1 (MQ=255)
                gtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGGGATGGTGGTTGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCt                                                                                                                              <  1:611657/150‑1 (MQ=255)
                 tgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCtt                                                                                                                             <  2:1083865/150‑1 (MQ=255)
                    cATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGt                                                                                                                          <  1:989283/150‑1 (MQ=255)
                       cGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGa                                                                                                                       <  2:392232/150‑1 (MQ=255)
                       cGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGa                                                                                                                       <  1:353262/150‑1 (MQ=255)
                       cGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                       >  1:786212/1‑150 (MQ=255)
                                ttAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                              <  2:690377/150‑1 (MQ=255)
                                   gcgcGGAGTAACCGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                           <  1:537716/147‑1 (MQ=255)
                                   gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                           >  2:981835/1‑150 (MQ=255)
                                    cTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTGGAAGCTTCTCATCCt                                                                                                          <  2:1340172/150‑1 (MQ=255)
                                      cGGAGTAACAGTTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                                        <  2:575612/150‑1 (MQ=255)
                                        gAGTAACTGGTTTTGATTATTTAAGGAATCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                      <  2:728953/150‑1 (MQ=255)
                                        gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTccc                                                                                                      <  2:786212/150‑1 (MQ=255)
                                              cAGGTTTTGATTATTTAAGGAAGCGCTGGTGGTGGTGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGAAACGCTCTCTCGCTGGCGCTCAAGTCGAACATTGGTCGAAGCTTCTCATCCTTCCCCGCTTg                                                                                                <  1:692313/150‑1 (MQ=255)
                                              cAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                <  2:567681/150‑1 (MQ=255)
                                                ggTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                              >  1:926642/1‑150 (MQ=255)
                                                 gTTTTGAGTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTGGTCGAAGCTGCTCATCCTTCCCCGCATGGGc                                                                                             >  2:167106/1‑150 (MQ=255)
                                                           tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCACGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                   <  2:1164178/150‑1 (MQ=255)
                                                             tAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGAt                                                                                 <  1:753883/150‑1 (MQ=255)
                                                                         ggtggtggGGGAAGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTGTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCGCCGCTTGGGCTGAATATTTGATTGCGGATTCGtt                                                                     <  2:1349349/150‑1 (MQ=11)
                                                                         gggggtggGGGATGGATTACTCAGCGGTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGGCGCTAAAGTCGAAGCTTGGTCGAAGCATCTCATCCTTCCCCGCTTGGGCAGAATATATGATTGCGGATTCGtt                                                                     <  2:1344203/147‑1 (MQ=11)
                                                                                        aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGc                                                      <  2:850344/150‑1 (MQ=255)
                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTg                                                 <  1:774578/150‑1 (MQ=255)
                                                                                             tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTg                                                 <  1:1349510/150‑1 (MQ=255)
                                                                                                               cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                               <  2:470727/150‑1 (MQ=255)
                                                                                                               cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTgggg                               >  1:722201/1‑150 (MQ=255)
                                                                                                                   tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAg                           >  1:303811/1‑150 (MQ=255)
                                                                                                                                       aCGCTCTCTCGCTGCCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGt                                            <  2:28190/113‑1 (MQ=14)
                                                                                                                                       aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGt                                            >  1:28190/1‑113 (MQ=11)
                                                                                                                                            ctctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCGTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc  >  1:1167813/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGC  >  NC_000913/1286638‑1286927

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