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

AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747829‑1747975
                                                                          |                                                                        
agGCTGATAATGTTCTTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                         >  2:344312/1‑76 (MQ=255)
 ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa                                                                        <  1:317821/76‑1 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:1035520/1‑76 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCAc                                                                       <  1:796516/76‑1 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                      >  1:760981/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:725594/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                     >  1:299940/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                     >  1:842334/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                     >  1:1264311/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                     >  1:659417/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                     >  1:504116/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  2:919809/1‑76 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     <  1:438912/76‑1 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                    >  1:161262/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAccga                                                                    >  1:1120483/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:1490345/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:433247/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:1632992/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:900981/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                    >  1:251156/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                    <  1:518919/76‑1 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                    >  1:1495453/1‑76 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                   >  1:829273/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac                                                                   >  1:1001938/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcca                                                                   >  1:178933/1‑75 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  2:1704931/76‑1 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   <  2:1068938/76‑1 (MQ=255)
       tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                  >  1:466851/1‑76 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                 >  2:345477/1‑76 (MQ=255)
        aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGcc                                                                 <  1:301830/76‑1 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                <  1:699121/76‑1 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                               <  1:1767674/76‑1 (MQ=255)
           tttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                              <  1:1187263/75‑1 (MQ=255)
            ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             >  1:676887/1‑76 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            >  1:948120/1‑76 (MQ=255)
              atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                           >  2:661699/1‑76 (MQ=255)
               tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          >  2:696076/1‑76 (MQ=255)
                 tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                        >  1:188718/1‑76 (MQ=255)
                  atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                       >  1:1450609/1‑76 (MQ=255)
                   tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      >  2:407477/1‑76 (MQ=255)
                   tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      <  1:189720/76‑1 (MQ=255)
                    atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                     >  2:192698/1‑76 (MQ=255)
                    atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                     >  1:706037/1‑76 (MQ=255)
                     ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                    >  2:1073480/1‑76 (MQ=255)
                      tttattttttatttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                   <  1:1405512/76‑1 (MQ=255)
                      tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                   <  1:1555410/76‑1 (MQ=255)
                         aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAca                                                <  2:146757/76‑2 (MQ=255)
                          ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                               <  2:1505553/76‑1 (MQ=255)
                          ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                               <  1:69508/76‑1 (MQ=255)
                           tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                              <  2:1685060/76‑1 (MQ=255)
                           tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                              >  1:917999/1‑76 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                            <  2:1706420/76‑1 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                            <  1:1631415/76‑1 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  2:50411/1‑76 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           <  2:1271057/76‑1 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  1:91827/1‑76 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  1:1168593/1‑76 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  1:1551411/1‑76 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         >  2:388546/1‑76 (MQ=255)
                                 tttCATGCAGCCATTACACTATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                        <  2:1349878/76‑1 (MQ=255)
                                   tCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                      >  1:330971/1‑76 (MQ=255)
                                     aTGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                    <  1:176449/76‑1 (MQ=255)
                                       gCAACCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                  >  1:1180223/1‑76 (MQ=255)
                                        cAGCCATTACACGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                 <  1:1572785/76‑1 (MQ=255)
                                        cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                 >  2:1641027/1‑76 (MQ=255)
                                           ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              >  2:1721821/1‑76 (MQ=255)
                                             aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                            >  2:1071631/1‑76 (MQ=255)
                                               tACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                          >  1:259972/1‑76 (MQ=255)
                                                acacGATGAACGCTTATGACCTCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                         <  1:1035072/76‑1 (MQ=255)
                                                acacGATGAACGCTTAGGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                         <  1:1035081/76‑1 (MQ=255)
                                                 cacGATGACCGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                        <  1:644980/76‑1 (MQ=255)
                                                   cGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGc                      >  1:1738977/1‑76 (MQ=255)
                                                     aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                    <  2:339417/76‑1 (MQ=255)
                                                     aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                    >  2:389605/1‑76 (MQ=255)
                                                      tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                   <  2:70763/76‑1 (MQ=255)
                                                      tGAACGCTTATGAACTCCAGGCACGCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                   <  2:1469871/76‑1 (MQ=255)
                                                          cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc               <  1:781374/76‑1 (MQ=255)
                                                                 gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa        <  1:476496/76‑1 (MQ=255)
                                                                  accTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg       <  1:1456693/74‑1 (MQ=255)
                                                                  aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg       >  2:954944/1‑76 (MQ=255)
                                                                      ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   >  2:1436996/1‑76 (MQ=255)
                                                                      ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:380644/76‑1 (MQ=255)
                                                                      ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:1446205/76‑1 (MQ=255)
                                                                      ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   <  1:1024568/76‑1 (MQ=255)
                                                                       cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa  <  1:1694341/76‑1 (MQ=255)
                                                                          |                                                                        
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747829‑1747975

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