Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,249 C→A 17.7% intergenic (‑150/+83) thiC ← / ← rsd phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2490CA17.7% 78.8 / 4.0 44intergenic (‑150/+83)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base A (4/4);  total (21/23)
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.58e-01

AAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCA  >  NC_000913/4196110‑4196390
                                                                                                                                           |                                                                                                                                             
aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTGTCa                                                                                                                                            >  1:511918/1‑143 (MQ=255)
    aGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCAGCCa                                                                                                                                        <  1:241168/143‑1 (MQ=255)
       tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCAGCCAGtt                                                                                                                                     >  1:40747/1‑143 (MQ=255)
         aGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGCTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTaa                                                                                                                                   <  1:496498/143‑1 (MQ=255)
          ggtgtcggCTTGTCGGGCGACGGATGAATCAGGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACGCTGGCTCCCTTCGCAGGTATTAGCCTGATGAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAg                                                                                                                                  <  2:669702/137‑1 (MQ=255)
          ggACCTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCAGCCAGTTAAg                                                                                                                                  <  1:592314/143‑1 (MQ=255)
           gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGAACCCCCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCAGCCAGTTAAgg                                                                                                                                 >  1:25310/1‑143 (MQ=255)
            aCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGAGGCGATGCGATAGCAGATTCACTCTTGTTCCCGTCGCAGGTATTAGACTGATCAGGTTCGGAGGGTCACGAAATAACGGAATCAGCCAGTTAAGGc                                                                                                                                >  2:53276/1‑143 (MQ=255)
                ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCGCTc                                                                                                                            >  1:86714/1‑143 (MQ=255)
                   ttGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGCTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGa                                                                                                                         >  1:651728/1‑143 (MQ=255)
                         gacgacGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGCCAAGaa                                                                                                                   >  2:606269/1‑143 (MQ=255)
                             acgGAAGAATCAACTGACGAGTGCCGCAGGGGCGATGCGATAGCAGATGATCTCTTGCTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCGGCCAGTTAAGGCACTCCGACAAGAAATTa                                                                                                               <  2:590757/143‑1 (MQ=255)
                                gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCGTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcc                                                                                                            <  1:202917/143‑1 (MQ=255)
                                  tGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCCGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcccc                                                                                                          <  2:40859/143‑1 (MQ=255)
                                   gAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATCAGCCTGAGCAGCTTCCGCGGGTCCCGACTAAACGGTCTCAGCCAGTTAAGGCACGCCGACAAGAAAAGAGCcgcg                                                                                                         >  2:429577/1‑143 (MQ=255)
                                       cAAGCGACCATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa                                                                                                     <  1:375549/143‑1 (MQ=255)
                                        aaGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCCCAGCCAGTTAAGGCACTCCGACAAGAAAATAGCCCCGCAAAt                                                                                                    >  2:513829/1‑143 (MQ=255)
                                          gagcCGATCGGCGCCGGCGCGATGCGATAGCAGATGAGCTCTTGCGCTCTTCGCAGGTATTAGCCTGATCACGTTCCGCGCATCCTGTATAAACGGTCGCAGCCAGTTAAGGCACTCAGACAAGAAATTAGCCCCGCAAATgg                                                                                                  <  2:101414/143‑1 (MQ=255)
                                           agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg                                                                                                 >  1:157944/1‑143 (MQ=255)
                                           agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg                                                                                                 >  1:133153/1‑143 (MQ=255)
                                             aCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGCTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc                                                                                               >  1:288773/1‑143 (MQ=255)
                                                    cccctgcggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAAt                                                                                        >  1:206756/8‑143 (MQ=255)
                                                      gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTATCAGCCAGTTAAGGCACTCCGCCAAGAAATTAGCGCCGCAAACGGGGCGTTGAATGt                                                                                      >  2:213593/1‑143 (MQ=255)
                                                                     caGCATATTAACTCTTGGTCCCTTCGCAGGTATTAGCCTGATCAGGGTCCGCGGATCGCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAAc                                                                       <  2:321761/142‑1 (MQ=255)
                                                                           aTTAACTCCTGTTACCTTCGTAGGTATGAGCCTGATCAGGTTCCGCGGACCCCGAATAAACGGTCTCAGCCAGTTAAGGCACGCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                                                                 <  2:599829/143‑1 (MQ=255)
                                                                                     gctCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGATCTCAGCCAGTTAAGGCACTCCGCCAAGAAATTAGCCCCGCAAATGGGGCATTGACTGTAAATTACGGGTTAACAGCGCAGATCTCAAGc                                                       >  1:408296/3‑143 (MQ=255)
                                                                                           ttCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCTGAATAATCGGTCTCAGGCAGTTAAGGCACTCCGACAAGAAATGAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATg                                                 <  1:592121/143‑1 (MQ=255)
                                                                                              gCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACCAGAAATTAGCCCCGCAAATGGGGCATTGAAAGCAAATTACCCGTTAACAGCGCAGCACTCAAGCAGGATGtct                                              >  2:265588/1‑141 (MQ=255)
                                                                                                     tgaGCCTCATCCGGTGCCGCGGCTCCCGAATAAACGTTCTCAGCCAGTTACGGGACGCCGACAAGAAATTAGCCCCGCAAATGGAGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgg                                       <  2:399123/141‑1 (MQ=255)
                                                                                                      tcgCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcggg                                      <  1:588288/141‑1 (MQ=255)
                                                                                                        cccTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcg                                    <  1:505190/142‑1 (MQ=255)
                                                                                                         ccTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcgg                                   >  2:570717/1‑143 (MQ=255)
                                                                                                          gtggtcgggTTCGGCGGGTCGCGAATAAACGGGCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc                                  <  2:384336/136‑1 (MQ=255)
                                                                                                          cggaccAGGTTCCGCGGACCCCGAACAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc                                  <  2:206756/138‑1 (MQ=255)
                                                                                                                aGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAg                            <  1:369248/143‑1 (MQ=255)
                                                                                                                 ggTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGCTTGACGCGGGCGGCGTCAAGc                           <  2:288773/143‑1 (MQ=255)
                                                                                                                  gtcccGCGGATCCCGACTAAACGGTCTCAGCCAGTTAAGTCACTCCGACAAGAAATTAGCCCCGCAAATGGGTCATTGAATGTAAATTACGCGTGAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                          <  1:405343/140‑1 (MQ=255)
                                                                                                                  gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                          <  1:644/143‑1 (MQ=255)
                                                                                                                  gTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCa                          >  1:163939/1‑143 (MQ=255)
                                                                                                                   ttCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc                         >  1:617096/1‑143 (MQ=255)
                                                                                                                      cgcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGAGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGGACcag                      >  2:357855/1‑143 (MQ=255)
                                                                                                                        cgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagca                    <  1:331029/143‑1 (MQ=255)
                                                                                                                               cGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCCTTGAATGTAAGTTAAGGGGCATCAACTCTAAACTCAACCAGGATGTTTGACGCGGGCGACGTCGAGCACCAGCAGAATCaa             >  1:491064/1‑142 (MQ=255)
                                                                                                                                   ggaaCGGTATCAGCCACGTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTa         <  1:53276/141‑1 (MQ=255)
                                                                                                                                          tctcAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCa  <  2:412014/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
AAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGCAGAATCAGCTTATCTTCCA  >  NC_000913/4196110‑4196390

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑