Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT22.1% 186.6 / 33.7 68A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/18);  new base T (15/0);  total (50/18)
Fisher's exact test for biased strand distribution p-value = 7.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747828‑1747976
                                                                           |                                                                         
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                           <  2:76714/76‑1 (MQ=255)
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                           >  2:422254/1‑76 (MQ=255)
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                           <  2:290148/76‑1 (MQ=255)
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                           >  2:256582/1‑76 (MQ=255)
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                           >  2:256571/1‑76 (MQ=255)
 agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                          >  2:908491/1‑76 (MQ=255)
 agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                          >  1:951713/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa                                                                         <  2:860837/76‑1 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:962228/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:1534917/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:85381/1‑76 (MQ=255)
    ntGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                       >  1:119230/2‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:588691/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:160889/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                      >  1:1482337/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                      >  1:864455/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                      >  1:490216/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                      >  1:1079776/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccag                                                                     >  1:1407029/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc                                                                     >  1:89493/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                     >  1:336416/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                     >  1:1546876/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                     >  2:139635/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                     <  2:468180/76‑1 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                    >  1:206910/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                    >  1:136133/1‑76 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                    <  2:1473202/76‑1 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                    >  2:215075/1‑76 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                    <  2:252084/76‑1 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   >  1:1222289/1‑76 (MQ=255)
        tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   >  2:823735/1‑76 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                  >  1:1012383/1‑73 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                  >  2:1075605/1‑76 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                  >  2:629096/1‑76 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                 <  2:773114/76‑1 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                 >  2:1445352/1‑76 (MQ=255)
            gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                               >  1:975219/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                              >  2:713831/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                            >  2:898625/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                            >  2:313577/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                            >  2:401057/1‑76 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCCCTCCGCCATAtttt                                                           <  2:1245039/76‑1 (MQ=255)
                  tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                         >  2:933315/1‑76 (MQ=255)
                  tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                         >  1:910735/1‑76 (MQ=255)
                   atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                        >  2:808799/1‑76 (MQ=255)
                   atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                        <  2:1392776/76‑1 (MQ=255)
                    tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                       >  1:1192969/1‑76 (MQ=255)
                     tttttattttttattTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                      <  1:664839/75‑1 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                  >  2:788258/1‑76 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                  >  2:1285229/1‑76 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                  >  1:1211262/1‑76 (MQ=255)
                          aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                 >  2:1552438/1‑76 (MQ=255)
                            ttttttttTCATGCAGCCATTACACGATGAACGCTTATGAACCCCAGGCCCTCCGCCATATTTTTGCCATGACTAt                                               <  1:761092/74‑1 (MQ=255)
                            tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                               >  1:839154/1‑76 (MQ=255)
                             gTTTATTTCATGCNGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                              >  1:1105772/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                             >  2:1090134/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                          <  2:289425/76‑1 (MQ=255)
                                  tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                         >  2:12843/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                        >  1:813538/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                        <  2:1036185/76‑1 (MQ=255)
                                    tCATGCAGCCATTACAAGATGACCGTTTATAACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                       <  1:741330/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGACCGTTTATGAACTCCGGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      <  1:83002/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      <  1:421296/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTAGGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      <  1:1485445/76‑1 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                     <  1:951405/76‑1 (MQ=255)
                                          aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                 <  2:643952/76‑1 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                              >  2:870788/1‑76 (MQ=255)
                                              aTTACAGGACGAACGCTTAATATCTCCGGCCATTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                             <  1:1498547/76‑1 (MQ=255)
                                              aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                             >  2:672577/1‑76 (MQ=255)
                                                 acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                          >  1:39322/1‑76 (MQ=255)
                                                   acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                        <  1:1481109/76‑1 (MQ=255)
                                                   acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                        >  1:934055/1‑76 (MQ=255)
                                                    cGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGc                       <  2:1332067/76‑1 (MQ=255)
                                                      agaaaCGCTTATAACCTCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                     <  1:802109/73‑1 (MQ=255)
                                                         aaCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                  <  2:1288236/76‑1 (MQ=255)
                                                           cGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                <  1:1089911/76‑1 (MQ=255)
                                                           cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                >  1:842872/1‑76 (MQ=255)
                                                             cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg              <  2:1136921/76‑1 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg      <  2:504924/76‑1 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg      >  1:350678/1‑76 (MQ=255)
                                                                     cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg      <  1:333166/76‑1 (MQ=255)
                                                                         aGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg  <  1:1144001/76‑1 (MQ=255)
                                                                         aGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg  <  2:183811/76‑1 (MQ=255)
                                                                           |                                                                         
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747828‑1747976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: