Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT30.9% 184.4 / 75.8 84A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/20);  new base T (26/0);  total (64/20)
Fisher's exact test for biased strand distribution p-value = 2.12e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747829‑1747976
                                                                          |                                                                         
agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                          <  1:312134/76‑1 (MQ=255)
agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                          >  1:963037/1‑76 (MQ=255)
 ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCTAGTCa                                                                         >  1:35054/1‑76 (MQ=255)
 ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                         >  1:753454/1‑76 (MQ=255)
 ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa                                                                         <  2:1513900/76‑1 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:926397/1‑76 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:820599/1‑76 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                        >  1:69077/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                       >  1:623164/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:1987666/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:1625751/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:158444/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                       >  1:1353439/1‑76 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgc                                                                      >  1:163507/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcga                                                                      >  1:1362166/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                      >  1:1334364/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccg                                                                      >  1:1825833/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                      >  1:342235/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                      >  1:706267/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                      >  1:1988481/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                      >  1:2023557/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                      >  1:1449166/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                      >  2:839886/1‑76 (MQ=255)
    ngATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                      >  1:1355500/2‑75 (MQ=255)
    cgATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                      >  1:1901844/2‑75 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcaca                                                                     >  1:1987994/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                     >  1:415044/1‑74 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                    >  1:6093/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac                                                                    >  1:1032144/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtagc                                                                    >  1:279954/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                    >  1:1224875/1‑76 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   >  1:1830381/1‑76 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   >  1:1709138/1‑76 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                   <  1:1873577/76‑1 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                  <  1:1428771/76‑1 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                  >  2:608382/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                 >  1:441887/1‑76 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                                >  2:1968537/1‑76 (MQ=255)
            ntatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                              >  1:31084/2‑76 (MQ=255)
              atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCATAttg                                                            >  1:1660089/1‑75 (MQ=255)
              atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                            <  2:1262759/76‑1 (MQ=255)
               tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                           >  2:1687342/1‑76 (MQ=255)
                 tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                         >  2:1825066/1‑76 (MQ=255)
                  atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                        >  1:227675/1‑76 (MQ=255)
                          ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                >  1:583256/1‑76 (MQ=255)
                           tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                               <  1:1775992/76‑1 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                             <  1:491224/76‑1 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                            >  1:32739/1‑76 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAATATtga                                            >  2:469112/1‑76 (MQ=255)
                               tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                           >  2:772143/1‑76 (MQ=255)
                               tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                           <  2:194641/76‑1 (MQ=255)
                               tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCAATTTTTGCCATGAATAatgatg                                          >  1:978479/1‑76 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                          <  2:1258166/76‑1 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                          >  1:1865060/1‑76 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                          >  2:749760/1‑76 (MQ=255)
                                 tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                         >  1:1055117/1‑76 (MQ=255)
                                 tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                         >  2:1342498/1‑76 (MQ=255)
                                  ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                        <  2:1610832/76‑1 (MQ=255)
                                   tCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                       >  1:356568/1‑76 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      >  1:458585/1‑76 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      >  2:1774574/1‑76 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      <  1:1666133/76‑1 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      >  2:1545047/1‑76 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                      >  1:1812589/1‑76 (MQ=255)
                                       gCAGCCATTACACGATGAACGCTTATGAACTCCAGGTACTCTGCCATATTTTTGCCATGACTATTGATGAATgcgc                                   >  2:1463043/1‑76 (MQ=255)
                                         aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                 >  2:246432/1‑76 (MQ=255)
                                           ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGaac                               <  1:1411697/76‑3 (MQ=255)
                                           ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                               <  1:1910429/76‑1 (MQ=255)
                                            cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                              >  1:8401/1‑76 (MQ=255)
                                             aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                             <  1:1263651/76‑1 (MQ=255)
                                             aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                             >  1:229940/1‑76 (MQ=255)
                                                  acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                        <  2:9976/76‑1 (MQ=255)
                                                  acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                        >  2:265163/1‑76 (MQ=255)
                                                      tGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                    <  1:289047/76‑1 (MQ=255)
                                                      tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                    >  1:1338807/1‑76 (MQ=255)
                                                       gagcGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                   <  1:1671558/73‑1 (MQ=255)
                                                       gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                   <  1:622242/76‑1 (MQ=255)
                                                        aaCGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                  <  2:935242/76‑1 (MQ=255)
                                                            cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg              >  2:2052801/1‑76 (MQ=255)
                                                            cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg              >  2:1135775/1‑76 (MQ=255)
                                                            cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATCCGCGACCTGGATTGCCCAGACCgg              >  2:1047015/1‑76 (MQ=255)
                                                             ttATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt             >  2:1965438/1‑76 (MQ=255)
                                                              tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg            >  1:1024019/1‑76 (MQ=255)
                                                                 gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa         >  1:1797653/1‑76 (MQ=255)
                                                                  aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg        >  1:1804013/1‑76 (MQ=255)
                                                                   aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt       <  1:436855/76‑1 (MQ=255)
                                                                   aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt       <  1:1524789/76‑1 (MQ=255)
                                                                     tCCAGGTACTCCGCCATATATTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa     <  1:1463043/76‑1 (MQ=255)
                                                                     tCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa     <  2:2012292/76‑1 (MQ=255)
                                                                      ccAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa    <  1:1919009/76‑1 (MQ=255)
                                                                        aGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg  <  2:877351/76‑1 (MQ=255)
                                                                          |                                                                         
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAG  >  REL606/1747829‑1747976

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