Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,027 C→G intergenic (‑151/+58) ypjA ← / ← b2648 putative ATP‑binding component of a transport system/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,781,0270CG98.2% 137.5 / ‑6.1 57intergenic (‑151/+58)ypjA/b2648putative ATP‑binding component of a transport system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (27/29);  minor base A (1/0);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTT  >  NC_000913/2780955‑2781098
                                                                         |                                                                       
tgtgGTTGGATAAAAAAGTCAAGCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCa                                                                       >  2:4130688/1‑76 (MQ=38)
 gtgGTTGGATAAAAAAGTCAATCCATTCAGGAACTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaa                                                                      <  1:4709750/76‑1 (MQ=38)
  tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                     >  2:1005755/1‑76 (MQ=38)
  tgGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaa                                                                     <  2:2633455/76‑1 (MQ=38)
   gtttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                    >  1:3688103/3‑76 (MQ=38)
   ggTTGGATAGAAAAGTCAAGCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                    <  2:1875448/76‑1 (MQ=37)
   ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                    <  1:3367493/76‑1 (MQ=38)
   ggTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTAGTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaa                                                                    <  2:1514314/76‑1 (MQ=37)
    gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                   <  1:1370716/76‑1 (MQ=38)
    gTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCaaaat                                                                   >  2:2024836/1‑76 (MQ=38)
     ttGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTCACaaaata                                                                  >  2:4779150/1‑76 (MQ=255)
          tAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATg                                                             <  2:1514801/76‑1 (MQ=37)
           aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                            <  1:2352269/76‑1 (MQ=37)
           aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                            <  1:2586844/76‑1 (MQ=37)
           aaaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGa                                                            >  1:1664886/1‑76 (MQ=37)
            aaaaaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGAt                                                           >  2:700313/1‑76 (MQ=37)
               aaGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTa                                                        >  2:150108/1‑76 (MQ=25)
                 gTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACa                                                      <  1:4317170/76‑1 (MQ=37)
                   cAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                    <  2:4567656/76‑1 (MQ=37)
                   cAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                    <  1:1621521/76‑1 (MQ=37)
                   cAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATa                                                    <  1:4062527/76‑1 (MQ=37)
                    aaTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATaa                                                   >  2:992557/1‑76 (MQ=37)
                     aTCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAt                                                  <  2:2647703/76‑1 (MQ=37)
                     aTCCATTCAGGAAATACGGGCGGATTCTTTTCTTTCGACAATGAGGCCGGTGGCAAAATAAAATGATTTACATAAt                                                  <  2:2657467/76‑1 (MQ=25)
                      tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                 >  1:1425451/1‑76 (MQ=37)
                      tCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATc                                                 >  1:4759706/1‑76 (MQ=37)
                       ccATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCg                                                <  1:4269728/76‑1 (MQ=37)
                        cATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGt                                               >  2:3304146/1‑76 (MQ=37)
                          ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                             <  2:357527/76‑1 (MQ=37)
                          ttCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGttt                                             <  2:3443548/76‑1 (MQ=37)
                          ttCAGGAAATACGGGAGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAAGCGgtt                                             >  1:794848/1‑76 (MQ=21)
                              ggAAATACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtga                                         <  1:2491662/76‑1 (MQ=25)
                                aaaTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCtgatg                                       <  1:1011427/76‑1 (MQ=25)
                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                     <  2:4136693/76‑1 (MQ=25)
                                  aTACGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGaa                                     >  1:2721064/1‑76 (MQ=25)
                                   ttcGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAat                                    <  1:241476/74‑1 (MQ=21)
                                     cGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                  >  1:2025962/1‑76 (MQ=25)
                                     cGGGCGTATTCTTTTCTTTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAatat                                  <  1:3468285/76‑1 (MQ=25)
                                     cGGGCGTATTCTTGTCTTTCGACAATGAGGCCGTGGGCAAAATAAAAGGATTTACATAATCGTTTCTGATGAatat                                  <  1:3490809/76‑1 (MQ=18)
                                             ttcttttcttTCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCt                          >  1:1767398/1‑76 (MQ=18)
                                              tcttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTc                         >  1:1968993/1‑76 (MQ=18)
                                               cttttctttCGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATAGTTTCTGATGAATATCTTCTGCTca                        >  1:625621/1‑76 (MQ=18)
                                                 tttctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTcaca                      <  1:2152662/76‑1 (MQ=17)
                                                    ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                   >  2:2610052/1‑76 (MQ=17)
                                                    ctttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaa                   >  2:2592587/1‑76 (MQ=17)
                                                     tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaa                  >  1:2292511/1‑76 (MQ=17)
                                                     tttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaa                  >  2:1843970/1‑76 (MQ=17)
                                                      ttcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaa                 >  1:4113533/1‑76 (MQ=17)
                                                       tcGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATaaaaa                >  2:3897203/1‑76 (MQ=17)
                                                        cGACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAAt               >  2:2436036/1‑76 (MQ=17)
                                                         gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc              <  1:4436922/76‑1 (MQ=17)
                                                         gACAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATc              <  2:940587/76‑1 (MQ=17)
                                                          aCAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATca             <  2:3076437/76‑1 (MQ=17)
                                                           cAATGAGGCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATcac            >  2:2961540/1‑76 (MQ=17)
                                                                aggCCGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAAt       >  2:4761082/3‑76 (MQ=255)
                                                                   gcGTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAAc    <  1:1518321/75‑1 (MQ=255)
                                                                     gTTGGCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACtt  >  2:4605441/1‑76 (MQ=255)
                                                                         |                                                                       
TGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGTATTCTTTTCTTTCGACAATGAGCCGTTGCCAAAATAAAATGATTTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTT  >  NC_000913/2780955‑2781098

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