Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,315,024 C→T 10.3% D213N (GAC→AAC)  alsA ← fused D‑allose transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,315,0240CT10.3% 205.8 / 11.3 68D213N (GAC→AAC) alsAfused D‑allose transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (5/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGA  >  W3110S.gb/4314956‑4315094
                                                                    |                                                                      
gTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCg                                                                       >  1:569776/1‑70 (MQ=255)
gTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCg                                                                       >  1:430579/1‑70 (MQ=255)
 tCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCg                                                                       >  1:4467620/1‑69 (MQ=255)
 tCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGc                                                                      >  1:1768145/1‑70 (MQ=255)
  cATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCa                                                                     >  1:210407/1‑70 (MQ=255)
   aTTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCaaa                                                                   >  1:822380/1‑71 (MQ=255)
   aTTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCaaa                                                                   <  1:3167123/71‑1 (MQ=255)
    tttGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACAGTATAGCGGTCGCAAAt                                                                  <  1:4100456/71‑1 (MQ=255)
     ttGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATa                                                                 >  1:4358420/1‑71 (MQ=255)
      tGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATa                                                                 >  1:834740/1‑70 (MQ=255)
      tGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATAc                                                                >  1:3746840/1‑71 (MQ=255)
       gACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATAc                                                                <  1:1556811/70‑1 (MQ=255)
          acaTCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGc                                                             >  1:194216/1‑70 (MQ=255)
          acaTCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCg                                                            <  1:4494577/71‑1 (MQ=255)
          acaTCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCg                                                            <  1:266859/71‑1 (MQ=255)
            aTCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACg                                                               <  1:2010450/66‑1 (MQ=255)
              cGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCAAATACGGCGaa                                                          >  1:4496880/1‑69 (MQ=255)
              cGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGaa                                                          >  1:3087837/1‑69 (MQ=255)
               gCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCAAATa                                                                 >  1:3338527/1‑61 (MQ=255)
               gCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCAAATACGGCGAAtt                                                        <  1:196706/70‑1 (MQ=255)
               gCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGTCGAAttt                                                       <  1:1002750/71‑1 (MQ=255)
               gCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAAttt                                                       <  1:2478883/71‑1 (MQ=255)
                cTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTc                                                      >  1:2600195/1‑71 (MQ=255)
                 ttACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGaa                                                          <  1:2666667/66‑1 (MQ=255)
                  tACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTc                                                      >  1:1910089/1‑69 (MQ=255)
                   aCTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGc                                                   >  1:2988395/1‑71 (MQ=255)
                     tATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATAc                                                                <  1:3108172/56‑1 (MQ=255)
                       tGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAAc                                               >  1:3676559/1‑71 (MQ=255)
                         ccGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACt                                              >  1:504469/1‑70 (MQ=255)
                         ccGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACt                                              >  1:3375991/1‑70 (MQ=255)
                          cGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTa                                            <  1:2628851/71‑1 (MQ=255)
                           gCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTAt                                           <  1:219076/71‑1 (MQ=255)
                           gCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTAt                                           <  1:1671910/71‑1 (MQ=255)
                            cTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCa                                                                     >  1:1141771/1‑44 (MQ=255)
                              gCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCAAATACGGCGAATTTCCGCCAACTTATGc                                         >  1:4021119/1‑70 (MQ=255)
                                aaaCGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                         >  1:3784089/1‑68 (MQ=255)
                                  aCGCTGCTGCCGTCTTTCATCACCGTATAGCGGTTGCAAATACGGCGAATTTCCGCCAACTTATGCGAGat                                    <  1:1448146/71‑1 (MQ=255)
                                   cGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGata                                   >  1:1299157/1‑71 (MQ=255)
                                     ctgctgCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGAtata                                 <  1:3051968/71‑1 (MQ=255)
                                       gctgcCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCgag                                      >  1:3862612/1‑64 (MQ=255)
                                       gctgcCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGa                               >  1:2601224/1‑71 (MQ=255)
                                        ctgcCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                              >  1:3954023/1‑71 (MQ=255)
                                        ctgcCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGAc                              <  1:3761130/71‑1 (MQ=255)
                                           ccGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGata                                   <  1:884723/63‑1 (MQ=255)
                                           ccGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGAt                           <  1:3195030/71‑1 (MQ=255)
                                           ccGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGAt                           >  1:683646/1‑71 (MQ=255)
                                           ccGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGAt                           <  1:3466485/71‑1 (MQ=255)
                                            cGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGAt                           <  1:1385600/70‑1 (MQ=255)
                                             gTCTTTCATCACCGTATAGCGGTTGCAAATACGGCGAATTTCCGCCAACTTATGc                                         >  1:2722837/1‑55 (MQ=255)
                                                tttCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATg                                          >  1:100746/1‑51 (MQ=255)
                                                  tcatcaCCGTATAGCGGTCGCAAATACGGCGAATTTcc                                                     >  1:2726068/1‑38 (MQ=255)
                                                  tcatcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCg                      <  1:3311451/69‑1 (MQ=255)
                                                  tcatcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGt                     <  1:1801731/70‑1 (MQ=255)
                                                  tcatcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTa                    >  1:3599708/1‑71 (MQ=255)
                                                  tcatcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTa                    >  1:333835/1‑71 (MQ=255)
                                                   catcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTAt                                           <  1:5781/47‑1 (MQ=255)
                                                   catcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                         <  1:1571642/49‑1 (MQ=255)
                                                   catcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTa                    <  1:2876631/70‑1 (MQ=255)
                                                    atcaCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTa                    >  1:2224682/1‑69 (MQ=255)
                                                      caCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGc                                         <  1:2578226/46‑1 (MQ=255)
                                                       aCCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCctc               <  1:3877420/71‑1 (MQ=255)
                                                        ccGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCctct              >  1:2585588/1‑71 (MQ=255)
                                                          gTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGAt                           >  1:4123922/1‑56 (MQ=255)
                                                             taGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACg         >  1:24093/1‑71 (MQ=255)
                                                               gCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGc        <  1:211013/70‑1 (MQ=255)
                                                                  gTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACt    <  1:1693072/71‑1 (MQ=255)
                                                                  gTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACt    <  1:3511720/71‑1 (MQ=255)
                                                                   tCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGa                            <  1:2311711/46‑1 (MQ=255)
                                                                    cGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGa  >  1:893790/1‑71 (MQ=255)
                                                                    |                                                                      
GTCATTTGACACATCGCTTACTATGCCGCTGCAAACGCTGCTGCCGTCTTTCATCACCGTATAGCGGTCGCAAATACGGCGAATTTCCGCCAACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGA  >  W3110S.gb/4314956‑4315094

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: