Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,787 G→A 14.3% 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,7870GA14.3% 90.4 / 4.1 37intergenic (+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/23);  new base A (3/2);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCC  >  NC_000913/1286638‑1286929
                                                                                                                                                     |                                                                                                                                              
aGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                                                                                                                                                >  1:356328/1‑150 (MQ=255)
 ggCCCTGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                               <  2:810447/150‑1 (MQ=255)
 ggCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                               <  2:1254802/150‑1 (MQ=255)
   cccGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                                             >  1:481191/1‑150 (MQ=255)
    ccGTTAAACGTGGTGCCTTCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                            <  2:395830/150‑1 (MQ=255)
    ccGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTTATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                                            <  2:788349/150‑1 (MQ=255)
     cGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                                           >  1:876105/1‑150 (MQ=255)
        tAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAg                                                                                                                                        >  2:30200/1‑150 (MQ=255)
         aaaCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGt                                                                                                                                       <  1:552492/150‑1 (MQ=255)
         aaaCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATAATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGt                                                                                                                                       <  1:669062/150‑1 (MQ=255)
                    cATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCCAGTCGAACCTTAGt                                                                                                                            >  1:222483/1‑150 (MQ=255)
                      tatggCAAAATTTGCTCGGAGTATCAGTTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCg                                                                                                                          <  2:222483/147‑1 (MQ=17)
                         gccAAAATTAGCTCGGAGTAACAGCATGTCATTCTTTAACGAAGCGATGCGGTGGGGGGATCGACTACTCCGCGCTGCGCGCGTCGCCCTTCGGGTCTTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCAAGTCGAACCTTGGTCGAAg                                                                                                                       <  1:191322/148‑1 (MQ=17)
                          gCAAACTTAGCTCGGAGTACCAGGTGTTGATTATTTAAGGAAGCGATGGGTGTGCGGGAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGGCGTGGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                                      <  2:984452/150‑1 (MQ=18)
                             aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTTGTCGAAGCTtc                                                                                                                   >  2:822195/1‑150 (MQ=255)
                             aaaTTAGCTCGGAGTAACAGGCTCTGACTATTTAAGGAAGCGATCGCGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtc                                                                                                                   <  2:331636/150‑1 (MQ=17)
                               aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAAGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                 <  1:93454/150‑1 (MQ=255)
                                   gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                             >  2:158263/1‑150 (MQ=255)
                                   gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                             <  2:1222678/150‑1 (MQ=255)
                                        gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                        <  2:1284903/150‑1 (MQ=255)
                                               aGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAgc                                                                                                                                       <  1:218140/112‑2 (MQ=255)
                                               aGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAgc                                                                                                                                       >  2:218140/1‑111 (MQ=255)
                                                ggTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTggg                                                                                                <  2:142077/150‑1 (MQ=255)
                                                      gATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                          <  1:158263/150‑1 (MQ=255)
                                                         tattTAAGGAAGCGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAtt                                                                                       <  2:244195/150‑1 (MQ=255)
                                                                ggAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGGCGAACCTTTGTCGAAGCTTCGCATCCATCCCCGCATGGGCAGAATATTTGAtttc                                                                                >  2:395690/1‑148 (MQ=11)
                                                                       aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCg                                                                         <  1:882470/150‑1 (MQ=37)
                                                                          gtggtggGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGttt                                                                      <  2:813489/150‑1 (MQ=255)
                                                                            ggggggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTTCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTga                                                                    <  1:1072521/147‑1 (MQ=255)
                                                                             gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgag                                                                   <  1:799816/150‑1 (MQ=40)
                                                                             gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgag                                                                   <  1:218981/150‑1 (MQ=40)
                                                                                      ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGACCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                          >  2:533502/1‑150 (MQ=255)
                                                                                      ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                          <  2:1267029/150‑1 (MQ=255)
                                                                                      ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                          <  1:97852/150‑1 (MQ=255)
                                                                                       gagtACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCTCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCgggg                                                         <  1:1193816/147‑1 (MQ=255)
                                                                                         ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCt                                                       >  1:837309/1‑150 (MQ=255)
                                                                                            cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCtttt                                                    <  1:864501/150‑1 (MQ=255)
                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtg                                          <  2:286450/150‑1 (MQ=255)
                                                                                                                      ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTGGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                          >  1:515265/1‑150 (MQ=255)
                                                                                                                       gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                         >  2:1115353/1‑150 (MQ=255)
                                                                                                                                              ctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc  <  1:728229/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AGGCCCGGTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCC  >  NC_000913/1286638‑1286929

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: