Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,342,922 C→T R89C (CGT→TGT)  yciH → putative translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,342,9220CT90.8% 206.8 / 10.5 65R89C (CGT→TGT) yciHputative translation factor
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (30/29);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

CTGCGGCGGAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGCGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTAA  >  NC_000913/1342873‑1342967
                                                 |                                             
ctgcggcGGAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGc                                               >  1:5778966/1‑50 (MQ=255)
  gcggcgGAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAgtgt                                             >  1:2218749/1‑50 (MQ=255)
   cggcggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAgtgt                                             >  1:1776530/1‑49 (MQ=255)
    ggcggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGa                                           <  1:2651714/50‑1 (MQ=255)
      cggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGAtt                                         <  1:5696612/50‑1 (MQ=255)
       ggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGAttt                                        <  1:2706157/50‑1 (MQ=255)
       ggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGAttt                                        <  1:5030356/50‑1 (MQ=255)
       ggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGAtt                                         >  1:5653620/1‑49 (MQ=255)
       ggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGAtt                                         >  1:1931879/1‑49 (MQ=255)
       ggAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGCGTGAtt                                         >  1:5385503/1‑49 (MQ=255)
        gAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATtta                                       <  1:4282545/50‑1 (MQ=255)
        gAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATtta                                       >  1:1530054/1‑50 (MQ=255)
        gAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATtta                                       <  1:4531881/50‑1 (MQ=255)
        gAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATtta                                       >  1:3752579/1‑50 (MQ=255)
         aGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttat                                      <  1:1249224/50‑1 (MQ=255)
         aGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttat                                      <  1:3814590/50‑1 (MQ=255)
         aGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttat                                      >  1:2343643/1‑50 (MQ=255)
         aGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGCGTGATttat                                      >  1:757046/1‑50 (MQ=255)
          gCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttatt                                     <  1:5760798/50‑1 (MQ=255)
          gCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttatt                                     <  1:3801822/50‑1 (MQ=255)
           cAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttatta                                    <  1:630484/50‑1 (MQ=255)
            aGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATttatta                                    >  1:674632/1‑49 (MQ=255)
             gTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTaa                                   >  1:4819087/1‑49 (MQ=255)
              ttAAAGATGAAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTaaaa                                 <  1:511174/50‑1 (MQ=255)
                aaaGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAAt                                <  1:4440983/49‑1 (MQ=255)
                 aaGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCa                              >  1:4431857/1‑50 (MQ=255)
                  aGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCac                             >  1:1975354/1‑50 (MQ=255)
                  aGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCac                             >  1:1620733/1‑50 (MQ=255)
                  aGATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCac                             >  1:2873505/1‑50 (MQ=255)
                  aGATGGAGTTATTGAAATCCAGGGCGATAAGCGTGATTTATTAAAATCac                             <  1:2508939/50‑1 (MQ=255)
                   gATGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCac                             >  1:940026/1‑49 (MQ=255)
                    aTGGAGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCacta                           <  1:1356500/50‑1 (MQ=255)
                    aTGGAGTTATTGAAATCCAGGGCGATAAGCGTGATTTATTAAAATCacta                           >  1:1276515/1‑50 (MQ=255)
                     tGGAGTTATTGAAATCCAGGGCGATAAGCGTGATTTATTAAAATCactac                          <  1:234630/50‑1 (MQ=255)
                        aGTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTg                       <  1:3558921/50‑1 (MQ=255)
                         gTTATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGa                      <  1:315044/50‑1 (MQ=255)
                          ttATTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTg                       >  1:3229723/1‑48 (MQ=255)
                            aTTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAg                    <  1:897001/49‑1 (MQ=255)
                            aTTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAg                    <  1:5714481/49‑1 (MQ=255)
                            aTTGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGc                   >  1:2760339/1‑50 (MQ=255)
                             ttGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGc                   <  1:1241787/49‑1 (MQ=255)
                             ttGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGc                   >  1:197731/1‑49 (MQ=255)
                             ttGAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGc                   <  1:3328140/49‑1 (MQ=255)
                               gAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGa                 <  1:5575799/49‑1 (MQ=255)
                               gAAATCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGa                 >  1:132587/1‑49 (MQ=255)
                                aaaTCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGaa                <  1:6712033/49‑1 (MQ=255)
                                  aTCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAgg             >  1:4674549/1‑50 (MQ=255)
                                   tCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAgg             >  1:6301530/1‑49 (MQ=255)
                                   tCCAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAgg             >  1:2254767/1‑49 (MQ=255)
                                    ccAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGa           >  1:3045014/1‑50 (MQ=255)
                                     cAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          >  1:3719444/1‑50 (MQ=255)
                                     cAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          <  1:2767042/50‑1 (MQ=255)
                                     cAGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          >  1:5308982/1‑50 (MQ=255)
                                     cAGGGCGATAAGCGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          >  1:5814000/1‑50 (MQ=255)
                                      aGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          >  1:5409451/1‑49 (MQ=255)
                                      aGGGCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGAt          >  1:1267134/1‑49 (MQ=255)
                                       gggCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATg         >  1:5502097/1‑49 (MQ=255)
                                        ggCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGaa       <  1:1621076/50‑1 (MQ=255)
                                         gCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAg      >  1:563896/1‑50 (MQ=255)
                                         gCGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAg      <  1:1885112/50‑1 (MQ=255)
                                          cGATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAgg     >  1:124841/1‑50 (MQ=255)
                                           gATAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGt    <  1:693411/50‑1 (MQ=255)
                                            aTAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGt    >  1:5686744/1‑49 (MQ=255)
                                            aTAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTa   <  1:909027/50‑1 (MQ=255)
                                            aTAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTa   <  1:1702667/50‑1 (MQ=255)
                                             tAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTaa  <  1:1995142/50‑1 (MQ=255)
                                             tAAGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTa   >  1:5680462/1‑49 (MQ=255)
                                              aaGTGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTaa  >  1:1219839/1‑49 (MQ=255)
                                                 |                                             
CTGCGGCGGAGCAGTTAAAGATGGAGTTATTGAAATCCAGGGCGATAAGCGTGATTTATTAAAATCACTACTTGAAGCGAAAGGGATGAAGGTAA  >  NC_000913/1342873‑1342967

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