Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT27.8% 228.0 / 68.4 90A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/23);  new base T (25/0);  total (67/23)
Fisher's exact test for biased strand distribution p-value = 2.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  REL606/1747828‑1747978
                                                                           |                                                                           
cagGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAgg                                                                             >  2:1546871/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                           >  1:310106/1‑76 (MQ=255)
  ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCa                                                                           >  1:1346311/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                          >  1:174514/1‑76 (MQ=255)
   gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCAc                                                                          <  1:1217688/76‑1 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                         >  1:470153/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                         >  1:371889/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                         >  1:694984/1‑75 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                         >  1:126599/1‑76 (MQ=255)
    cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                         <  1:2143450/76‑1 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                        >  1:137880/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                        >  1:139512/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                        >  1:1505427/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                        >  1:1879921/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                        >  1:528153/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                        >  1:908397/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                        >  1:1874465/1‑74 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtt                                                                        >  1:1457952/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtg                                                                        >  1:380357/1‑75 (MQ=255)
     tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                        >  1:2080158/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                       >  1:921433/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcatc                                                                       >  1:137191/1‑73 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                       >  1:2011244/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                       >  1:282924/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                       >  1:635372/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                       >  1:1492277/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                       >  1:1430583/1‑74 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                       >  1:2143639/1‑76 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                       <  1:1569787/76‑1 (MQ=255)
      gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                       >  1:30378/1‑76 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                      >  1:1796471/1‑73 (MQ=255)
       aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                      >  2:1122391/1‑76 (MQ=255)
          aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                   >  1:937657/1‑76 (MQ=255)
           tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                                  <  2:550231/76‑1 (MQ=255)
             ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                                >  1:1547834/1‑76 (MQ=255)
              tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                               >  1:465121/1‑76 (MQ=255)
               atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                              >  1:432062/1‑76 (MQ=255)
                 ttatATTTTATTTTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                            <  1:89707/76‑1 (MQ=255)
                 ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                            >  1:1525104/1‑76 (MQ=255)
                  tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                           >  1:1346592/1‑76 (MQ=255)
                     atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                        >  2:1869632/1‑76 (MQ=255)
                      ttttttttttTATTTCATGCAGCCATTACACGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCAt                                                       <  1:968732/76‑1 (MQ=255)
                       tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                      >  1:942080/1‑76 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                    >  1:513919/1‑76 (MQ=255)
                         tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                    >  1:485472/1‑76 (MQ=255)
                          aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                   >  1:365856/1‑76 (MQ=255)
                             tttttttTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                <  1:441720/75‑1 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                >  2:564260/1‑76 (MQ=255)
                             gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                >  1:1896823/1‑76 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               <  2:1729385/76‑1 (MQ=255)
                              tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               <  1:1932575/76‑1 (MQ=255)
                                tATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                             >  2:832974/1‑76 (MQ=255)
                                 ttttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                            <  1:1243097/75‑1 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                            >  2:1352565/1‑76 (MQ=255)
                                 aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                            >  1:681751/1‑76 (MQ=255)
                                  tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                           <  2:2067600/76‑1 (MQ=255)
                                  tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                           >  2:1752234/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                          <  1:545931/76‑1 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                          >  2:1483497/1‑76 (MQ=255)
                                   ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                          <  2:574055/76‑1 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                        >  1:1154793/1‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                        >  1:332394/1‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                        >  2:1120700/1‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                        >  2:1782836/1‑76 (MQ=255)
                                     cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                        >  2:433708/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                    >  2:496360/1‑76 (MQ=255)
                                         cAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                    >  2:1605639/1‑76 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                <  1:387060/76‑1 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                >  1:1924817/1‑76 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                >  1:15452/1‑76 (MQ=255)
                                             cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                >  1:555813/1‑76 (MQ=255)
                                               ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                              >  1:823274/1‑76 (MQ=255)
                                                tACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                             >  1:1615638/1‑76 (MQ=255)
                                                 acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                            >  2:2142203/1‑76 (MQ=255)
                                                  cacGATGAACGCTTATGAACTTCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                           <  1:1329058/76‑1 (MQ=255)
                                                  cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                           <  1:1329056/76‑1 (MQ=255)
                                                    cGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGc                         <  1:427224/76‑1 (MQ=255)
                                                      aTGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                       <  2:1699607/76‑1 (MQ=255)
                                                        gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGGATGCGCGACCTGGATTGCCCAg                     <  1:1525127/76‑1 (MQ=255)
                                                           cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                  <  2:2122733/76‑1 (MQ=255)
                                                           cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                  >  1:1372186/1‑76 (MQ=255)
                                                               tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg              >  1:2032941/1‑76 (MQ=255)
                                                                 tGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc            <  2:1791716/76‑1 (MQ=255)
                                                                 tGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc            <  2:550060/76‑1 (MQ=255)
                                                                   aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg          >  1:404022/1‑76 (MQ=255)
                                                                   aaCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg          >  1:2068378/1‑76 (MQ=255)
                                                                    aCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt         <  2:2163923/76‑1 (MQ=255)
                                                                      tCCGGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa       >  2:1630884/1‑76 (MQ=255)
                                                                      tCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa       >  2:1630872/1‑76 (MQ=255)
                                                                      cccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa       <  2:897857/75‑1 (MQ=255)
                                                                       ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa      <  1:601950/76‑1 (MQ=255)
                                                                       ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa      >  2:773119/1‑76 (MQ=255)
                                                                       ccAGCCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa      <  1:601970/76‑1 (MQ=255)
                                                                         aGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg    >  1:403103/1‑76 (MQ=255)
                                                                           gCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  <  1:1743320/76‑1 (MQ=255)
                                                                           |                                                                           
CAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  REL606/1747828‑1747978

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