Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,852,3240AC40.0% 33.6 / 100.3 70I283S (ATC→AGC) rsgAribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base C (13/15);  total (32/38)
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.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCC  >  minE/2852258‑2852391
                                                                  |                                                                   
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:3420681/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:3645799/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:4082575/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:4601851/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:2140367/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:1904762/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:185387/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:1811397/61‑1 (MQ=37)
          aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATACGACAAATCACTCAGGGCTTTg                                                                 <  1:3391353/61‑1 (MQ=25)
           aTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg                                                                 <  1:4467940/60‑1 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:4371975/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:1790294/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:4498201/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:3579138/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:4780622/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:3184812/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:1173593/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:1989108/1‑60 (MQ=25)
            tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc                                                               >  1:2894957/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:2644703/1‑60 (MQ=25)
                 atatTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGcc                                                <  1:1665753/71‑1 (MQ=255)
                      tACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCagag                                            <  1:1606548/70‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                                     <  1:4195010/71‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                                     <  1:4121507/71‑1 (MQ=255)
                            gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc                                     <  1:4012701/71‑1 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGAt                                                                    >  1:3109736/1‑35 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                                >  1:411856/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                                >  1:414741/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                                >  1:2535318/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                                >  1:3149772/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg                                >  1:1324699/1‑71 (MQ=255)
                                 aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCAGGCTCCAGGTGCCAGAGGCCGAACTCAc                                 >  1:4782047/1‑70 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:2265221/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:510481/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:4613349/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:4538645/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:3373312/71‑1 (MQ=255)
                                         aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:2578753/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:4408771/71‑1 (MQ=255)
                                                      aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac           <  1:2174942/71‑1 (MQ=255)
                                                      aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac           <  1:2471259/71‑1 (MQ=255)
                                                         cccTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg                        <  1:3075735/55‑1 (MQ=255)
                                                               gTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCAccgcc  >  1:316564/1‑71 (MQ=255)
                                                                  |                                                                   
TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCC  >  minE/2852258‑2852391

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: