Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,179,532 T→C 51.9% *328Q (TAA→CAA)  ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,179,5320TC51.9% 2.5 / 68.5 54*328Q (TAA→CAA) ansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base C (8/20);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACT  >  minE/1179464‑1179598
                                                                    |                                                                  
cGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaa                                                                  <  1:2653063/71‑1 (MQ=255)
 gCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGCGCCaaa                                                                 <  1:12243/71‑1 (MQ=255)
 gCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:4842924/71‑1 (MQ=255)
  cTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaaa                                                                <  1:106570/71‑1 (MQ=255)
     atcAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAAc                                                               <  1:2486513/67‑1 (MQ=255)
     aCCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCt                                                             <  1:3792839/71‑1 (MQ=255)
      ccAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTg                                                            >  1:898835/1‑71 (MQ=255)
      ccAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTg                                                            >  1:2892605/1‑71 (MQ=255)
       cAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCttc                                                           <  1:4496678/71‑3 (MQ=255)
       cAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgc                                                           >  1:3271803/1‑71 (MQ=255)
       cAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgc                                                           >  1:858623/1‑71 (MQ=255)
         aaCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgcgc                                                         <  1:44065/71‑1 (MQ=255)
          aCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGcgcg                                                        >  1:1460595/1‑71 (MQ=255)
             gCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCgcggc                                                      >  1:5178588/1‑70 (MQ=255)
              cATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCgcggcg                                                     <  1:2140221/70‑1 (MQ=255)
              cATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCgcggcg                                                     <  1:4744465/70‑1 (MQ=255)
               aTTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:165190/57‑1 (MQ=255)
               aTTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:3727288/57‑1 (MQ=255)
               aTTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:3893390/57‑1 (MQ=255)
               aTTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:3971997/57‑1 (MQ=255)
                 tACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACt                                                 >  1:4923195/1‑71 (MQ=255)
                  aCCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACt                                                 >  1:746055/1‑70 (MQ=255)
                  aCCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACt                                                 <  1:2948543/70‑1 (MQ=255)
                   ccTGCTTAGCCAGTAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgc                                                           <  1:79608/59‑1 (MQ=255)
                   ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTgc                                                           <  1:4251452/59‑1 (MQ=255)
                   ccTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCTCGGCGAACTg                                                <  1:1472010/70‑1 (MQ=255)
                        ttAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGccg                                          >  1:3939705/1‑69 (MQ=255)
                           gCCAGGAACTGTATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGc                                       <  1:4285207/71‑1 (MQ=255)
                            ccAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:2729597/44‑1 (MQ=255)
                                  aCTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAAcc                                                              >  1:2774454/1‑41 (MQ=255)
                                   cTGGATACTGAAACCATTCGCAAGGCCATGAGCCaaa                                                                 <  1:3934465/37‑1 (MQ=255)
                                    tGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCt                                      <  1:3785432/63‑1 (MQ=255)
                                     ggATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCa                             >  1:919498/1‑71 (MQ=255)
                                     ggATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCa                             >  1:5183983/1‑71 (MQ=255)
                                      gATACTGAAACCATTCGCAAGGCCATGAGCTAAAAc                                                               >  1:1504160/1‑36 (MQ=255)
                                       aTACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGAc                                              <  1:1191379/52‑1 (MQ=255)
                                       aTACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGAc                                              <  1:548966/52‑1 (MQ=255)
                                          cTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACttt                         <  1:3540266/70‑1 (MQ=255)
                                            gAAACCATTTGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCt                                      <  1:3121499/55‑1 (MQ=255)
                                            gAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTaa                       >  1:901442/1‑70 (MQ=255)
                                            gAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAAt                      >  1:2287663/1‑71 (MQ=255)
                                            gAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGc                                            <  1:2351412/49‑1 (MQ=255)
                                             aaaCCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATc                     >  1:4550020/1‑71 (MQ=255)
                                                 cATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGc                 >  1:224134/1‑71 (MQ=255)
                                                 cATTCGCAAGGCCATGAGCTAAAACCTGCGAGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGc                 >  1:2001226/1‑71 (MQ=255)
                                                  aTTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTa                        <  1:2576795/63‑1 (MQ=255)
                                                  aTTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCt                <  1:3881346/71‑1 (MQ=255)
                                                          ggCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCATATCTCACTTTAATCGTGCTCACATTAc        >  1:2865008/1‑71 (MQ=255)
                                                          ggCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACATTa         >  1:1571229/1‑70 (MQ=255)
                                                            ccATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACATTACGt      <  1:216491/71‑1 (MQ=255)
                                                              aTGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca              <  1:1624897/61‑1 (MQ=255)
                                                              aTGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTca              <  1:365010/61‑1 (MQ=255)
                                                               tGAGCTAAAACCTGCGCGGAGAACTGACGCTTTGCTca                                    >  1:5331619/1‑38 (MQ=255)
                                                                 aGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACt  <  1:4188420/70‑1 (MQ=255)
                                                                    |                                                                  
CGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCTTTGCTCACATCTCACTTTAATCGTGCTCACATTACGTGACT  >  minE/1179464‑1179598

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: