Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,805,151 C→A 9.4% I162I (ATC→ATA proW → glycine betaine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,805,1510CA 9.4% 162.7 / 2.7 53I162I (ATC→ATAproWglycine betaine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (22/26);  new base A (2/3);  total (24/29)
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

GGTGCCTGGTCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGTCCACTGCTTGA  >  W3110S.gb/2805083‑2805219
                                                                    |                                                                    
ggTGCCTGGTCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCgg                                                                    <  1:2346802/71‑1 (MQ=255)
     cTGGTCACAGGCAATGGTTACCCTGGCGCTGGTGTTAACAGCCCTGCTGTTCTGTATCGTCATAGGTTTGc                                                               >  1:3097453/1‑71 (MQ=255)
     cTGGTCACAGGCAATGGTTACCCTGGCGCTGGTGTTAACAGCCCTGCTGTTCTGTATCGTCATAGGTTTGc                                                               >  1:593629/1‑71 (MQ=255)
      tGGTCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGcc                                                              >  1:2676185/1‑71 (MQ=255)
        ggcGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGcc                                                              <  1:2610273/67‑1 (MQ=255)
        gTCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            >  1:1224767/1‑71 (MQ=255)
         tCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGtt                                                           <  1:1315910/71‑1 (MQ=255)
          cGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTg                                                          >  1:2378821/1‑71 (MQ=255)
              ggCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGa                                                      <  1:2746974/71‑1 (MQ=255)
              ggCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGa                                                      <  1:1974463/71‑1 (MQ=255)
               gCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGat                                                     <  1:3756413/71‑1 (MQ=255)
                cAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCg                                                                     >  1:3073316/1‑54 (MQ=255)
                cAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCg                                                                     >  1:2538021/1‑54 (MQ=255)
                cAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCg                                                                     >  1:2326736/1‑54 (MQ=255)
                cAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCg                                                                     >  1:905703/1‑54 (MQ=255)
                    ggTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATAtggc                                                >  1:729603/1‑71 (MQ=255)
                    ggTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATAtggc                                                >  1:159436/1‑71 (MQ=255)
                    ggTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATAtggc                                                >  1:764415/1‑71 (MQ=255)
                     gTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGata                                                    >  1:1957882/1‑66 (MQ=255)
                     gTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATAtggct                                               >  1:1391371/1‑71 (MQ=255)
                      tGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCgg                                                                    <  1:942682/49‑1 (MQ=255)
                      tGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCgg                                                                    <  1:583832/49‑1 (MQ=255)
                              gcgcTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCgag                                         <  1:207300/68‑1 (MQ=255)
                              gcgcTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCgag                                         <  1:3754564/68‑1 (MQ=255)
                              gcgcTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCgag                                         <  1:3210225/68‑1 (MQ=255)
                                       ttAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            <  1:411308/40‑1 (MQ=255)
                                       ttAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            <  1:2592299/40‑1 (MQ=255)
                                       ttAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            <  1:1862511/40‑1 (MQ=255)
                                       ttAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            <  1:3564719/40‑1 (MQ=255)
                                       ttAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGt                                                            <  1:1350696/40‑1 (MQ=255)
                                       ttAACCGCCATGCTGTTCTGTATCGTCATCGGTTTGCCGTTg                                                          >  1:144446/1‑42 (MQ=255)
                                       ttAACAGCCCTGCTGTTCTGTATCGTCATAGGTTTGCCGTTGGGGATCTGGCTGGCgaga                                        <  1:3263000/60‑1 (MQ=255)
                                       ttAACAGCCCTGCTGTTCTGTATCGTCATAGGTTTGCCGTTGGGGATCTGGCTGGCgaga                                        <  1:1173027/60‑1 (MQ=255)
                                       ttAACAGCCCTGCTGTTCTGTATCGTCATAGGTTTGCCGTTGGGGATCTGGCTGGCgaga                                        <  1:1476289/60‑1 (MQ=255)
                                        tAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAgcg                            >  1:431410/1‑71 (MQ=255)
                                         aaCCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTg                                                          >  1:1244979/1‑40 (MQ=255)
                                         aaCCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAgcg                            >  1:450542/1‑70 (MQ=255)
                                           ccGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGa                                                      <  1:3833671/42‑1 (MQ=255)
                                           ccGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGa                                                      <  1:2066487/42‑1 (MQ=255)
                                           ccGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGAGAgcggc                          >  1:1881843/1‑70 (MQ=255)
                                             gCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGaa                       <  1:2835435/71‑1 (MQ=255)
                                                ctgctgTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAat                    <  1:1304887/71‑1 (MQ=255)
                                                 tgctgTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAat                    >  1:2392689/1‑70 (MQ=255)
                                                   ctgTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAatt                   <  1:2852590/69‑1 (MQ=255)
                                                       tCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGt             >  1:2905928/1‑71 (MQ=255)
                                                        cTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGt             <  1:852759/70‑1 (MQ=255)
                                                          gTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCgag                                         <  1:563466/40‑1 (MQ=255)
                                                           tATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCgag                                         <  1:3071063/39‑1 (MQ=255)
                                                               gTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGTCCa          >  1:365149/1‑66 (MQ=255)
                                                               gTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGTCCACTGCt     >  1:3843858/1‑71 (MQ=255)
                                                                  aTCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGTCCACTGCTTGa  >  1:1728294/1‑71 (MQ=255)
                                                                   tCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAat                    <  1:3564120/52‑1 (MQ=255)
                                                                   tCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAat                    <  1:1029968/52‑1 (MQ=255)
                                                                    cGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGt             >  1:3757678/1‑58 (MQ=255)
                                                                    |                                                                    
GGTGCCTGGTCGCAGGCAATGGTGACTCTGGCGCTGGTGTTAACCGCCCTGCTGTTCTGTATCGTCATCGGTTTGCCGTTGGGGATATGGCTGGCGAGAAGTCCGCGAGCGGCGAAAATTATTCGTCCACTGCTTGA  >  W3110S.gb/2805083‑2805219

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: