Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 318,760 G→C 7.1% T174T (ACC→ACG priC ← primosomal replication protein N''

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE318,7600GC 7.1% 215.2 / 5.7 70T174T (ACC→ACGpriCprimosomal replication protein N''
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base C (2/3);  total (33/37)
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.91e-01

AATCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGCTTCATAGG  >  minE/318690‑318830
                                                                      |                                                                      
aaTCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCgg                                                                        >  1:2326364/1‑71 (MQ=255)
  tCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCggg                                                                       <  1:2564352/70‑1 (MQ=255)
  tCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCggg                                                                       <  1:2463052/70‑1 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     <  1:352840/70‑1 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     <  1:3182961/70‑1 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     >  1:1846686/1‑70 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     >  1:55408/1‑70 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     >  1:741730/1‑70 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                                     >  1:900728/1‑70 (MQ=255)
    aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGt                                                                      >  1:49311/1‑69 (MQ=255)
       tCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAc                                                                 >  1:2207046/1‑71 (MQ=255)
         gACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcgc                                                               >  1:712357/1‑71 (MQ=255)
         gACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcgc                                                               >  1:2247664/1‑71 (MQ=255)
              ccAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTaa                                                          >  1:161489/1‑71 (MQ=255)
               cAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGTTTAACCCGCTag                                                          <  1:2131568/70‑2 (MQ=38)
                aGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAAcc                                                        >  1:3291133/1‑71 (MQ=255)
                aGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAAcc                                                        >  1:3314044/1‑71 (MQ=255)
                  tttGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTg                                                      <  1:3119934/71‑1 (MQ=255)
                  tttGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTg                                                      <  1:1116303/71‑1 (MQ=255)
                  tttGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTg                                                      <  1:2024418/71‑1 (MQ=255)
                    tGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGtt                                                    >  1:1967934/1‑71 (MQ=255)
                     gACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTa                                                           <  1:555045/63‑1 (MQ=255)
                     gACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTa                                                           <  1:2225937/63‑1 (MQ=255)
                      aCATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcg                                                                <  1:1544715/57‑1 (MQ=255)
                       cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcg                                                                <  1:2377913/56‑1 (MQ=255)
                       cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                                 <  1:1379801/71‑1 (MQ=255)
                       cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                                 >  1:1722771/1‑71 (MQ=255)
                        aTCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCg                                                <  1:3530937/71‑1 (MQ=255)
                        aTCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCg                                                <  1:1079709/71‑1 (MQ=255)
                         tCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGttta                                                                   <  1:3405436/51‑3 (MQ=255)
                          cGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGttta                                                                   <  1:491615/50‑3 (MQ=255)
                           gTCAGGGGGATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTaaa                                                                  >  1:1040370/1‑50 (MQ=255)
                           gTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGTTTAAc                                                                 >  1:925117/1‑51 (MQ=255)
                           gTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGTTTAAc                                                                 >  1:3280941/1‑51 (MQ=255)
                           gTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTaaa                                                                  >  1:2484134/1‑50 (MQ=255)
                           gTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTaa                                                                   >  1:1567417/1‑49 (MQ=255)
                            tCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGAttt                                            >  1:591725/1‑71 (MQ=255)
                                gggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        <  1:3264525/71‑1 (MQ=255)
                                 ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTa                                                                    <  1:2523465/42‑1 (MQ=255)
                                 ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTa                                                                    <  1:62452/42‑1 (MQ=255)
                                 ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                                 <  1:1683747/61‑1 (MQ=255)
                                 ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                                 <  1:2771113/61‑1 (MQ=255)
                                  gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        >  1:3424086/1‑69 (MQ=255)
                                  gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        >  1:2181712/1‑69 (MQ=255)
                                  gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        >  1:1551120/1‑69 (MQ=255)
                                   cATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGa                                               <  1:1931461/61‑1 (MQ=255)
                                    aTTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAgc                                    >  1:2765640/1‑71 (MQ=255)
                                     ttttCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        <  1:842226/66‑1 (MQ=255)
                                     ttttCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAgcg                                   <  1:1422568/71‑1 (MQ=255)
                                          cAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTTTTTTCg                                                <  1:480709/53‑1 (MQ=255)
                                          cAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCg                                                <  1:1457404/53‑1 (MQ=255)
                                          cAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAAGCGCTAACCTGTTTTCg                                                <  1:3487688/53‑1 (MQ=255)
                                           aGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggc                             <  1:3516549/71‑1 (MQ=255)
                                            gTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggcg                            >  1:2890849/1‑71 (MQ=255)
                                            gTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggcg                            <  1:312023/71‑1 (MQ=255)
                                             tGATATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggcgg                           >  1:2489485/1‑71 (MQ=255)
                                             tGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggcgg                           >  1:2456053/1‑71 (MQ=255)
                                             tGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATggcg                            >  1:1421924/1‑70 (MQ=255)
                                                 atatTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGc                       <  1:3174203/71‑1 (MQ=255)
                                                    ttCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGt                     <  1:3376003/70‑1 (MQ=255)
                                                    ttCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGt                     <  1:2952291/70‑1 (MQ=255)
                                                    ttCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTg                    <  1:2285052/71‑1 (MQ=255)
                                                     tctctcCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTcc                                       >  1:3551980/1‑51 (MQ=255)
                                                        ctcCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAg                >  1:517386/1‑71 (MQ=255)
                                                         tcCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCa                 <  1:2338513/69‑1 (MQ=255)
                                                          ccATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGgc              >  1:1563195/1‑71 (MQ=255)
                                                           cATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAg                                     >  1:3249197/1‑47 (MQ=255)
                                                           cATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAg                                     >  1:612055/1‑47 (MQ=255)
                                                             ttGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        <  1:1204668/42‑1 (MQ=255)
                                                             ttGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                        <  1:3642031/42‑1 (MQ=255)
                                                              tGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGgcgcgc          >  1:2667413/1‑71 (MQ=255)
                                                               gCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGgcgcgc          <  1:1225058/70‑1 (MQ=255)
                                                                   gCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGCTTc       <  1:2950062/69‑1 (MQ=255)
                                                                     gggTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGCTTCATAg   >  1:3041466/1‑71 (MQ=255)
                                                                      ggTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGCTTCATAgg  >  1:450522/1‑71 (MQ=255)
                                                                      |                                                                      
AATCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGCTTCATAGG  >  minE/318690‑318830

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: