Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,666,801 T→A 12.3% intergenic (‑45/‑117) guaB ← / → xseA IMP dehydrogenase/exonuclease VII, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,666,8010TA12.3% 152.3 / 13.0 57intergenic (‑45/‑117)guaB/xseAIMP dehydrogenase/exonuclease VII, large subunit
Reads supporting (aligned to +/- strand):  ref base T (19/31);  new base A (2/5);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAGGG  >  minE/1666732‑1666865
                                                                     |                                                                
gagCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCAtt                                                                 <  1:184881/71‑1 (MQ=255)
 agCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATta                                                                <  1:1638960/71‑1 (MQ=255)
 agCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGAGGCATta                                                                <  1:2934472/71‑1 (MQ=255)
   cttcttTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATtat                                                               <  1:2037513/70‑1 (MQ=255)
          aGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGc                                                        <  1:164628/70‑1 (MQ=255)
            cGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTa                                                     <  1:807919/71‑1 (MQ=255)
              aTACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAAc                                                   >  1:614868/1‑71 (MQ=255)
              aTACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAAc                                                   >  1:682474/1‑71 (MQ=255)
                     gCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTg                                             >  1:831446/1‑70 (MQ=255)
                     gCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTg                                             >  1:2043424/1‑70 (MQ=255)
                     gCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGc                                            <  1:3657760/71‑1 (MQ=255)
                      cATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAAcc                                                  <  1:1654117/64‑1 (MQ=255)
                      cATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAAcc                                                  <  1:3658454/64‑1 (MQ=255)
                      cATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCa                                           <  1:1532550/71‑1 (MQ=255)
                       aTGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCg                                                 <  1:2399114/64‑1 (MQ=255)
                       aTGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCa                                           <  1:1331108/70‑1 (MQ=255)
                        ttggCAATATCTCGACCAGAGTTGTTAATAAATATTGCCGCGGTATTATACAGAGCGTAACCg                                                 <  1:1582001/61‑1 (MQ=255)
                        tGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGa                                                >  1:452251/1‑64 (MQ=255)
                        tGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATc                                         <  1:2799161/71‑1 (MQ=255)
                        tGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATc                                         >  1:1964503/1‑71 (MQ=255)
                        tGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATc                                         >  1:1968858/1‑71 (MQ=255)
                        tGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATc                                         <  1:850906/71‑1 (MQ=255)
                                atCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCtt                                 >  1:1487736/1‑71 (MQ=255)
                                atCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCtt                                 >  1:1440857/1‑71 (MQ=255)
                                 tctcGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCttt                                >  1:1289851/1‑71 (MQ=255)
                                     gACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCtttt                               <  1:2408549/68‑1 (MQ=255)
                                     gACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCttt                                <  1:2131018/67‑1 (MQ=255)
                                     gACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCttt                                <  1:795423/67‑1 (MQ=255)
                                      aCCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGc                            >  1:2488014/1‑70 (MQ=255)
                                       ccAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCa                           >  1:1623031/1‑70 (MQ=255)
                                          gactGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCaaaaa                       <  1:393692/68‑1 (MQ=255)
                                             tGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAt                     <  1:2165728/70‑1 (MQ=255)
                                               gTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtgc                   <  1:2391967/70‑1 (MQ=255)
                                               gTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtgc                   <  1:805972/70‑1 (MQ=255)
                                                 taataaATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtgcttg                >  1:2397126/1‑71 (MQ=255)
                                                  aataaataTTGCCGCGGCATTATACATAGCGTAACCGAtt                                              <  1:1553172/40‑1 (MQ=255)
                                                    taaataTTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTa             <  1:266661/71‑1 (MQ=255)
                                                    taaataTTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTa             >  1:485256/1‑71 (MQ=255)
                                                    taaataTTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTa             >  1:3728598/1‑71 (MQ=255)
                                                     aaataTTGCCGCGGCATTATACAGAGCTTAACCGATTGc                                            <  1:1770645/39‑1 (MQ=255)
                                                     aaataTTGCCGCGGCATTATACAGAGCGTAACCGATTGc                                            <  1:3277194/39‑1 (MQ=255)
                                                       atatTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATcc          >  1:3151975/1‑71 (MQ=255)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATc                                         <  1:1854972/40‑1 (MQ=37)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATc                                         <  1:1349962/40‑1 (MQ=37)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATc                                         <  1:3319228/40‑1 (MQ=37)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATc                                         <  1:699419/40‑1 (MQ=37)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATc                                         <  1:2830050/40‑1 (MQ=37)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATCTACCCCTTTTTg                             >  1:239947/1‑52 (MQ=38)
                                                       ataaTGCCGCGGCAATATACAGAGCGTAACCGATTGCATCTACCCCTTTTTg                             >  1:588050/1‑52 (MQ=38)
                                                          ttGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCg       <  1:1767571/71‑1 (MQ=255)
                                                            gCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCg       >  1:2075715/1‑69 (MQ=255)
                                                             ccGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCttt                                <  1:3013768/43‑1 (MQ=255)
                                                             ccGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtg                    <  1:975878/55‑1 (MQ=255)
                                                              cgcgGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtgcttgc               <  1:2001755/59‑1 (MQ=255)
                                                              cgcgGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtg                    <  1:1022790/54‑1 (MQ=255)
                                                                cgGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAAtgc                   <  1:2291744/53‑1 (MQ=255)
                                                                cgGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAggg  >  1:621022/1‑70 (MQ=255)
                                                                cgGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAggg  >  1:2476061/1‑70 (MQ=255)
                                                                     ttATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAg    >  1:3452985/1‑63 (MQ=255)
                                                                     ttATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAg    >  1:3046235/1‑63 (MQ=255)
                                                                     ttATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAg    >  1:2976103/1‑63 (MQ=255)
                                                                     ttATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAg    >  1:2692666/1‑63 (MQ=255)
                                                                     |                                                                
GAGCTTCTTTAGCGATACGTAGCATGGGCAATATCTCGACCAGAGTGGTTAATAAATATTGCCGCGGCATTATACAGAGCGTAACCGATTGCATCTACCCCTTTTTGCAAAAAATGCTTGCTATCCCCGAAGGG  >  minE/1666732‑1666865

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: