Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,395,375 C→G 17.4% T27T (ACC→ACG setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,395,3750CG17.4% 158.2 / 25.3 69T27T (ACC→ACGsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base C (35/22);  new base G (3/9);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTC  >  Exported/1395307‑1395442
                                                                    |                                                                   
cccGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                                   <  1:2579060/71‑1 (MQ=255)
 ccGCAGTCTCCAGCGCGAAATCGTATGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                                   >  1:635835/1‑70 (MQ=255)
  cGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                                 >  1:5225024/1‑71 (MQ=255)
   gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                                >  1:660005/1‑71 (MQ=255)
   gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                                >  1:5306107/1‑71 (MQ=255)
   gCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGGTTTTAATCGTTGCCTTTCTCACCGGTAt                                                                >  1:1389234/1‑71 (MQ=255)
        ctcCAGCGCGAAATCGTTTGACCGGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCgg                                                           >  1:3612387/1‑71 (MQ=255)
         tcCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCg                                                                    <  1:968773/61‑1 (MQ=255)
         tcCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCggg                                                          >  1:4405991/1‑71 (MQ=255)
           cAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                                        >  1:4176297/1‑71 (MQ=255)
           cAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                                        >  1:729124/1‑71 (MQ=255)
                cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                                   >  1:1719305/1‑71 (MQ=255)
                cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                                   >  1:7046/1‑71 (MQ=255)
                cgAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                                   >  1:1707370/1‑71 (MQ=255)
                 gAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                                  <  1:4810592/55‑1 (MQ=255)
                 gAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                                  <  1:2350916/55‑1 (MQ=255)
                  aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                                     <  1:1375379/51‑1 (MQ=255)
                  aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                                     <  1:2395392/51‑1 (MQ=255)
                  aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCAcc                                                                     <  1:1220346/51‑1 (MQ=255)
                  aaaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGt                                                                  <  1:3022663/54‑1 (MQ=255)
                   aaTCGTTTGACCTGACCTCGAGGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                                >  1:5423496/1‑55 (MQ=255)
                   aaTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTAt                                                                >  1:1540246/1‑55 (MQ=255)
                    aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                                   <  1:2479554/51‑1 (MQ=255)
                    aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                                   <  1:902423/51‑1 (MQ=255)
                    aTCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCgg                                                                   <  1:1561261/51‑1 (MQ=255)
                       gTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                            >  1:2179918/1‑71 (MQ=255)
                         ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                                 >  1:3335461/1‑48 (MQ=255)
                         ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                                 >  1:459019/1‑48 (MQ=255)
                         ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTa                                                                 >  1:4508323/1‑48 (MQ=255)
                         ttGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                                           >  1:1240193/1‑70 (MQ=255)
                                gacctCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                                        <  1:3104764/50‑1 (MQ=255)
                                    tCGGCGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                                           <  1:5507853/59‑1 (MQ=255)
                                    tCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCaa                                                <  1:3094563/54‑1 (MQ=255)
                                     ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCTCTCTTCAGACTCCTACCCTAAGTa                                <  1:1912440/66‑1 (MQ=255)
                                     ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                <  1:961379/66‑1 (MQ=255)
                                     ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                <  1:5461045/66‑1 (MQ=255)
                                     ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                <  1:5606809/66‑1 (MQ=255)
                                     cGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGgcg                                                        >  1:4075428/1‑45 (MQ=255)
                                     cGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAccc                                            >  1:570409/1‑57 (MQ=255)
                                      gACGGCGTTTATAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTAttt                             >  1:345318/1‑71 (MQ=255)
                                        cGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTAtattc                           <  1:1329297/71‑2 (MQ=255)
                                          gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt                         <  1:5433845/68‑1 (MQ=255)
                                          gcatttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCt                         <  1:433957/68‑1 (MQ=255)
                                            gTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                       >  1:3885576/1‑71 (MQ=255)
                                            atttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCtcg                       >  1:3109744/2‑69 (MQ=255)
                                             tttttACTTGTCGCCTTTCTGACGGGTATTGCGGgcg                                                        <  1:1611621/37‑1 (MQ=12)
                                             tttttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCa                                                 <  1:4844851/44‑1 (MQ=16)
                                             tttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGa                                         <  1:4913978/52‑1 (MQ=255)
                                              ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                     >  1:5049856/1‑71 (MQ=255)
                                              ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                     >  1:2495962/1‑71 (MQ=255)
                                               tttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCtaa                                   >  1:3363584/1‑54 (MQ=255)
                                               tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTg                                                   >  1:4477010/1‑40 (MQ=255)
                                               tttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                     >  1:4926033/1‑70 (MQ=255)
                                                ttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGa                   >  1:4581485/1‑71 (MQ=255)
                                                   aTCGTTGCCTTTCTCACCGGTATTGCGGGTGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                 <  1:2188262/70‑1 (MQ=255)
                                                   aTCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                <  1:4649244/71‑1 (MQ=255)
                                                    tCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                                           <  1:4514188/43‑1 (MQ=255)
                                                    tCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAAcccc                                           <  1:5166761/43‑1 (MQ=255)
                                                    tCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCtt                        <  1:5064535/62‑1 (MQ=255)
                                                       ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa            <  1:5577743/71‑1 (MQ=255)
                                                       ttGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTa            <  1:3752649/71‑1 (MQ=255)
                                                         gCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGt             >  1:5173362/1‑68 (MQ=255)
                                                          ccTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACAt         <  1:3398826/71‑1 (MQ=255)
                                                             ttCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGaa               >  1:1215099/1‑62 (MQ=255)
                                                             ttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc      >  1:4008181/1‑71 (MQ=255)
                                                             ttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc      >  1:796732/1‑71 (MQ=255)
                                                              tctcACCGGTATTGCGGGCGCTCTGCAAACCCCGac                                        >  1:1206445/1‑36 (MQ=255)
                                                              tctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTGTTACCGATGAAGTACATGccc     >  1:2692158/1‑71 (MQ=255)
                                                              tctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGccc     >  1:1987298/1‑71 (MQ=255)
                                                               ctcACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGccc     >  1:557475/1‑70 (MQ=255)
                                                                  aCCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTa                                <  1:2298029/40‑1 (MQ=255)
                                                                  aCCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTc  <  1:4786186/70‑1 (MQ=255)
                                                                    |                                                                   
CCCGCAGTCTCCAGCGCGAAATCGTTTGACCTGACCTCGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTC  >  Exported/1395307‑1395442

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: