Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,477,041 Δ1 bp intergenic (‑151/‑269) b4256 ← / → yjgN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,477,0410G.98.8% 361.9 / ‑3.1 81intergenic (‑151/‑269)b4256/yjgNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (33/47);  total (33/48)
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.99e-01

GACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCC  >  NC_000913/4476972‑4477110
                                                                     |                                                                       
gACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTa                                                                  >  1:1651177/1‑76 (MQ=255)
    gtgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCt                                                              <  1:702440/76‑1 (MQ=255)
     tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGGCGCGGGCGATAGCAGGATTATCCTg                                                             >  2:2464436/1‑76 (MQ=255)
     tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTg                                                             <  1:3636107/76‑1 (MQ=255)
     tgtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTg                                                             <  1:1590341/76‑1 (MQ=255)
      gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCCGGATTATCCTGc                                                            >  2:4321434/1‑76 (MQ=255)
      gtAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGc                                                            <  1:1847279/76‑1 (MQ=255)
       tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCa                                                           <  2:1972087/76‑1 (MQ=255)
       tAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCa                                                           <  2:4601004/76‑1 (MQ=255)
         gCCTTTATCAGCGGTAAGCCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGc                                                         <  2:2368793/76‑1 (MQ=255)
         gCCGTTGTCAGCGGTGAGACCGTATTCGGCGGATGCCTGACGCATGACGCGGGCGATAGCAGGCTTATCCTGCAGc                                                         <  2:808659/76‑1 (MQ=255)
           cTTTATCAGCGGTAAGACCGTATGCGGCGGATACCTGACGAATGCCGCGGGCGATCTCAGGCTTATCCTGCAGCGt                                                       <  1:556254/76‑1 (MQ=255)
              tATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                    <  1:4321434/76‑1 (MQ=255)
              tATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAg                                                    <  2:2382660/76‑1 (MQ=255)
               aTCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAgg                                                   <  1:2504283/76‑1 (MQ=255)
                tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                  >  2:3837435/1‑76 (MQ=255)
                tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                  >  2:3108474/1‑76 (MQ=255)
                tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                  >  1:1057778/1‑76 (MQ=255)
                tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                  <  2:4233097/76‑1 (MQ=255)
                tCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGc                                                  >  1:1533975/1‑76 (MQ=255)
                 cAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcg                                                 <  1:2353991/76‑1 (MQ=255)
                   gCGGTGAGGCCGTGTGCGGCGGATGCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCGGCGTCAGGcgac                                               <  2:2653779/76‑1 (MQ=255)
                   gCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgac                                               <  1:1556703/76‑1 (MQ=255)
                   gCGGGAAGGCCGTGTTCGGGGGGTGCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgac                                               <  2:332780/76‑1 (MQ=255)
                    cGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacg                                             >  1:3437306/1‑76 (MQ=255)
                     ggTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGcgacgc                                            >  1:3004857/1‑75 (MQ=255)
                          gACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTAc                                       >  1:1564502/1‑76 (MQ=255)
                           aCCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACa                                      <  1:2257034/76‑1 (MQ=255)
                            ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                     <  1:3481326/76‑1 (MQ=255)
                            ccGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAg                                     <  2:1761506/76‑1 (MQ=255)
                             cggaTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAgg                                    <  2:3704070/73‑1 (MQ=255)
                             cGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGTGACGCATTACAgg                                    <  2:4023304/76‑1 (MQ=255)
                              gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAggg                                   >  1:634606/1‑75 (MQ=255)
                              gTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGt                                   >  2:2972790/1‑76 (MQ=255)
                                aTTCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGa                                 >  2:3960893/1‑76 (MQ=255)
                                  tCGGCGGATACCTGGCGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                               <  2:446383/76‑1 (MQ=255)
                                  tCGGCGGATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGAtt                               <  1:3837435/76‑1 (MQ=255)
                                   cggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTg                              <  2:3024876/76‑1 (MQ=255)
                                    ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             >  2:2343336/1‑76 (MQ=255)
                                    ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             <  1:3298575/76‑1 (MQ=255)
                                    ggcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgc                             <  2:1478775/76‑1 (MQ=255)
                                     gcggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcg                            >  1:1747763/1‑76 (MQ=255)
                                      cggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                           <  1:1159088/76‑1 (MQ=255)
                                      cggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcgg                           >  1:875595/1‑76 (MQ=255)
                                       ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                          <  2:4809603/76‑1 (MQ=255)
                                       ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTgcggc                          <  2:3605509/76‑1 (MQ=255)
                                       ggATACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTccggc                          >  2:176809/1‑76 (MQ=255)
                                         aTACCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGgcgc                        <  1:2685434/76‑1 (MQ=255)
                                           aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                      >  1:393854/1‑76 (MQ=255)
                                           aCCTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCaa                      >  1:2972853/1‑76 (MQ=255)
                                            ccTGACGAGTCAGGCGGGCGAGAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGATGATTGCGGCGCAAt                     <  1:3320466/76‑1 (MQ=255)
                                             cTGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATg                    >  2:2731245/1‑76 (MQ=255)
                                              tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGTGACGCATTACAGGTGATTGCGGCGCAATGt                   >  2:3424175/1‑76 (MQ=255)
                                              tGACGAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGt                   <  2:2725327/76‑1 (MQ=255)
                                                  gactGACGCGGGCGATAGGTGGATTATCCTGCAGCGTCAGGCGACGCATTCCAGGAGATTGCGGCGCAATGTTCtt               <  2:3000782/73‑1 (MQ=255)
                                                  gAATGGCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt               <  2:4700800/76‑1 (MQ=255)
                                                  gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt               <  2:1636352/76‑1 (MQ=255)
                                                  gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt               <  1:2963121/76‑1 (MQ=255)
                                                  gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt               >  1:1819241/1‑76 (MQ=255)
                                                  gAATGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTAtt               >  1:4534887/1‑76 (MQ=255)
                                                   aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc              <  1:2576690/76‑1 (MQ=255)
                                                   aaTGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTc              >  2:2520471/1‑76 (MQ=255)
                                                    aTGACGCGGGCGACAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCa             >  2:1396884/1‑76 (MQ=255)
                                                     tGACGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCAt            >  2:829264/1‑76 (MQ=255)
                                                       aCGCGGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAg          <  2:2572941/76‑1 (MQ=255)
                                                        cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  2:3088664/76‑1 (MQ=255)
                                                        cgcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGt         <  1:3656851/76‑1 (MQ=255)
                                                         gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt        >  2:3755647/1‑76 (MQ=255)
                                                         gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt        >  1:3697519/1‑76 (MQ=255)
                                                         gcgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGtt        >  1:1711655/1‑76 (MQ=255)
                                                          cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt       >  2:706424/1‑76 (MQ=255)
                                                          cgGGCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGttt       <  1:3367764/76‑1 (MQ=255)
                                                          catgCGATAGCAGGATTATCCGGCAGCGTCAGGCGACGCATGACAGGTGATTGCGGCGCAAGGTTATTCATAGttt       <  1:3392416/73‑1 (MQ=255)
                                                           gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa      >  1:1845650/1‑76 (MQ=255)
                                                           gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa      >  1:2323139/1‑76 (MQ=255)
                                                           gggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTa      <  2:3517267/76‑1 (MQ=255)
                                                            ggCGATAGCAGGGTTATCCTGCAGCGTCAGGCGACGGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa     <  2:4580832/76‑1 (MQ=255)
                                                            ggCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTaa     >  2:248973/1‑76 (MQ=255)
                                                             gCGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAc    <  1:1564600/76‑1 (MQ=255)
                                                              cGATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTGCAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAcc   <  1:3843536/76‑1 (MQ=255)
                                                               gATAGCAGGATTATCCTGCAGCGTCAGGCGACGCATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAAccc  <  1:537851/76‑1 (MQ=255)
                                                                     |                                                                       
GACGGTGTAGCCTTTATCAGCGGTAAGACCGTATTCGGCGGATACCTGACGAATGACGCGGGCGATAGCGAGGATTATCCTGCAGCGTCAGGCGACGATTACAGGTGATTGCGGCGCAATGTTATTCATAGTTTAACCC  >  NC_000913/4476972‑4477110

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