Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,395,405 A→C 19.7% T37T (ACA→ACC setB → lactose/glucose efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,395,4050AC19.7% 105.7 / 33.1 66T37T (ACA→ACCsetBlactose/glucose efflux system
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base C (3/10);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCG  >  Exported/1395344‑1395475
                                                             |                                                                      
ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCTCTCTTCAGACTCCTACCCTAAGTa                                                                 <  1:1912440/66‑1 (MQ=255)
ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                                                 <  1:961379/66‑1 (MQ=255)
ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                                                 <  1:5606809/66‑1 (MQ=255)
ccgcGGCATTTTTACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCTAAGTa                                                                 <  1:5461045/66‑1 (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)
         ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                                                      >  1:5049856/1‑71 (MQ=255)
         ttttAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTAcc                                                      >  1:2495962/1‑71 (MQ=255)
          tttACTTGTCGCCTTTCTGACGGGTATTGCGGGCGCTCTTCAGACTCCTACCCtaa                                                                    >  1:3363584/1‑54 (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)
               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:796732/1‑71 (MQ=255)
                        ttCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc                                       >  1:4008181/1‑71 (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)
                                  tATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                                                 <  1:2040793/51‑1 (MQ=255)
                                  tATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGa                                                 <  1:5633364/51‑1 (MQ=255)
                                   aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                               <  1:5156775/68‑4 (MQ=255)
                                   aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                               <  1:4691839/68‑4 (MQ=255)
                                   aTTGCGGGCGCTCTTCAGACTCCTACCCTAAGTATATTCCTCGCAGATGAACTGAAAGCCCGTcctat                               <  1:1824549/68‑4 (MQ=255)
                                    ttGCGGGCGCTCTGCAAACCCCGACACTCAGTAttt                                                              <  1:1616264/36‑1 (MQ=255)
                                    ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAg                                               <  1:4564215/51‑1 (MQ=255)
                                    ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtg                           <  1:4153850/71‑1 (MQ=255)
                                    ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtg                           <  1:3216434/71‑1 (MQ=255)
                                    ttGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACGATGAAGTAc                                            >  1:261351/1‑53 (MQ=255)
                                      gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                        >  1:4262473/1‑40 (MQ=255)
                                      gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTa                                                        >  1:3901986/1‑40 (MQ=255)
                                      gCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggt                         >  1:4023571/1‑71 (MQ=255)
                                       cGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggtg                        >  1:4463386/1‑71 (MQ=255)
                                        gggCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtggtgg                       <  1:2389281/71‑1 (MQ=255)
                                           cgcTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAt                    >  1:5515954/1‑71 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                  <  1:2538815/37‑1 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                  <  1:3362466/37‑1 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATg                                                  <  1:4276862/37‑1 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGcc                                       >  1:4824744/1‑48 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg                          <  1:5537199/61‑1 (MQ=255)
                                               ctGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGAtgg                          <  1:2533110/61‑1 (MQ=255)
                                                tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc               >  1:459481/1‑71 (MQ=255)
                                                tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc               >  1:2153550/1‑71 (MQ=255)
                                                tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttcttc               >  1:5036779/1‑71 (MQ=255)
                                                tGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttct                 <  1:5132943/69‑1 (MQ=255)
                                                 gCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAtt                   >  1:2390881/1‑66 (MQ=255)
                                                  cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACa                                           >  1:2365508/1‑41 (MQ=255)
                                                  cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACa                                           >  1:2636017/1‑41 (MQ=255)
                                                  cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttctt                <  1:3241071/68‑1 (MQ=255)
                                                  cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTcttctt             >  1:696417/1‑71 (MQ=255)
                                                  cAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTcttctt             >  1:2719088/1‑71 (MQ=255)
                                                   aaaCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATtcttct              <  1:1268916/69‑1 (MQ=255)
                                                     aCCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGAttc                  >  1:4026626/1‑63 (MQ=255)
                                                       cccGACACTCAGTATTTTTCTTACCGATGAAGTACa                                           >  1:1671558/1‑36 (MQ=255)
                                                        ccGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCgg       <  1:1553053/71‑1 (MQ=255)
                                                             acTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGa                             >  1:2832268/1‑44 (MQ=255)
                                                             acTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAgcg  <  1:4657633/71‑1 (MQ=255)
                                                             |                                                                      
CGACGGCGTTTTTAATCGTTGCCTTTCTCACCGGTATTGCGGGCGCTCTGCAAACCCCGACACTCAGTATTTTTCTTACCGATGAAGTACATGCCCGTCCGGCGATGGTGGGATTCTTCTTTACCGGCAGCG  >  Exported/1395344‑1395475

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