Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,634,018 T→C V46A (GTA→GCA)  yigF → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,634,0180TC100.0% 266.4 / NA 78V46A (GTA→GCA) yigFconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/49);  total (29/49)

TCAACATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTG  >  W3110S.gb/3633950‑3634083
                                                                    |                                                                 
tCAACATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCaa                                                                 <  1:4645651/71‑1 (MQ=255)
   aCATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:2039924/71‑1 (MQ=255)
   aCATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:1935742/71‑1 (MQ=255)
   aCATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:3774461/71‑1 (MQ=255)
   aCATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:625222/71‑1 (MQ=255)
           ttCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       >  1:1102292/1‑70 (MQ=255)
           ttCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:3076721/70‑1 (MQ=255)
           ttCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:1561183/70‑1 (MQ=255)
           ttCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:429986/70‑1 (MQ=255)
            tCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTa                                                      <  1:2563417/70‑1 (MQ=255)
            tCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTa                                                      <  1:2795363/70‑1 (MQ=255)
                   ttCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCa                                              <  1:794783/71‑1 (MQ=255)
                   ttCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCa                                              <  1:1118822/71‑1 (MQ=255)
                     cTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTa                                                  >  1:3420490/1‑65 (MQ=255)
                        gggggCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:3225294/70‑1 (MQ=255)
                        gggggCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:1759737/70‑1 (MQ=255)
                        gggggCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:2256575/70‑1 (MQ=255)
                        gggggCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:4665158/70‑1 (MQ=255)
                        gggggCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:2455603/70‑1 (MQ=255)
                            gCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATg                                          <  1:678295/66‑1 (MQ=255)
                            gCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAAttt                                     <  1:124189/71‑1 (MQ=255)
                            gCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAAttt                                     <  1:3352719/71‑1 (MQ=255)
                             cAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGTATTTc                                    <  1:4621043/71‑1 (MQ=255)
                               acacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGAc                                                              <  1:2835274/43‑1 (MQ=255)
                               acacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:2897920/50‑1 (MQ=255)
                                cacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:225526/49‑1 (MQ=255)
                                cacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:615888/49‑1 (MQ=255)
                                cacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:1667775/49‑1 (MQ=255)
                                cacCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:1613805/49‑1 (MQ=255)
                                 acCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTg                                >  1:92069/1‑71 (MQ=255)
                                 acCAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTg                                >  1:3025607/1‑71 (MQ=255)
                                  ccAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTTTTCAAATGAATTTCATTGc                               <  1:35990/71‑1 (MQ=255)
                                  ccAGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGc                               >  1:1707753/1‑71 (MQ=255)
                                    aGAATATAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCtt                             >  1:2211476/1‑71 (MQ=255)
                                       atatAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTttct                          >  1:840479/1‑71 (MQ=255)
                                       atatAAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTttct                          >  1:699389/1‑71 (MQ=255)
                                        tataAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAGCTATTCAAATGAATTTCATTGCTttctt                         >  1:4239292/1‑71 (MQ=255)
                                         ataAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTttctt                         <  1:4111303/70‑1 (MQ=255)
                                          taAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCAtt                                 >  1:4224003/1‑61 (MQ=255)
                                          taAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTtcttct                       >  1:897052/1‑71 (MQ=255)
                                          taAAGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTtcttct                       >  1:696977/1‑71 (MQ=255)
                                            aaGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     >  1:1130207/1‑71 (MQ=255)
                                            aaGCAGCATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:780272/71‑1 (MQ=255)
                                             agcagcATTTAAAGCATTAAAAGCAAGACAACGCtt                                                       <  1:625975/36‑1 (MQ=255)
                                             agcagcATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:1479288/70‑1 (MQ=255)
                                             agcagcATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:3914722/70‑1 (MQ=255)
                                             agcagcATTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:3970418/70‑1 (MQ=255)
                                                   aTTTAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGa              <  1:4059673/71‑1 (MQ=255)
                                                    tttAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:548926/63‑1 (MQ=255)
                                                    tttAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:1825760/63‑1 (MQ=255)
                                                      tAAAGCATTAAACGCAAGACACCGCTTAACTATTCAAATGAATTTCATTGCt                              >  1:815092/1‑52 (MQ=37)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCAtt                                 >  1:1596682/1‑49 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTtcttct                       >  1:284416/1‑59 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1708758/1‑70 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:1952202/1‑70 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:4169427/1‑70 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:3013413/1‑70 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:379255/1‑70 (MQ=255)
                                                      tAAAGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGata            >  1:3127909/1‑70 (MQ=255)
                                                      tAAAGCATTAAAACCAAGACAACGCTTAACTATTCAAATGAATTTCAATGCTTTCTACTGTTCCt                 >  1:3558088/1‑65 (MQ=255)
                                                       aaaGCATTAAAGGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGatata          >  1:1800831/1‑71 (MQ=255)
                                                       aaaGCATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGatat           <  1:663978/70‑1 (MQ=255)
                                                           cATTAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtattta      >  1:3358953/1‑71 (MQ=255)
                                                             ttAAAAGCAAGACAACTCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtattta      <  1:2368863/69‑1 (MQ=255)
                                                             ttAAAAGCAAGACAACTCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATAtattta      <  1:2363580/69‑1 (MQ=255)
                                                             ttAAAAGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGt                     <  1:2340053/54‑1 (MQ=255)
                                                               aaaaGCAAGACAACGCTTAACTATTCAAATGAATTTCa                                   >  1:3476330/1‑38 (MQ=255)
                                                               aaaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTttctt                         >  1:2951621/1‑48 (MQ=255)
                                                               aaaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  >  1:55223/1‑71 (MQ=255)
                                                               aaaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  >  1:2799910/1‑71 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:3298332/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:2368277/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:3518432/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:541957/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:2171433/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:891246/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:2969043/70‑1 (MQ=255)
                                                                aaaGCAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTg  <  1:3897962/70‑1 (MQ=255)
                                                                    |                                                                 
TCAACATGGTTTTCCTGGATTCTGGGGGGCAACACCAGAATATAAAGCAGCATTTAAAGCATTAAAAGTAAGACAACGCTTAACTATTCAAATGAATTTCATTGCTTTCTTCTGTTCCTGGATATATTTATTTG  >  W3110S.gb/3633950‑3634083

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: