Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 16,749 T→A intergenic (+192/+2) insL1 → / ← hokC IS186/IS421 family transposase/protein HokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091316,7490TA89.7% 75.1 / 2.9 29intergenic (+192/+2)insL1/hokCIS186/IS421 family transposase/protein HokC
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (21/5);  total (22/7)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTT  >  NC_000913/16512‑16977
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
gCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTc                                                                                                                                                                                                                          >  1:52727/1‑251 (MQ=255)
      tttCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGc                                                                                                                                                                                                                    >  1:204937/1‑251 (MQ=255)
                         ggATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTcc                                                                                                                                                                                                 >  1:236544/1‑251 (MQ=255)
                              cGAAACGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGAACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTg                                                                                                                                                                                            >  1:303611/1‑251 (MQ=255)
                              cGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTg                                                                                                                                                                                            >  1:129456/1‑251 (MQ=255)
                                                 gTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAAcc                                                                                                                                                                         >  1:413155/1‑251 (MQ=255)
                                                      ggAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGCCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATGGAAGAACGCGAGGCATCGTCTTAACGAGGCCCCGAGGCGTCGCATTCTTCAGATGGTTCACCCCTTAATTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTcaca                                                                                                                                                                    <  1:26517/251‑1 (MQ=255)
                                                          aatcccAAAAATGGGCACCTGGAGCTTACAGCAGGCCATACGAGGGACAGCATCCCTGACCGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGCGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTCAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACACGGATTCGCACGCCGTGAAAACAGCAACCTCCGTCTGGCCAGCTCGGATGTGAACCTCAAAGAgg                                                                                                                                                                >  2:223063/6‑251 (MQ=18)
                                                                      ggTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCCACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGtt                                                                                                                                                    >  1:498713/1‑251 (MQ=255)
                                                                          aTCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTAcc                                                                                                                                                 >  1:465089/1‑250 (MQ=255)
                                                                            gtgggccttccgggTGGCCATTCGGGGGACAGTACCCCTGACAGCCTAAAAAACGCAATTGAAGAACGCGAGGCATCGTCTAACCGAGGCCCGGGGGCGCCGCTTTCTTCAGATGGTTCAACCCTAAAGTTAGCGCTTATGGGATCACCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAgc                                                                                                                                              <  1:289355/239‑1 (MQ=25)
                                                                                  cTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGGCTTTTCTCGTTACCAGCgccgcc                                                                                                                                        >  1:401790/1‑251 (MQ=255)
                                                                                                  cGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATAc                                                                                                                        >  2:137474/1‑251 (MQ=255)
                                                                                                                 cTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGCTCAggg                                                                                                         >  1:34132/1‑251 (MQ=255)
                                                                                                                            cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGCAACAGATGACGGTCCGGGGGCCAatcatc                                                                                              >  1:247329/1‑251 (MQ=37)
                                                                                                                                                 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGAAGACGATCAGGGCGACAAAACTCGCCTTAtggagcttttttgct                                                                         >  2:99659/1‑238 (MQ=38)
                                                                                                                                                 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGAAACAGATGACGAACAGGGCGGCAATAATAGACATATGCTGCTTCATTGct                                                                         >  1:62093/1‑251 (MQ=38)
                                                                                                                                                                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCCTCTTCATTGCTCTCTTCTCCTTGAACTTACGGTc                                                 >  1:329830/1‑251 (MQ=255)
                                                                                                                                                                         gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTtc                                                                                                                                                         >  2:180734/1‑147 (MQ=255)
                                                                                                                                                                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAgc                                                                                                                                              >  2:363271/1‑148 (MQ=255)
                                                                                                                                                                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCCGGGCGACAATCATCGCCCTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACAtt                             >  1:99928/1‑250 (MQ=255)
                                                                                                                                                                                                          atgtagtcgtttatggaactccTCCCCGCCGTTGCACTTACTCGGATGCGTAAGCCGTGAAAACACCAACCGACGTCTGGCCATTTCGGATGTGAACCCCACAGAGTTCTTTTCTCTTTACCAGGGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATag               <  1:339829/230‑1 (MQ=255)
                                                                                                                                                                                                              tAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGAGATACAGATGACGATCACGGCGAAAATAATCACCTTTTGCGGCTTCATTGATCTCTTCTCCTTGGACTTTCGGTCCGTAAGAGGCCCTCTAACTGTGtttcgcattatggggg           >  1:291919/1‑237 (MQ=255)
                                                                                                                                                                                                                gcgttttcggcctccccccccGCCGTTGCACTTCCTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGGCTGGCCAGTTCGGATGGGAACCTCACAGAGGTCTTTTCTCGTAACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGcc         <  1:302872/234‑1 (MQ=255)
                                                                                                                                                                                                                gcgcTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGGTCAAGATTGCGATCAGAGAGGCAAACACCCCCCTATGCTTCGTCATTGCTCTCTTCTCCTTTGACTTCCGGTCAGTAAGAGGCACTCTACATGTGTTTAGCGTATAGGGGGcc         >  2:381748/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtttgttggcctccccccccGCCGTGGCATTTACTCGGATTCTTAGGCCGTTAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCt        <  1:462563/235‑1 (MQ=255)
                                                                                                                                                                                                                   cTTATGGGATTACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGTTAAAGATGGAGGTCAGGGGGGCCATCATCTCCTTTTGCTGCTTCATTTCTCTCTTTTCCTTGACCTTTCGGTCAGTTAAAGGCACTTTACATGTGTTCCACATATAAGGGGCCTCg      >  2:100885/1‑251 (MQ=255)
                                                                                                                                                                                                                      ggggcacccccccccGCCGTTCCACTTACTCGGATTCGTAAGCCGGAAAAACAGCAACCCCCGTCTGGCCGGTTCGGATGTGACCCTCACAGAGGTCTTTTCCCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGt   <  1:381748/240‑1 (MQ=255)
                                                                                                                                                                                                                       gggcccccccccccGCCGTTGCACTTACTCGGTTTCGTAAGCCGTAAAAACAGCCACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGtt  <  2:34132/241‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
GCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTT  >  NC_000913/16512‑16977

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