Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,964,464 C→A 10.2% A390S (GCG→TCG)  yhjE ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,964,4640CA10.2% 207.1 / 5.8 68A390S (GCG→TCG) yhjEpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (5/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

TAGCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTAACAGCGCACC  >  W3110S.gb/3964399‑3964533
                                                                 |                                                                     
tAGCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGt                                                                  >  1:2786503/1‑71 (MQ=255)
tAGCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGt                                                                  >  1:1415037/1‑71 (MQ=255)
tAGCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGt                                                                  >  1:820871/1‑71 (MQ=255)
  gCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTg                                                                <  1:2870093/71‑1 (MQ=255)
    ccGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGt                                                               <  1:3675503/70‑1 (MQ=255)
    ccGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGt                                                               <  1:771753/70‑1 (MQ=255)
    ccGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTa                                                              <  1:1351461/71‑1 (MQ=255)
    ccGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTa                                                              <  1:1778715/71‑1 (MQ=255)
     cGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTAc                                                                    >  1:3068125/1‑64 (MQ=255)
     cGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTAc                                                                    >  1:2377944/1‑64 (MQ=255)
     cGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTAc                                                                    >  1:1858783/1‑64 (MQ=255)
      gTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAag                                                            <  1:455043/71‑1 (MQ=255)
      gTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAag                                                            <  1:728304/71‑1 (MQ=255)
         gttagtCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAagag                                                          <  1:443561/70‑1 (MQ=255)
            agtCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAAc                                                       <  1:2832931/70‑1 (MQ=255)
             gtCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGa                                                     <  1:2536097/71‑1 (MQ=255)
              tCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAAc                                                       <  1:441752/68‑1 (MQ=255)
                tGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATg                                                   >  1:2644843/1‑70 (MQ=255)
                       aGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtg                                           >  1:1167585/1‑71 (MQ=255)
                       aGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtg                                           <  1:1933053/71‑1 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACg                                                                   >  1:2327652/1‑46 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtgt                                          <  1:1001968/71‑1 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtgt                                          <  1:3345575/71‑1 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGATa                                                                     >  1:2156593/1‑44 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGATa                                                                     >  1:559536/1‑44 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGATa                                                                     >  1:2351962/1‑44 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGATa                                                                     >  1:2891431/1‑44 (MQ=255)
                        ggCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGATa                                                                     >  1:118385/1‑44 (MQ=255)
                         gCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAAc                                                       <  1:2220243/57‑1 (MQ=255)
                         gCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAAc                                                       <  1:2998177/57‑1 (MQ=255)
                               atatatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCAc                                   <  1:3524665/71‑1 (MQ=255)
                                 atatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGt                                                                  >  1:3777582/1‑38 (MQ=255)
                                 atatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACtt                                 >  1:1104499/1‑71 (MQ=255)
                                 atatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACtt                                 >  1:3817218/1‑71 (MQ=255)
                                 atatGGCGCAACGGAAGCCCCGAGAATCGACGATACGTTGTAAGAGAACGATGCTCCGGTGTAACGCAc                                   <  1:901911/69‑1 (MQ=255)
                                  tatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGaa                                                        >  1:818405/1‑47 (MQ=255)
                                  tatGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGa                                                     >  1:2216129/1‑50 (MQ=255)
                                    tGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTg                              >  1:1982037/1‑71 (MQ=255)
                                    tGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTg                              >  1:3136942/1‑71 (MQ=255)
                                     ggCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTaa                                        >  1:1636130/1‑60 (MQ=255)
                                     ggCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTaa                                        >  1:836623/1‑60 (MQ=255)
                                      gcgcAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGt                             >  1:2010048/1‑70 (MQ=255)
                                      gcgcAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGt                             >  1:1350417/1‑70 (MQ=255)
                                       cgcAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGa                                                     >  1:1576999/1‑45 (MQ=255)
                                       cgcAACGGAAGCCCCGAGAATCGACGATACGTTGTAAGAGAACGATGCTCCGGTGTAACGCAc                                   <  1:2828158/63‑1 (MQ=255)
                                          aaCGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtgt                                          <  1:3812266/53‑1 (MQ=255)
                                          aaCGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCgg                          >  1:3814536/1‑69 (MQ=255)
                                          aaCGGAAGCCCCGAGAATCGACGCTACGTTGAAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCg                           >  1:910719/1‑68 (MQ=255)
                                            cGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGaaa                       >  1:1244673/1‑70 (MQ=255)
                                              gAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACACACTTCTGTCGGAAACAg                    >  1:2487341/1‑71 (MQ=255)
                                                aGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGct                  <  1:923764/71‑1 (MQ=255)
                                                aGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGct                  <  1:2109918/71‑1 (MQ=255)
                                                    ccGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGaaa                       <  1:3767811/62‑1 (MQ=255)
                                                    ccGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGaaa                       <  1:1426061/62‑1 (MQ=255)
                                                    ccGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGaaa                       <  1:2865477/62‑1 (MQ=255)
                                                      gagaATCGACGCTACGTTGTAAGAGAACGATGCTCCGgtgt                                          <  1:658052/41‑1 (MQ=255)
                                                      gagaATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGt             >  1:2322558/1‑70 (MQ=255)
                                                      gagaATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGATCTGGt             >  1:3695155/1‑70 (MQ=255)
                                                      gagaATCGACGCTACATTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTa            <  1:1975772/71‑1 (MQ=255)
                                                      gagaATAGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGt             >  1:494326/1‑70 (MQ=255)
                                                        gaATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTaa           >  1:11842/1‑70 (MQ=255)
                                                        gaATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTAAc          <  1:3107722/71‑1 (MQ=255)
                                                          aTCGACGCTACGTTGTAAGAGAACGATGCTCCGgtgt                                          <  1:1527601/37‑1 (MQ=255)
                                                           tCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAAc                      <  1:360371/56‑1 (MQ=255)
                                                           tCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTAACAgc       <  1:1748842/71‑1 (MQ=255)
                                                           tCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTCTCGGAAAc                      <  1:1715349/56‑1 (MQ=255)
                                                             gACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAg                    >  1:1016850/1‑56 (MQ=255)
                                                                gCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTAACAGCGCAcc  <  1:338649/71‑1 (MQ=255)
                                                                 cTACGTTGTAAGAGAACGATGCTCCGGTGTAACGc                                     <  1:2502996/35‑1 (MQ=255)
                                                                 cTACGTTGTAAGAGAACGATGCTCCGGTGTAACGc                                     <  1:1540603/35‑1 (MQ=255)
                                                                 cTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGt             <  1:1979584/59‑1 (MQ=255)
                                                                 cTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGt             <  1:1358907/59‑1 (MQ=255)
                                                                 |                                                                     
TAGCCCGTAGTTAGTCTGCAACCAGGCTGCGATATATGGCGCAACGGAAGCCCCGAGAATCGACGCTACGTTGTAAGAGAACGATGCTCCGGTGTAACGCACTTCTGTCGGAAACAGCTCTGGTAACAGCGCACC  >  W3110S.gb/3964399‑3964533

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: ‑