Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 318,754 C→T 6.4% *176* (TAG→TAA priC ← primosomal replication protein N''

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE318,7540CT 6.4% 252.6 / 5.5 78*176* (TAG→TAApriCprimosomal replication protein N''
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base T (2/3);  total (39/39)
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.39e-01

CACAATCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGC  >  minE/318687‑318822
                                                                   |                                                                    
cacaATCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAg                                                                   >  1:1003719/1‑71 (MQ=255)
   aaTCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAg                                                                   <  1:3445008/68‑1 (MQ=255)
   aaTCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCgg                                                                >  1:2326364/1‑71 (MQ=255)
     tCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCggg                                                               <  1:2564352/70‑1 (MQ=255)
     tCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCggg                                                               <  1:2463052/70‑1 (MQ=255)
     tCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCg                                                                 >  1:1600033/1‑68 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             >  1:900728/1‑70 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             <  1:3182961/70‑1 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             >  1:741730/1‑70 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             <  1:352840/70‑1 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             >  1:1846686/1‑70 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGtt                                                             >  1:55408/1‑70 (MQ=255)
       aaaTCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGt                                                              >  1:49311/1‑69 (MQ=255)
          tCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAc                                                         >  1:2207046/1‑71 (MQ=255)
            gACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcgc                                                       >  1:2247664/1‑71 (MQ=255)
            gACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcgc                                                       >  1:712357/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:2024418/71‑1 (MQ=255)
                     tttGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTg                                              <  1:1116303/71‑1 (MQ=255)
                     tttGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTg                                              <  1:3119934/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)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGc                                                                  >  1:2996163/1‑46 (MQ=255)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGc                                                                  >  1:3442789/1‑46 (MQ=255)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGc                                                                  >  1:2245314/1‑46 (MQ=255)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAAcg                                                        <  1:2377913/56‑1 (MQ=255)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                         >  1:1722771/1‑71 (MQ=255)
                          cATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                         <  1:1379801/71‑1 (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:3280941/1‑51 (MQ=255)
                              gTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGTTAGCGCGTTTAAc                                                         >  1:925117/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:62452/42‑1 (MQ=255)
                                    ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTa                                                            <  1:2523465/42‑1 (MQ=255)
                                    ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                         <  1:2771113/61‑1 (MQ=255)
                                    ggCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTc                                         <  1:1683747/61‑1 (MQ=255)
                                     gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                >  1:1551120/1‑69 (MQ=255)
                                     gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                >  1:2181712/1‑69 (MQ=255)
                                     gCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                >  1:3424086/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:2952291/70‑1 (MQ=255)
                                                       ttCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGt             <  1:3376003/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:3642031/42‑1 (MQ=255)
                                                                ttGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTc                                <  1:1204668/42‑1 (MQ=255)
                                                                 tGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGgcgcgc  >  1:2667413/1‑71 (MQ=255)
                                                                  gCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGgcgcgc  <  1:1225058/70‑1 (MQ=255)
                                                                   |                                                                    
CACAATCAAATCGACAGCCAGTTTGACATCGTCAGGGGCATTTTCCAGTGACATATTCTCTCCATTGCTAGCGGGTTAAACGCGCTAACCTGTTTTCGATTTTTTCCAGCGCATGGCGGCAGCGTGCCAGGCGCGC  >  minE/318687‑318822

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: