Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	488,810	(A)5→4	coding (3051/3363 nt)	kefA →	fused mechanosensitive channel proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	488,806	0	A	.	100.0%	316.3 / NA	77	coding (3047/3363 nt)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (38/39);  total (38/39)

ACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAAC  >  W3110S.gb/488737‑488865                                                                      |                                                            aCTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaa                                                            <  1:2799800/71‑1 (MQ=255) aCTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaa                                                            <  1:3617237/71‑1 (MQ=255)  cTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaa                                                           <  1:3429588/71‑1 (MQ=255)  cTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGACTCGGCGTGGCCTATGGCTCCGATCTGGaa                                                            >  1:355534/1‑70 (MQ=255)      tCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                       <  1:1160351/71‑4 (MQ=255)        gTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaaagtg                                                       <  1:2900764/69‑4 (MQ=255)           gactgatACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:2119711/70‑1 (MQ=255)           gactgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:4628203/70‑1 (MQ=255)           gactgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:3186803/70‑1 (MQ=255)           gactgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:2982761/70‑1 (MQ=255)            actgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTaa                                                 <  1:676331/70‑1 (MQ=255)              tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgt                                              >  1:4038651/1‑71 (MQ=255)              tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgt                                              >  1:1638303/1‑71 (MQ=255)              tgacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgt                                              >  1:2418299/1‑71 (MQ=255)               gacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaa                                                            <  1:3976734/57‑1 (MQ=255)               gacACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaa                                                            <  1:3228805/57‑1 (MQ=255)                   ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgtgt                                            <  1:296987/68‑1 (MQ=255)                   ctactaCGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgtgt                                            <  1:2476582/68‑1 (MQ=255)                     actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACt                                        >  1:111626/1‑70 (MQ=255)                     actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACt                                        >  1:3150939/1‑70 (MQ=255)                     actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACt                                        >  1:823317/1‑70 (MQ=255)                     actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACt                                        >  1:1826489/1‑70 (MQ=255)                     actacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACt                                        >  1:1232054/1‑70 (MQ=255)                      ctacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTg                                       <  1:749852/70‑1 (MQ=255)                      ctacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGa                                      <  1:2754523/71‑1 (MQ=255)                       tacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgtgt                                            <  1:3936271/64‑1 (MQ=255)                       tacGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgtgt                                            <  1:1772363/64‑1 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAg                                               >  1:1931935/1‑59 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTa                                          >  1:4211513/1‑64 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAgg                                   >  1:1672247/1‑71 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAgg                                   >  1:3581060/1‑71 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAgg                                   >  1:1402534/1‑71 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAgg                                   >  1:458228/1‑71 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAgg                                   >  1:4385433/1‑71 (MQ=255)                         cgcgTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAg                                    >  1:4002729/1‑70 (MQ=255)                             tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAggcggc                               <  1:3806275/71‑1 (MQ=255)                             tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAggcggc                               <  1:1993783/71‑1 (MQ=255)                             tCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAggcggc                               <  1:1254295/71‑1 (MQ=255)                                 gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGgcggcg                              <  1:2925323/68‑1 (MQ=255)                                 gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGgcggcg                              <  1:163006/68‑1 (MQ=255)                                 gTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACAGAAGgcggcg                              <  1:3509650/68‑1 (MQ=255)                                  tGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTg                          <  1:3434784/71‑1 (MQ=255)                                  tGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTg                          <  1:2001236/71‑1 (MQ=255)                                    aTCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGaaa                                                           <  1:766877/37‑1 (MQ=255)                                     tCCGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGc                       <  1:3881409/71‑1 (MQ=255)                                      ccGTCTCGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCa                      >  1:3496813/1‑71 (MQ=255)                                         tctcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:3230686/40‑1 (MQ=39)                                         tctcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTa                                                  <  1:2754083/40‑1 (MQ=39)                                         tctcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGtt                                           <  1:2396549/47‑1 (MQ=255)                                         tctcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAggcggc                               <  1:2302173/59‑1 (MQ=255)                                         tctcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAggcggc                               <  1:2558521/59‑1 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCa                      >  1:4454118/1‑66 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                   >  1:3859082/1‑69 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                   >  1:3819881/1‑69 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                   >  1:3357269/1‑69 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                   >  1:400891/1‑69 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAccc                   >  1:4667616/1‑69 (MQ=255)                                           tcGGCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCAc                     >  1:313191/1‑67 (MQ=255)                                             ggCGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAgg               >  1:1450146/1‑71 (MQ=255)                                               cGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:207174/71‑1 (MQ=255)                                               cGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:1186429/71‑1 (MQ=255)                                               cGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:2875761/71‑1 (MQ=255)                                               cGTGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:194346/71‑1 (MQ=255)                                                 tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTg                          <  1:2634748/56‑1 (MQ=255)                                                 tGGCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTg                          <  1:2640264/56‑1 (MQ=255)                                                  ggCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAgtgt                                            <  1:68495/37‑1 (MQ=38)                                                   gCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATg         >  1:3462769/1‑71 (MQ=255)                                                   gCCTATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATg         >  1:1994515/1‑71 (MQ=255)                                                      tATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGa           >  1:3810459/1‑66 (MQ=255)                                                      tATGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGa           >  1:3862460/1‑66 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:3118277/62‑1 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     <  1:2201052/70‑1 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     >  1:2622744/1‑70 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     >  1:588625/1‑70 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     >  1:4004363/1‑70 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     <  1:3959630/70‑1 (MQ=255)                                                        tGGCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACg     >  1:3941871/1‑70 (MQ=255)                                                          gCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGa           >  1:168766/1‑62 (MQ=255)                                                          gCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGaa   >  1:511053/1‑70 (MQ=255)                                                          gCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAAc  >  1:3226263/1‑71 (MQ=255)                                                          gCTCCGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAAc  >  1:3511296/1‑71 (MQ=255)                                                          gCTACGATCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAAc  >  1:1505102/1‑71 (MQ=255)                                                                aTCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:302810/54‑1 (MQ=255)                                                                aTCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:2817636/54‑1 (MQ=255)                                                                aTCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:2045402/54‑1 (MQ=255)                                                                aTCTGG‑AAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGt             <  1:222411/54‑1 (MQ=255)                                                                      |                                                            ACTGGTCGTTGACTGACACTACTACGCGTCTGGTGATCCGTCTCGGCGTGGCCTATGGCTCCGATCTGGAAAAAGTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAAC  >  W3110S.gb/488737‑488865

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