Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT21.7% 175.3 / 28.0 60A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/19);  new base T (13/0);  total (41/19)
Fisher's exact test for biased strand distribution p-value = 5.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  REL606/1747828‑1747978
                                                                           |                                                                           
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTGCAgg                                                                             <  2:1401301/76‑1 (MQ=255)
 agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                            <  1:898827/76‑1 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCCGGCa                                                                           >  2:59720/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                           >  1:401562/1‑76 (MQ=255)
   gCTGATAATGTTNTTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                          >  1:12933/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                         >  1:935024/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                        >  1:684629/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                        >  1:62440/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                        >  1:1206555/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                        >  1:772732/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                        >  1:1152683/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                        >  1:540211/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                        >  1:1188461/1‑76 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                        >  2:574189/1‑76 (MQ=255)
     tGATAATGTNATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtg                                                                        >  1:30144/1‑75 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCGGCCATTACACGATGAACGCTTATGAACTCCNGTCAcatc                                                                       >  1:1257390/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                       <  1:824989/76‑1 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                       >  2:1709/1‑76 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                      >  2:1133902/1‑76 (MQ=255)
       aTAATGTTANTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                      >  1:33683/1‑73 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                    >  2:696754/1‑76 (MQ=255)
         aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                    <  2:160105/76‑1 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTGTGAACTCCAGGCACTCCGCCa                                                                   >  1:1338011/1‑76 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                   >  1:524784/1‑76 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                                <  2:109606/76‑1 (MQ=255)
              tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatatc                                                               >  1:814487/1‑75 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                             >  1:611954/1‑76 (MQ=255)
                tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                             <  2:729945/76‑1 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                            >  2:759462/1‑76 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                            >  1:482804/1‑76 (MQ=255)
                 ttatATTTTATTCTTTATTTCATGCAGCCATTNCAAGATGANCGCTTATGAACTCCAGGCANTCCGCCATAttttt                                                            >  2:1182771/1‑76 (MQ=255)
                   atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                          >  1:597203/1‑76 (MQ=255)
                    tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                         >  1:1057691/1‑76 (MQ=255)
                    tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                         >  1:756653/1‑76 (MQ=255)
                      ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                       >  1:215074/1‑76 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                      >  2:444139/1‑76 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                      <  1:1140806/76‑1 (MQ=255)
                           tttttttttTCATGCAGCCATTACACGATGAACGCTTATGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                  <  1:1233444/73‑1 (MQ=255)
                           ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                  <  1:778684/76‑1 (MQ=255)
                           ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                  <  1:605664/76‑1 (MQ=255)
                            tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                                 <  1:361712/76‑1 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATTAACTGCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                <  2:1225321/76‑1 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               >  1:1076483/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               >  2:1268266/1‑76 (MQ=255)
                               ttATTTCANGAAGCCATTANACGATGAACGCTTATGAACTCCNGGCACTCCGCCATATTTTTGCCATGACTATtga                                              <  2:1188461/76‑1 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                            >  1:578435/1‑76 (MQ=255)
                                  tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                           >  1:1161812/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                          >  1:301048/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                          >  1:331639/1‑76 (MQ=255)
                                    tCATGCAGCCATTACACGATGAACGCTTATGCACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                         <  2:227906/76‑1 (MQ=255)
                                      aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                       <  1:383198/76‑1 (MQ=255)
                                        gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                     <  1:215187/76‑1 (MQ=255)
                                           gCCATTACACGATGAACGCTTATGACCTCCGGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAc                                  <  1:557915/76‑1 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                <  1:584075/76‑1 (MQ=255)
                                              aTTACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                               <  1:519313/76‑1 (MQ=255)
                                              aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                               >  1:674870/1‑76 (MQ=255)
                                                tACACGATGAACNCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                             >  1:15258/1‑76 (MQ=255)
                                                  cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                           <  2:674287/76‑1 (MQ=255)
                                                   acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                          >  1:1082819/1‑76 (MQ=255)
                                                           cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                  <  1:828371/76‑1 (MQ=255)
                                                            gCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 <  2:773466/76‑1 (MQ=255)
                                                             tttATGAACTCCAGCCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg                <  1:1016272/75‑1 (MQ=255)
                                                                aTGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGa             <  1:444332/76‑1 (MQ=255)
                                                                aTGAACTCCAGGCACTCCGCCATATTTNTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGa             >  1:1418766/1‑76 (MQ=255)
                                                                     ccccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg        <  1:943689/74‑1 (MQ=255)
                                                                      cccAGGCACTCCGCCATATTTTTGACATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa       <  1:631918/75‑1 (MQ=255)
                                                                       ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa      >  2:463616/1‑76 (MQ=255)
                                                                           gCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  <  1:1139185/76‑1 (MQ=255)
                                                                           gCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  <  1:649886/76‑1 (MQ=255)
                                                                           |                                                                           
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  REL606/1747828‑1747978

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