Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,034,547 T→C 8.8% L40L (TTA→TTG yhhN ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,034,5470TC 8.8% 175.9 / 8.2 68L40L (TTA→TTGyhhNconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base C (2/4);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCAGTGAGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCACGATACG  >  W3110S.gb/4034477‑4034611
                                                                      |                                                                
gcagTGAGGCGCACAGCCCTGCCAGCACGAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGt                                                                  <  1:2552899/71‑1 (MQ=255)
gcagTGAGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGt                                                                  <  1:4426396/71‑1 (MQ=255)
gcagTGAGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGt                                                                  <  1:3360476/71‑1 (MQ=255)
   gTGAGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTaa                                                                >  1:3992887/1‑70 (MQ=255)
      aGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAAcagc                                                            <  1:4319658/71‑1 (MQ=255)
      aGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAAcag                                                             >  1:1801929/1‑70 (MQ=255)
      aGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAAcag                                                             >  1:843451/1‑70 (MQ=255)
             caGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGgag                                                       >  1:486369/1‑69 (MQ=255)
                  cTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAgg                                                <  1:3891402/71‑1 (MQ=255)
                   tGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGga                                                        <  1:884532/62‑1 (MQ=255)
                    gCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGt                                              >  1:4213239/1‑71 (MQ=255)
                        gCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTaa                                                                >  1:407953/1‑49 (MQ=255)
                        gCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAgg                                                         <  1:791933/56‑1 (MQ=255)
                         cACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACg                                         >  1:2017680/1‑71 (MQ=255)
                          aCCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACgg                                        >  1:2385027/1‑71 (MQ=255)
                           ccAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTaa                                                                >  1:1978653/1‑46 (MQ=255)
                           ccAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTaa                                                                >  1:3532208/1‑46 (MQ=255)
                           ccAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAAcagcag                                                          <  1:1867296/52‑1 (MQ=255)
                              gATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTaa                                                                >  1:1444622/1‑43 (MQ=255)
                              gATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:51950/71‑1 (MQ=255)
                               aTAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGttaa                                    >  1:4058906/1‑68 (MQ=255)
                                tAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTaaa                                  >  1:3145579/1‑71 (MQ=255)
                                 aGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTaaa                                  >  1:3379419/1‑70 (MQ=255)
                                 aGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTaaa                                  >  1:1612749/1‑70 (MQ=255)
                                  gCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGt                                              <  1:1637121/57‑1 (MQ=255)
                                  gCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGt                                              <  1:4430408/57‑1 (MQ=255)
                                  gCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAac                                 <  1:3874384/70‑1 (MQ=255)
                                   cTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAacac                               >  1:4050222/1‑71 (MQ=255)
                                   cTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGCAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:2548057/66‑1 (MQ=255)
                                   cTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGCAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:3275297/66‑1 (MQ=255)
                                   cTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGCAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:2862067/66‑1 (MQ=255)
                                   cTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGCAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:2415941/66‑1 (MQ=255)
                                    tAATGGCGGCGAACATCGGCGCTTGCCAGGCCAGTAAca                                                              <  1:3555396/39‑1 (MQ=255)
                                       tGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTa                                    <  1:1955171/62‑1 (MQ=255)
                                         gCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTaa                                   >  1:2232726/1‑61 (MQ=255)
                                         gCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCg                         >  1:651816/1‑71 (MQ=255)
                                            tCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAgaagaa                                                  <  1:1233340/43‑1 (MQ=255)
                                            tCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGttt                                     <  1:2528295/56‑1 (MQ=255)
                                              gAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGcc                    <  1:1263017/71‑1 (MQ=255)
                                              gAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGc                     >  1:1795012/1‑70 (MQ=255)
                                              gAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGc                     >  1:3875835/1‑70 (MQ=255)
                                               aaCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAAc                                          <  1:1750234/48‑1 (MQ=255)
                                               aaCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGttt                                     <  1:699063/53‑1 (MQ=255)
                                               aaCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCa                   <  1:2588567/71‑1 (MQ=255)
                                               aaCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCa                   <  1:3950726/71‑1 (MQ=255)
                                                aCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAg                  >  1:933070/1‑71 (MQ=255)
                                                aCATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAg                  >  1:3413027/1‑71 (MQ=255)
                                                 cATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAgg                 >  1:774135/1‑71 (MQ=255)
                                                 cATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAg                  >  1:779450/1‑70 (MQ=255)
                                                  ttCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAAc                                          <  1:722311/44‑1 (MQ=255)
                                                  aTCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGt                <  1:3090087/71‑1 (MQ=255)
                                                  aTCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGt                <  1:2927800/71‑1 (MQ=255)
                                                     ggCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAgg                                                >  1:513627/1‑36 (MQ=255)
                                                     ggCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAgg                                                >  1:2321896/1‑36 (MQ=255)
                                                     ggCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTgg             <  1:3078759/71‑1 (MQ=255)
                                                       cgcTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGtt                                             >  1:409761/1‑37 (MQ=255)
                                                        gcTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGccc          <  1:18820/71‑1 (MQ=255)
                                                            gccaggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCac        <  1:3465777/69‑1 (MQ=255)
                                                             ccaggccagTAACAGCAGGATAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgat     <  1:4399369/71‑1 (MQ=255)
                                                             ccaggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCGAGCGTTGCCAGGTTGGCCCacgat     <  1:2175101/71‑1 (MQ=255)
                                                              caggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgata    >  1:3587162/1‑71 (MQ=255)
                                                               aggccagTAACAGCAGTAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgatac   <  1:1295389/71‑1 (MQ=255)
                                                               aggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgatac   <  1:4281455/71‑1 (MQ=255)
                                                               aggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgatac   <  1:1790140/71‑1 (MQ=255)
                                                               aggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgatac   <  1:4037899/71‑1 (MQ=255)
                                                               aggccagCAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCg                         >  1:4400182/1‑49 (MQ=255)
                                                               aggccagCAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCg                         >  1:4293322/1‑49 (MQ=255)
                                                                ggccagTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacgatacg  <  1:1690370/71‑1 (MQ=255)
                                                                    agTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCacga      <  1:3696121/63‑1 (MQ=255)
                                                                     gTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGc                          <  1:77079/42‑1 (MQ=255)
                                                                     gTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAg                  <  1:1735316/50‑1 (MQ=255)
                                                                      |                                                                
GCAGTGAGGCGCACAGCCCTGCCAGCACCAGATAGCTAATGGCGTCGAACATCGGCGCTTGCCAGGCCAGTAACAGCAGGAGAAGAAGGGTTAACGGTTTAAACACCCAGCGTTGCCAGGTTGGCCCACGATACG  >  W3110S.gb/4034477‑4034611

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: