Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,179,007 A→T 13.2% N153Y (AAC→TAC)  ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,179,0070AT13.2% 236.8 / 27.6 91N153Y (AAC→TAC) ansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base A (46/33);  new base T (8/4);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGCC  >  minE/1178938‑1179069
                                                                     |                                                              
aTTAATCTGCTGAATGCGTTGTAGGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTAATAAc                                                              >  1:1513353/1‑71 (MQ=255)
aTTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATaa                                                               >  1:795734/1‑71 (MQ=255)
 ttAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:3915788/69‑1 (MQ=255)
 ttAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:3183194/69‑1 (MQ=255)
 ttAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:2769790/69‑1 (MQ=255)
 ttAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:4116194/69‑1 (MQ=255)
 ttAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:3599800/69‑1 (MQ=255)
  tAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             >  1:2297348/1‑71 (MQ=255)
  tAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             >  1:2956763/1‑71 (MQ=255)
   aaTCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             >  1:1229299/1‑70 (MQ=255)
   aaTCTGATGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             >  1:1525473/1‑70 (MQ=255)
    ttCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCg                                                            <  1:2802561/69‑1 (MQ=255)
    aTCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             <  1:2364612/69‑1 (MQ=255)
    aTCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             <  1:1892859/69‑1 (MQ=255)
     tCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             <  1:874624/68‑1 (MQ=255)
       tgctgAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTg                                                        >  1:2286952/1‑71 (MQ=255)
         ctgAATGCGT‑‑TACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     <  1:2047058/70‑1 (MQ=255)
            aaTGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     <  1:1389161/69‑1 (MQ=255)
            aaTGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     <  1:876272/69‑1 (MQ=255)
            aaTGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     <  1:377216/69‑1 (MQ=255)
            aaTGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATc                                                    <  1:2747075/70‑1 (MQ=255)
            aaTGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATc                                                    <  1:3629484/70‑1 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 >  1:3089551/1‑71 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 <  1:2047543/71‑1 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 >  1:3763190/1‑71 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 >  1:3632822/1‑71 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 <  1:2470867/71‑1 (MQ=255)
              tGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                                                 <  1:1345119/71‑1 (MQ=255)
                 gTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                                                <  1:3199290/53‑1 (MQ=255)
                 gTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:3964609/1‑70 (MQ=255)
                 gTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:1320892/1‑70 (MQ=255)
                 gTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:3356041/1‑70 (MQ=255)
                 gTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:2685330/1‑70 (MQ=255)
                    gTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGa                                                           >  1:168655/1‑55 (MQ=255)
                    gTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGa                                                           >  1:317933/1‑55 (MQ=255)
                     tACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                                          >  1:3453219/1‑71 (MQ=255)
                     tACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                                          <  1:3396130/71‑1 (MQ=255)
                     tACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                                          <  1:3865602/71‑1 (MQ=255)
                     tACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                                          >  1:4086511/1‑71 (MQ=255)
                     tACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg                                          >  1:814275/1‑71 (MQ=255)
                       cGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGc                                         <  1:3697090/70‑1 (MQ=255)
                       cGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCa                                        >  1:2842337/1‑71 (MQ=255)
                         ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     >  1:1234820/1‑56 (MQ=255)
                         ttGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTAt                                                     >  1:3419944/1‑56 (MQ=255)
                              gcgAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaa                                 <  1:840735/71‑1 (MQ=255)
                                  aTTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGcc                             >  1:2405343/1‑71 (MQ=255)
                                   ttATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                                             <  1:878824/38‑1 (MQ=255)
                                     aTCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCaa                                             >  1:581527/1‑52 (MQ=255)
                                      tCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATg                         >  1:3862971/1‑71 (MQ=255)
                                      tCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATg                         >  1:106255/1‑71 (MQ=255)
                                        cGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGcc                       <  1:856121/71‑1 (MQ=255)
                                        cGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGcc                       <  1:1008792/71‑1 (MQ=255)
                                            cAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:741652/1‑43 (MQ=255)
                                            cAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGc                                               >  1:4186684/1‑43 (MQ=255)
                                            cAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATg                   >  1:279638/1‑71 (MQ=255)
                                                gAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaaa                                <  1:995845/54‑1 (MQ=255)
                                                gAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaaa                                <  1:3154332/54‑1 (MQ=255)
                                                gAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaaa                                <  1:2125464/54‑1 (MQ=255)
                                                 aaGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCTCACTACCAAAGCCCATGCCGATGGtttt              <  1:1486129/71‑1 (MQ=255)
                                                 aaGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGtttt              <  1:3233328/71‑1 (MQ=255)
                                                  aGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCa                                        >  1:3455228/1‑44 (MQ=255)
                                                  aGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGtttt              >  1:3296177/1‑70 (MQ=255)
                                                     aaCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt           >  1:138097/1‑70 (MQ=255)
                                                     aaCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt           >  1:2898045/1‑70 (MQ=255)
                                                     aaCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGAt           >  1:3894634/1‑70 (MQ=255)
                                                     aaCGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATg          >  1:3862813/1‑71 (MQ=255)
                                                       cGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCg        <  1:375722/71‑1 (MQ=255)
                                                        gCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGt       <  1:1763059/71‑1 (MQ=255)
                                                        gCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGt       <  1:1716479/71‑1 (MQ=255)
                                                        gCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGt       <  1:1774647/71‑1 (MQ=255)
                                                        gCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGt       <  1:3961064/71‑1 (MQ=255)
                                                         ttcttTTTCAATTACCGgtgtgtcatctgtattcagta                                       >  1:1709466‑M2/2‑17 (MQ=255)
                                                         ttcttTTTCAATTACCGgtgtgtcatctgtattcagta                                       >  1:2928702‑M2/2‑17 (MQ=255)
                                                          tctTTTTCAATTACCGgtgtgtcatctgtattcagta                                       <  1:378878‑M2/37‑22 (MQ=255)
                                                          tctTTTTCAATTACCGgtgtgtcatctgtattcagta                                       >  1:2473875‑M2/1‑16 (MQ=255)
                                                          tctTTTTCAATTACCGgtgtgtcatctgtattcagta                                       <  1:2163955‑M2/37‑22 (MQ=255)
                                                          tctTTTTCAATTACCGgtgtgtcatctgtattcagta                                       <  1:1590322‑M2/37‑22 (MQ=255)
                                                          tctTTTTCAATTACCGgtgtgtcatctgtattcagta                                       <  1:2920206‑M2/37‑22 (MQ=255)
                                                          tctTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGttt     >  1:869037/1‑71 (MQ=255)
                                                           ctTTTTCAATTACCGgtgtgtcatctgtattcagtac                                      >  1:1445550‑M2/1‑15 (MQ=255)
                                                           ctTTTTCAATTACCGgtgtgtcatctgtattcagtac                                      >  1:1231206‑M2/1‑15 (MQ=255)
                                                           ctTTTTCAATTACCGgtgtgtcatctgtattcagtac                                      >  1:3946086‑M2/1‑15 (MQ=255)
                                                           ctTTTTCAATTACCGgtgtgtcatctgtattcagtac                                      >  1:2887637‑M2/1‑15 (MQ=255)
                                                           ctTTTTCAATTACCGgtgtgtcatctgtattcagtac                                      >  1:3024537‑M2/1‑15 (MQ=255)
                                                           ctTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGtttt              >  1:3026120/1‑61 (MQ=255)
                                                             ttttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGcc  >  1:4252897/1‑71 (MQ=255)
                                                             ttttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGcc  >  1:1894822/1‑71 (MQ=255)
                                                             ttttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGcc  >  1:908728/1‑71 (MQ=255)
                                                              tttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGa            >  1:1528700/1‑60 (MQ=255)
                                                              tttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGa            >  1:3648987/1‑60 (MQ=255)
                                                              tttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGc   <  1:1738383/69‑1 (MQ=255)
                                                              tttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGc   <  1:1132191/69‑1 (MQ=255)
                                                               ttCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGc   <  1:401419/68‑1 (MQ=255)
                                                                tCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCa                           >  1:661089/1‑43 (MQ=255)
                                                                tCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCa                           >  1:1475440/1‑43 (MQ=255)
                                                                tCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCa                           >  1:1466874/1‑43 (MQ=255)
                                                                tCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCa                           >  1:292125/1‑43 (MQ=255)
                                                                 cAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCg        >  1:2678475/1‑61 (MQ=255)
                                                                     |                                                              
ATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGCCCATGCCGATGGTTTTGATGCGTTTGCC  >  minE/1178938‑1179069

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: