Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 475,972 C→G 14.6% intergenic (+21/+10) ybaY → / ← atl PF09619 family lipoprotein YbaY/DNA base‑flipping protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913475,9720CG14.6% 109.7 / 12.8 47intergenic (+21/+10)ybaY/atlPF09619 family lipoprotein YbaY/DNA base‑flipping protein
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base G (5/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

TACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAA  >  NC_000913/475836‑476112
                                                                                                                                        |                                                                                                                                            
tACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGgaga                                                                                                                                        <  1:792016/143‑1 (MQ=255)
  cAGCAAACCGCCGGGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgg                                                                                                                                      <  1:156807/143‑1 (MQ=255)
   aggaacccgccgTGCCGGTTCAGGCCCGCGGGGGCGCAACGCCTACCGTCCCTTCGACTTCCCCAACTCAGGTGACTCCGTCTTCGGCAGTTTCAGCTCCTACGGAATATTGGGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                     <  2:248766/137‑1 (MQ=255)
   aGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAggg                                                                                                                                     <  2:434745/143‑1 (MQ=255)
    gCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGCATCCGTCTTCGGCAGTTACAGCTCCTACGCAACATTAAGCGACGTTAACCCTCTCCGACCGGCGAGGGt                                                                                                                                    >  2:585108/1‑143 (MQ=255)
       aaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAgggccag                                                                                                                                 >  1:143345/1‑139 (MQ=255)
         acgccgTGCCGTCTGAGGCCAGCGGGGGCGCAACGACTACCGGACCTTCGGCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAGTATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagta                                                                                                                               <  1:164476/142‑1 (MQ=255)
                            cAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTACGTCGGTTCCAGATCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGc                                                                                                            >  1:636986/1‑143 (MQ=255)
                                   ggCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGa                                                                                                     <  1:547038/143‑1 (MQ=255)
                                      gcAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGCGGGTTAGTAGTTCCAGCGATGACGCGGCAAGTCGAttt                                                                                                  >  1:276596/1‑143 (MQ=255)
                                       cAACGACTACCGTACCTTCGACTTCACCACCTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTg                                                                                                 <  2:347973/143‑1 (MQ=255)
                                            aCTACCGTCCCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGc                                                                                            <  2:268394/143‑1 (MQ=255)
                                              tACCGTCCCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGcct                                                                                          <  1:529655/143‑3 (MQ=255)
                                                ccGTCCCTTCGACTTCCCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTcc                                                                                        <  1:800334/143‑1 (MQ=255)
                                                ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCAGCTTcc                                                                                        >  2:229830/1‑143 (MQ=255)
                                                       ttggACCTCACCAACTCAGGTGAATCCGTCTGCGGCGGTTTCAGCGCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                 <  2:90685/140‑1 (MQ=255)
                                                       ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTTCAGATTCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCTGGTGAGGGTTAGTATTTTCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTACCGATAcc                                                                                 >  1:194177/1‑143 (MQ=255)
                                                       ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCCGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGCGGGTTCGTAGTTCCAGCGATACCCCTGCAAGTCGGTGTGCCCGCGTCCCGACAcc                                                                                 >  2:334283/1‑143 (MQ=255)
                                                       ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATAcc                                                                                 >  1:351000/1‑143 (MQ=255)
                                                         cgccTTCACCAACTCAGGTGAAGCCGTCCTCCGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTGAACTCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                               <  1:45233/140‑1 (MQ=255)
                                                         cGACTTCACCCACTCCGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCAt                                                                               <  1:59975/143‑1 (MQ=255)
                                                                aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCt                                                                        >  1:39089/1‑143 (MQ=255)
                                                                 ccAACTCAGGTGAATCCGTCTTCGGCAGTTCCATGTCCTACGCAAGATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                       >  1:346357/1‑143 (MQ=255)
                                                                     cTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                   >  1:546388/1‑143 (MQ=255)
                                                                        aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACTCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAg                                                                <  1:380405/143‑1 (MQ=255)
                                                                                   tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGCTTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTACTGCCTGTc                                                     >  2:643693/1‑143 (MQ=255)
                                                                                   tCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGCTACCACCACACCTTCTGCCAGACATTCCTGta                                                     >  1:145447/1‑142 (MQ=255)
                                                                                          cAGTTCCACCTCCTACGCAATATTAAGCGACGCCACCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACg                                              <  2:596016/143‑1 (MQ=255)
                                                                                               ccACCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTa                                         <  1:802952/143‑1 (MQ=255)
                                                                                                cAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCCGTCGCTGACGCTGTaa                                        >  1:504620/1‑143 (MQ=255)
                                                                                                     ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCACGTCGATTTGCCCGCTTCCCGATACCATCGCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccg                                   <  2:215812/143‑1 (MQ=255)
                                                                                                      cTACGCGAGATTGGGCGACGTGAACACTCGCCGCCCGGAGAGGGTTAGTTGGTCCAGCGACAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccgg                                  <  1:374284/143‑1 (MQ=255)
                                                                                                          gCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccg                              >  1:387221/1‑143 (MQ=255)
                                                                                                             gtagtAAGCGACGTTAACCCGCGCCGCGCGGAGAGGGTTAGTAGTTCCAGCTATAACGCTGCAAGTCGATGTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGtt                           <  2:634748/139‑1 (MQ=255)
                                                                                                             atatTGAGCGACGTTAACGCGCTCCGCCGGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGtt                           <  1:318264/143‑1 (MQ=255)
                                                                                                                   aGCGTCGTTAACCCTCTCCGCCCGGAGGGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGaa                     <  1:793540/143‑1 (MQ=255)
                                                                                                                       aCGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAATTTCGGTCCGGTTAGCGTAATtg                 >  1:214597/1‑143 (MQ=255)
                                                                                                                        cGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGACAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATtgt                <  1:788089/143‑1 (MQ=255)
                                                                                                                            aaCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCg            >  2:95967/1‑143 (MQ=255)
                                                                                                                             aCCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGCTACCATCACACCTTCTGCCAGTACTGCCTGCCGCTCACGCTGTAAACCCGGCCCGGCTGGCGAGATCGCGccct           >  2:466021/1‑141 (MQ=255)
                                                                                                                              cccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATGCCATCACACCTTCTGCCAGTAATGTGTGTGGGGGACCCCGCAACTACGGTCTTCGTGGCGAAATTGTGCCGTg          >  1:806793/1‑143 (MQ=255)
                                                                                                                              cccTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAACTTGTGCGGTg          >  1:143945/1‑143 (MQ=255)
                                                                                                                                ctctCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTCAGCGAAATTGTGCCGTGGc        >  1:485848/1‑143 (MQ=255)
                                                                                                                                ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGc        <  2:346357/143‑1 (MQ=255)
                                                                                                                                  ctcAGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa      <  2:276596/143‑1 (MQ=255)
                                                                                                                                   tttgCCCGGAGAGGGTGAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAACTCCGGTCCGGTTAGCGAAATTGTGCCGTGCCGAt     >  2:670388/4‑143 (MQ=255)
                                                                                                                                      gCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTaa  <  1:332546/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
TACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAA  >  NC_000913/475836‑476112

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: ‑