Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,113 Δ1 bp coding (372/396 nt) yhgKJ ← RNA phosphage cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,1130C.94.9% 308.0 / 9.8 80coding (372/396 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base . (44/30);  minor base G (2/2);  total (48/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCA  >  NC_000913/3554043‑3554183
                                                                       |                                                                       
ccTATTTTCGCCAGCAGCGGCTCCAGTACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                                   <  1:4394423/76‑1 (MQ=255)
 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                  <  1:1807443/76‑2 (MQ=255)
 cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                  <  1:4604619/76‑2 (MQ=255)
  tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                                                 <  1:2747906/76‑3 (MQ=255)
   aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGTCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                                                >  1:3859406/1‑73 (MQ=255)
       tCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                             <  1:4099206/76‑1 (MQ=255)
       tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                             >  1:3617561/1‑76 (MQ=255)
       tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCCACGGGGTATTGGTGGCGC‑CGCTccc                                                             >  1:475139/1‑74 (MQ=255)
        cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                            <  1:2630716/76‑1 (MQ=255)
        cGCCAGCAGCGGCTCCAGCACCCCGCGGAGAAAATCAGCAGGCGGCGCCCCCGGGTTATCGGTCCCCC‑CGCCCACt                                                            >  2:1718441/1‑76 (MQ=255)
         cccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                           <  1:700046/75‑1 (MQ=255)
           cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                                         >  1:2950094/1‑76 (MQ=255)
            agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATAGGT‑GCCGCCGCTCACTTCaa                                                        >  1:4850711/1‑76 (MQ=255)
              cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGGCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCCacc                                                      >  2:4316987/1‑75 (MQ=255)
              cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                                      >  1:4030622/1‑76 (MQ=255)
               agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAacac                                                     >  1:2387443/1‑76 (MQ=255)
               agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGGCCCCC‑CGCCCCCCTCAaaac                                                     >  2:916362/1‑76 (MQ=255)
                gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                                    <  2:3019909/76‑1 (MQ=255)
                  ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                                  >  1:3458170/1‑76 (MQ=255)
                   gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                                 >  1:1516797/1‑76 (MQ=255)
                     tCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCCC‑AGCTCACTGCAACACGCCAAg                                               >  2:3098578/1‑76 (MQ=255)
                     cccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                               <  1:3094290/75‑1 (MQ=255)
                       cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                             <  1:598005/76‑1 (MQ=255)
                       cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                             <  1:247615/76‑1 (MQ=255)
                        aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTc                                            >  1:3708937/1‑76 (MQ=255)
                         ggaCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                           <  2:2430666/74‑1 (MQ=255)
                         gCACCCGGCGGATAAAATCCGCGGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                           <  1:4245024/76‑1 (MQ=255)
                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCCACACGCGAAGGTcc                                           >  1:1152671/1‑76 (MQ=255)
                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                           <  1:2674116/76‑1 (MQ=255)
                         gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTcc                                           <  2:2585498/76‑1 (MQ=255)
                          cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                          >  2:2679099/1‑76 (MQ=255)
                          cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGGCCa                                          >  2:3695984/1‑76 (MQ=255)
                            cccGGCGGATAGAATCCGCGGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                        <  2:1149211/76‑1 (MQ=255)
                            cccGGCGGATAAAATCCGCGGGCGGCACCGAGGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                        <  1:601504/76‑1 (MQ=255)
                            cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                        <  2:1255810/76‑1 (MQ=255)
                            cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                        <  1:4690778/76‑1 (MQ=255)
                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                     >  1:706242/1‑76 (MQ=255)
                                gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                                    >  1:1905830/1‑76 (MQ=255)
                                gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                                    >  1:3994208/1‑76 (MQ=255)
                                gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                                    >  1:3364904/1‑76 (MQ=255)
                                 cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATAGGCaa                                   >  1:1313397/1‑76 (MQ=255)
                                 cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCACCGGCaa                                   >  2:3412423/1‑76 (MQ=255)
                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                  <  1:2835090/76‑1 (MQ=255)
                                  ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                  >  1:4496734/1‑76 (MQ=255)
                                   gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                 >  2:1637919/1‑76 (MQ=255)
                                    aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAACACGCGAAGGTCCATCGGCACAcc                                >  1:4762489/1‑76 (MQ=255)
                                    aTAAAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                <  2:4287181/76‑1 (MQ=255)
                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                              >  2:4113008/1‑76 (MQ=255)
                                      aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGCCCATCGGCAAACcac                              >  2:1161618/1‑76 (MQ=255)
                                            ctcaGGCGGCGCCGACTGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                        <  1:280840/73‑1 (MQ=255)
                                            cgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                        <  1:1550361/73‑1 (MQ=255)
                                             gcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                       >  1:2969177/3‑76 (MQ=255)
                                               aGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                     >  2:1694704/1‑76 (MQ=255)
                                                  gggcgCCGGCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagc                  <  2:3729118/75‑1 (MQ=255)
                                                  cggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCCTCGGCAAACCACAGCGCGGgccgc                  >  2:1586912/1‑76 (MQ=255)
                                                   ggcgCCGGCGGGTTATCGG‑GGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                 <  2:4672949/76‑1 (MQ=255)
                                                   ggcgCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCCTCGGCAAACCACAGCGCGGccagcc                 >  2:1534720/1‑75 (MQ=255)
                                                     cgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc               <  2:1516796/76‑1 (MQ=255)
                                                     cgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc               <  1:2560613/76‑1 (MQ=255)
                                                     cgcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc               <  1:2523800/76‑1 (MQ=255)
                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCCGCAccc              >  2:3842954/1‑75 (MQ=255)
                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg              >  1:2330399/1‑76 (MQ=255)
                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg              >  2:1851063/1‑76 (MQ=255)
                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg              >  1:176520/1‑76 (MQ=255)
                                                      gcCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg              >  1:2430744/1‑76 (MQ=255)
                                                       ccGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCCCCGt             >  2:2603160/1‑76 (MQ=255)
                                                       ccGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt             >  2:2617655/1‑76 (MQ=255)
                                                       ccGACGGGTGATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt             <  2:939794/76‑1 (MQ=255)
                                                        cGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGCGCCGCACCGtt            >  2:3250702/1‑76 (MQ=255)
                                                         gACGGGTTATCGGT‑GCCGCCGCTCATTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt           >  1:3957235/1‑76 (MQ=255)
                                                         gACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt           <  1:882252/76‑1 (MQ=255)
                                                         gACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt           >  1:2486988/1‑76 (MQ=255)
                                                         gACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt           <  2:1988943/76‑1 (MQ=255)
                                                         gACGGGTTATCGGT‑GCCGCCGATCACTTCAACACGCGAAGGACCATAGGCAAACCACAGCGACGGCAGCACCGatt           >  1:3015646/1‑76 (MQ=255)
                                                           ggggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgc         <  2:1463133/75‑1 (MQ=255)
                                                           cGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgc         >  1:726229/1‑76 (MQ=255)
                                                            gggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca        <  1:663137/76‑1 (MQ=255)
                                                             ggTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGCTTgcag       >  2:3072149/1‑76 (MQ=255)
                                                              gTTATCGGT‑GCCGCCGCTCACTTCAACACGCGGAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc      <  1:2718866/76‑1 (MQ=255)
                                                               ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca     <  2:1684167/76‑1 (MQ=255)
                                                               ttATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca     >  1:1055129/1‑76 (MQ=255)
                                                               ttATCGGT‑GCCGCCGATCACTTCCACACGCGAAGGTCCATCGGGAAACCACAGCGCGGGCAGCACCGTTTgcagca     >  1:3674348/1‑76 (MQ=255)
                                                                 aTCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc   >  2:2243152/1‑76 (MQ=255)
                                                                 aTCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc   >  2:4371008/1‑76 (MQ=255)
                                                                  tCGGTGCCCC‑CGCGCACCTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa  <  2:91153/76‑1 (MQ=255)
                                                                       |                                                                       
CCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCA  >  NC_000913/3554043‑3554183

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