Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,275,144 C→T intergenic (‑3/+197) yjcH ← / ← acs DUF485 family inner membrane protein/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,275,1440CT98.2% 189.6 / ‑4.6 57intergenic (‑3/+197)yjcH/acsDUF485 family inner membrane protein/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/30);  total (26/31)
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.14e-01

TTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAC  >  CP009273/4275079‑4275207
                                                                 |                                                               
tttttCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAg                                                               <  1:3830927/68‑1 (MQ=255)
   ttCGACTAACTCACTGAAATGCGCATAGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGAtt                                                            >  1:3199960/1‑68 (MQ=255)
     cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATtct                                                          <  1:834237/68‑1 (MQ=255)
     cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATtct                                                          <  1:79444/68‑1 (MQ=255)
     cGACTAACACCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATtct                                                          <  1:3325835/68‑1 (MQ=255)
     cGACAAACACCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATtct                                                          <  1:4424507/68‑1 (MQ=255)
          aaCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCaga                                                     <  1:3649035/68‑1 (MQ=255)
          aaCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCaga                                                     <  1:3765369/68‑1 (MQ=255)
           aCTCCCTGTAATGCGCATTGTCTTCTATCCGCTGATATATAGTGCCATTCATCATAGCTTCTCCagag                                                    <  1:4098233/68‑1 (MQ=255)
           aCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCagag                                                    <  1:3252966/68‑1 (MQ=255)
           aCTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCagag                                                    <  1:1708054/68‑1 (MQ=255)
           aCTCCCTGAAATGCGCATTGTCATCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCagag                                                    <  1:3903827/68‑1 (MQ=255)
             tCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGt                                                  >  1:297972/1‑68 (MQ=255)
                cTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAAt                                               <  1:693504/68‑1 (MQ=255)
                 tGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATg                                              <  1:318615/68‑1 (MQ=255)
                 tGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATg                                              <  1:2505873/68‑1 (MQ=255)
                  gAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGt                                             >  1:1266242/1‑68 (MQ=255)
                   aaaTGCGCATTGTCTTCTATCCGCTGATATATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTa                                            >  1:2563631/1‑68 (MQ=255)
                   aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTATTGTa                                            >  1:1756412/1‑68 (MQ=255)
                   aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTa                                            <  1:213114/68‑1 (MQ=255)
                   aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTa                                            >  1:1762436/1‑68 (MQ=255)
                         gcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGAtt                                      >  1:2810207/1‑68 (MQ=255)
                         gcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGAtt                                      >  1:3498627/1‑68 (MQ=255)
                           aTTGTCTTCTATCCGCTTATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGt                                    <  1:1304755/68‑1 (MQ=255)
                           aTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGt                                    >  1:2914814/1‑68 (MQ=255)
                            ttGTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGtt                                   <  1:1017809/68‑1 (MQ=255)
                            tagTCTTCTATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGAATGtt                                   <  1:4268155/66‑1 (MQ=255)
                                 ttctATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAt                              >  1:3296622/1‑68 (MQ=255)
                                  tctATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAtt                             >  1:1214398/1‑68 (MQ=255)
                                   ctATCCGCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTc                            >  1:285290/1‑68 (MQ=255)
                                         gCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTc                      >  1:442331/1‑68 (MQ=255)
                                         gCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTc                      <  1:923754/68‑1 (MQ=255)
                                         gCTGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTc                      <  1:1097412/68‑1 (MQ=255)
                                           tGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                    <  1:449062/68‑1 (MQ=255)
                                           tGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                    <  1:1380969/68‑1 (MQ=255)
                                           tGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                    <  1:3602730/68‑1 (MQ=255)
                                           tGATAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                    <  1:2356512/68‑1 (MQ=255)
                                           tGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                    <  1:1257893/68‑1 (MQ=255)
                                           ggATAAATAGTGCCATTCATCATAGATTCTCCAGCGGTAATGTAGGGCTTGTTTTAATTCCCGCTccc                    <  1:2320854/67‑1 (MQ=255)
                                             aTAAATAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCtt                  <  1:2135603/68‑1 (MQ=255)
                                                 aTAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATgg              >  1:2859621/1‑68 (MQ=255)
                                                 aTAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATgg              >  1:1103960/1‑68 (MQ=255)
                                                 aTAGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCACCCTTATgg              >  1:314977/1‑68 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCGTTATGGga            >  1:47648/1‑68 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATgggc            >  1:2558272/1‑67 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            >  1:4173845/1‑68 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            <  1:1201710/68‑1 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            >  1:3350949/1‑68 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            >  1:1377883/1‑68 (MQ=255)
                                                   aGTGCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            >  1:2954420/1‑68 (MQ=255)
                                                   aGTGACATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGga            >  1:2686789/1‑68 (MQ=255)
                                                      gCCATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGaa         >  1:4423833/1‑68 (MQ=255)
                                                       ccATTCATCATAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAg        <  1:1457071/68‑1 (MQ=255)
                                                           tcatcaTAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa    >  1:2969652/1‑68 (MQ=255)
                                                           tcatcaTAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa    >  1:4082799/1‑68 (MQ=255)
                                                             atcaTAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAc  <  1:4457256/68‑1 (MQ=255)
                                                             atcaTAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAc  <  1:3058107/68‑1 (MQ=255)
                                                                 |                                                               
TTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAC  >  CP009273/4275079‑4275207

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