Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,747,9030GT28.4% 223.0 / 68.8 88A8S (GCA→TCA) ydiLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/25);  new base T (25/0);  total (63/25)
Fisher's exact test for biased strand distribution p-value = 4.93e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA  >  REL606/1747831‑1747974
                                                                        |                                                                       
gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:545459/1‑76 (MQ=255)
gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:2070394/1‑76 (MQ=255)
gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:2018408/1‑76 (MQ=255)
gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                      >  1:1981300/1‑76 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                     >  1:361648/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcg                                                                     >  1:321676/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                     >  1:434624/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                     >  1:2016684/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                     >  1:1607781/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                     >  1:769827/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                     >  1:1875386/1‑75 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:935365/1‑76 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:747244/1‑76 (MQ=255)
 cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                     >  1:331642/1‑76 (MQ=255)
  tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgc                                                                    >  1:1053030/1‑74 (MQ=255)
  tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                    >  1:1398085/1‑74 (MQ=255)
  tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                    >  1:387615/1‑74 (MQ=255)
  tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                    >  1:1807960/1‑76 (MQ=255)
  tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATAAACGCTTATGAACTCCAGTCAcca                                                                    >  1:1219763/1‑74 (MQ=255)
   gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcgat                                                                   >  1:1848667/1‑73 (MQ=255)
   gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcact                                                                   >  1:967114/1‑73 (MQ=255)
   gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtta                                                                   >  1:1203207/1‑74 (MQ=255)
   gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                   >  1:1377929/1‑76 (MQ=255)
   gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTcc                                                                   >  1:898201/1‑76 (MQ=255)
    aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCCTTACACGATGAACGCTTATGAACTACAGTCACTCCg                                                                  >  2:1502513/1‑76 (MQ=255)
    aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcca                                                                  >  1:283772/1‑75 (MQ=255)
    aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                  >  2:106753/1‑76 (MQ=255)
     tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaa                                                                 >  1:280031/1‑74 (MQ=255)
     tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                 >  2:1177603/1‑76 (MQ=255)
     tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                 >  2:917942/1‑76 (MQ=255)
     tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGc                                                                 <  2:1816351/76‑1 (MQ=255)
      aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                >  1:1832016/1‑73 (MQ=255)
        tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                              <  1:141614/76‑1 (MQ=255)
        tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                              <  1:200236/76‑1 (MQ=255)
        tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCat                                                              <  2:561321/76‑1 (MQ=255)
          ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                            >  2:2068639/1‑76 (MQ=255)
          ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                            >  1:75877/1‑76 (MQ=255)
           ttgttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                           <  1:1171435/73‑1 (MQ=255)
            atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                          >  1:90527/1‑76 (MQ=255)
            atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttt                                                          >  1:1779510/1‑76 (MQ=255)
             tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                         <  2:215645/76‑1 (MQ=255)
              ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                        >  1:1324056/1‑76 (MQ=255)
                ntatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                      >  1:5928/2‑76 (MQ=255)
                atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                      >  2:1793708/1‑76 (MQ=255)
                atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                      >  2:1984676/1‑76 (MQ=255)
                atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                      >  2:1557016/1‑76 (MQ=255)
                 tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                     <  2:941622/76‑1 (MQ=255)
                 tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                     <  2:1214695/76‑1 (MQ=255)
                  atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                    >  1:1909456/1‑76 (MQ=255)
                    tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                  <  1:1363716/76‑1 (MQ=255)
                    tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                  >  1:816639/1‑76 (MQ=255)
                      tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                >  2:475545/1‑76 (MQ=255)
                      tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                >  2:1188009/1‑76 (MQ=255)
                      tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                >  2:1837692/1‑76 (MQ=255)
                       aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                               >  2:1180272/1‑76 (MQ=255)
                        ttttttATTTCTTGCAGCCATTACACGATGAACGCTTATGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTa                                              <  1:1710552/73‑1 (MQ=255)
                        ttGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                              >  1:771236/1‑76 (MQ=255)
                         tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             >  1:1166133/1‑76 (MQ=255)
                         tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             <  2:1017059/76‑1 (MQ=255)
                         tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             >  2:767881/1‑76 (MQ=255)
                         tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             >  1:341572/1‑76 (MQ=255)
                         tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                             >  2:1858327/1‑76 (MQ=255)
                          gTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                            >  1:329987/1‑76 (MQ=255)
                           tttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                           >  2:1102157/1‑76 (MQ=255)
                            ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                          >  2:1748238/1‑76 (MQ=255)
                            ttATTTCATGCAGCCATCACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                          >  1:161621/1‑76 (MQ=255)
                              aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                        >  2:9742/1‑76 (MQ=255)
                              aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                        >  2:830459/1‑76 (MQ=255)
                              aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                        >  1:1395729/1‑76 (MQ=255)
                                ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                      <  2:966245/76‑1 (MQ=255)
                                ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                      >  1:1536792/1‑76 (MQ=255)
                                ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                      <  1:1807546/76‑1 (MQ=255)
                                ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                      >  1:1967652/1‑76 (MQ=255)
                                   atccAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                   <  1:1205854/73‑1 (MQ=255)
                                    tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                  >  1:1670330/1‑76 (MQ=255)
                                          cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                            >  1:95989/1‑76 (MQ=255)
                                          cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                            <  1:367104/76‑1 (MQ=255)
                                          cATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                            <  1:1953294/76‑1 (MQ=255)
                                            ttACACGATGAACGTTTATAAACTCCGGGCTCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                          <  1:467559/76‑1 (MQ=255)
                                            ttACACGATGAACGCTTATGACCCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                          <  1:1137650/76‑1 (MQ=255)
                                            ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                          <  1:1405177/76‑1 (MQ=255)
                                            ttACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTgg                          <  1:1432971/76‑1 (MQ=255)
                                             tACACGATGAACGCTTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                         <  1:1633082/76‑1 (MQ=255)
                                             tACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGa                         <  1:1591915/76‑1 (MQ=255)
                                               cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                       <  2:1395729/76‑1 (MQ=255)
                                               cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                       <  2:496559/76‑1 (MQ=255)
                                                acGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                      >  1:1936678/1‑76 (MQ=255)
                                                 cGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGc                     <  2:1685980/76‑1 (MQ=255)
                                                    tGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                  <  2:1214407/76‑1 (MQ=255)
                                                     gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                 <  1:49823/76‑1 (MQ=255)
                                                       aCGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAc               <  2:1795817/76‑1 (MQ=255)
                                                        cGCTTATGATCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc              <  1:1355147/76‑1 (MQ=255)
                                                              tGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc        <  2:1346517/76‑1 (MQ=255)
                                                              tGAACCCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc        <  1:2007984/76‑1 (MQ=255)
                                                               gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa       <  2:1575060/76‑1 (MQ=255)
                                                                    ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa  <  2:1943253/76‑1 (MQ=255)
                                                                    ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa  >  1:1323074/1‑76 (MQ=255)
                                                                        |                                                                       
GCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAA  >  REL606/1747831‑1747974

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