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

AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA  >  REL606/1747829‑1747974
                                                                          |                                                                       
agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGc                                                                        >  2:896619/1‑76 (MQ=255)
 ggNTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAATTCCATNCa                                                                       <  2:1092123/76‑1 (MQ=255)
 ggCTGATAATGTTATTTACATNTTNTTGTTTATTTCATGCAGCCATNACACGATNANCGCTTATGAACTCCAGGCa                                                                       >  2:193210/1‑76 (MQ=255)
 ggCTGATAAGGATATTTATATTTTATTGTTTATATCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                       >  1:78820/1‑76 (MQ=255)
  gCTGATAATGTTATTTATCTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:1265318/1‑76 (MQ=255)
  gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:370245/1‑76 (MQ=255)
  gCTGATAATGATATTTACATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:109094/1‑76 (MQ=255)
   cTGATAATNTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACTCTTATGAACTCCAGTCACt                                                                     >  1:1292402/1‑76 (MQ=255)
   cTGATAATGTTATTTATNTTTTATTGTTTATGTCATGCAGCCATTACACGATGAACGCTTATGCACACCAGTCACt                                                                     >  1:170359/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                     >  1:373596/1‑75 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:573539/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:1073649/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:679858/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:693768/1‑76 (MQ=255)
   cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:723522/1‑76 (MQ=255)
   cTGATAATGCTCTTTATATCTTATTGTTTATNTGATTCCCCCATTACACGATGAACGCTTATGCNCTCCAGTCACt                                                                     >  1:163720/1‑76 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCNTTACACGATGAACGCTTANGANCTCCAGTCAcca                                                                    >  1:1086810/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCCTTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                    >  1:1107777/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                    >  1:384901/1‑74 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                    >  1:678011/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                    >  1:617756/1‑75 (MQ=255)
    tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                    >  1:499056/1‑75 (MQ=255)
    tGATAATGTTAGTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACta                                                                    >  1:39955/1‑75 (MQ=255)
     gatgatATTATGTATATTTTATTGTTTATTTCACGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                   >  1:1224214/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACGCGATGAACGCTTATGAACTCcactcactta                                                                   >  1:27454/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc                                                                   >  1:771687/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc                                                                   >  1:149945/1‑73 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                   >  1:945876/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                   >  1:393529/1‑74 (MQ=255)
     gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                   >  1:696734/1‑74 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                  >  1:1054650/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                  >  1:221728/1‑73 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                  >  2:671798/1‑76 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                  >  1:145982/1‑76 (MQ=255)
      aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                  >  2:1273507/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                               >  1:966479/1‑76 (MQ=255)
         aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                               >  1:966480/1‑76 (MQ=255)
          tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                              >  1:958463/1‑76 (MQ=255)
          tGTTATTTATATGCTGTGACTGATTCCATGCAGCCATTACGCGCTGAACGCTTATGAACTCCAGACACTTAGCCat                                                              >  1:15213/1‑76 (MQ=255)
           gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCcta                                                             >  2:586608/1‑76 (MQ=255)
           gTTATTTATATTTTATTGTTTATTTCATCCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                             >  2:1204149/1‑76 (MQ=255)
            ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                            >  2:1001992/1‑76 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                           >  2:344479/1‑76 (MQ=255)
             tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                           <  2:454703/76‑1 (MQ=255)
              atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                          >  2:184457/1‑76 (MQ=255)
                 tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                       >  2:183258/1‑76 (MQ=255)
                   tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                     >  2:1114500/1‑76 (MQ=255)
                     ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                   <  1:312404/76‑1 (MQ=255)
                       ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                 <  1:347246/76‑1 (MQ=255)
                       ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                 >  2:1148503/1‑76 (MQ=255)
                       ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                 >  1:390932/1‑76 (MQ=255)
                       ttATTGTTCATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                 >  1:1246156/1‑76 (MQ=255)
                         tttGTTTTTTTCAGGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                               <  1:782061/75‑1 (MQ=255)
                         aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                               >  2:622771/1‑76 (MQ=255)
                          ttttttATTCCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                              <  1:1078277/73‑1 (MQ=255)
                          ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                              <  2:1228564/76‑1 (MQ=255)
                           tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             >  2:451743/1‑76 (MQ=255)
                            gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                            >  2:1198969/1‑76 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                           >  2:351239/1‑76 (MQ=255)
                             tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                           >  1:952222/1‑76 (MQ=255)
                              ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                          >  1:764711/1‑76 (MQ=255)
                               tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                         >  1:850462/1‑76 (MQ=255)
                               tATNTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                         >  1:1205188/1‑76 (MQ=255)
                                aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                        >  1:795228/1‑76 (MQ=255)
                                 tttCATGCNGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATagatga                                       >  1:1282151/1‑76 (MQ=255)
                                    cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGACTg                                    >  2:1251298/1‑76 (MQ=255)
                                     atccAGCCATTACACGATGAACGCTTATGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                   <  1:697942/73‑1 (MQ=255)
                                      tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                  >  2:75833/1‑76 (MQ=255)
                                      tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                  >  2:750546/1‑76 (MQ=255)
                                       gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                 >  2:1275995/1‑76 (MQ=255)
                                       gCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                 >  2:895947/1‑76 (MQ=255)
                                           ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                             <  2:424699/76‑1 (MQ=255)
                                            cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                            >  2:78621/1‑76 (MQ=255)
                                                 cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                       <  2:807063/76‑1 (MQ=255)
                                                  acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATNTTGGCCATGACTATTGATGAATGCGCGACCTGGATTg                      <  1:1309190/76‑1 (MQ=255)
                                                     aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                   <  2:488808/76‑1 (MQ=255)
                                                      agAACGCTTAAGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                  <  2:566443/75‑1 (MQ=255)
                                                        aaCGCTTATGAACTCCAGGCACTCCGCCATATTTTTTCCATGACTATTNATGAATGCGCGACCTGGATTGCCCAGa                <  1:1176649/76‑1 (MQ=255)
                                                          cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGCCCTGGATTGCCCAGAcc              >  1:431382/1‑76 (MQ=255)
                                                          cGCTTATGAACNCCNAGCACTCCGCCTTATTTTTGCCATGACTATTGATGAATGCGCTACCTGGATGGCCCAGAcc              <  2:203884/76‑1 (MQ=255)
                                                            cTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg            <  1:517798/76‑1 (MQ=255)
                                                                  aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg      <  1:1029361/76‑1 (MQ=255)
                                                                   aCTCCAGGCACTCCGCCCTATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt     >  2:471250/1‑76 (MQ=255)
                                                                   aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt     >  2:471260/1‑76 (MQ=255)
                                                                      ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa  >  2:497215/1‑76 (MQ=255)
                                                                          |                                                                       
AGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA  >  REL606/1747829‑1747974

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