Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,852,3140GT41.2% 89.2 / 68.4 68G286G (GGC→GGArsgAribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base T (13/15);  total (31/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 = 6.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCAC  >  minE/2852258‑2852382
                                                        |                                                                    
tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg                                                        >  1:4886460/1‑71 (MQ=255)
tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg                                                        >  1:1473326/1‑71 (MQ=255)
tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg                                                        >  1:1971755/1‑71 (MQ=255)
tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg                                                        >  1:1813510/1‑71 (MQ=255)
  gtgtTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGtt                                                      >  1:1250706/1‑71 (MQ=255)
     tttGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCg                                                   >  1:1174973/1‑71 (MQ=255)
     tttGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCg                                                   >  1:4412269/1‑71 (MQ=255)
       tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 >  1:845769/1‑71 (MQ=255)
       tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 >  1:4317865/1‑71 (MQ=255)
       tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 >  1:2349188/1‑71 (MQ=255)
       tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 >  1:2195896/1‑71 (MQ=255)
        gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 <  1:3059882/70‑1 (MQ=255)
        gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 <  1:3732542/70‑1 (MQ=255)
        gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 <  1:2325207/70‑1 (MQ=255)
        gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc                                                 <  1:1600312/70‑1 (MQ=255)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:4293240/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:3048949/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:296000/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:2948145/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:2855677/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:1811397/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:3420681/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:3645799/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:4082575/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:2140367/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:4601851/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:1904762/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:185387/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATACGACAAATCACTCAGGGCTTTg                                                        <  1:3391353/61‑1 (MQ=25)
           aTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                        <  1:4467940/60‑1 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:3184812/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:4371975/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:3579138/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:4498201/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:4780622/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:1173593/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:2894957/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:1790294/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:2762747/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:2761603/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:2742086/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:1989108/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                      >  1:2644703/1‑60 (MQ=25)
                 atatTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGcc                                       <  1:1665753/71‑1 (MQ=255)
                      tACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCagag                                   <  1:1606548/70‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                            <  1:4012701/71‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                            <  1:4195010/71‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                            <  1:4121507/71‑1 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGAt                                                           >  1:3109736/1‑35 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                       >  1:414741/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                       >  1:3149772/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                       >  1:1324699/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                       >  1:2535318/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                       >  1:411856/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCAGGCTCCAGGTGCCAGAGGCCGAACTCAc                        >  1:4782047/1‑70 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:2265221/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:2578753/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:4538645/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:4613349/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:3373312/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:510481/71‑1 (MQ=255)
                                          tGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg               <  1:1514603/70‑1 (MQ=255)
                                            gAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTgggg            <  1:72508/71‑1 (MQ=255)
                                               ttCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGcc                                       <  1:12341/41‑1 (MQ=255)
                                                   aCAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAat     <  1:2741247/71‑1 (MQ=255)
                                                      aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac  <  1:2471259/71‑1 (MQ=255)
                                                      aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac  <  1:4408771/71‑1 (MQ=255)
                                                      aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac  <  1:2174942/71‑1 (MQ=255)
                                                        |                                                                    
TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCAC  >  minE/2852258‑2852382

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: