Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,296,3290AC37.1% 131.0 / 31.8 70intergenic (‑149/‑66)inaA/yfaHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/1);  new base C (0/26);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 2.42e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTT  >  REL606/2296254‑2296403
                                                                           |                                                                          
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                            >  1:1404874/1‑76 (MQ=255)
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                            >  1:316049/1‑76 (MQ=255)
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                            >  1:668769/1‑76 (MQ=255)
   gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca                                                                         >  1:315679/1‑76 (MQ=255)
   gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca                                                                         >  2:343181/1‑76 (MQ=255)
    tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg                                                                        >  1:1020668/1‑76 (MQ=255)
     aTTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGc                                                                       >  2:215358/1‑76 (MQ=255)
      ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc                                                                      >  1:27890/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTCCCCTTCCTTTGCATCGCATAGATTCCACCAGccc                                                                     >  1:1527603/1‑75 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                     >  1:1527586/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                     >  1:1425180/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                     >  1:1115735/1‑76 (MQ=255)
       naaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                     >  1:229292/2‑76 (MQ=255)
        aaTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccac                                                                    >  2:262050/1‑76 (MQ=255)
         aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact                                                                   >  1:1818624/1‑76 (MQ=255)
               ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact                                                             >  2:904299/1‑76 (MQ=255)
               ttCTTTCGCCACAACTAGCGTGACCCCGTTTTTGCCCTTCATGTGCAGCGAATGGTTCCCACCAGccactaccact                                                             <  2:644189/76‑1 (MQ=255)
                tCTTTTGCCACTACCAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccacta                                                            >  1:1048006/1‑76 (MQ=255)
                 cTTTTGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTaa                                                           <  1:1191/76‑1 (MQ=255)
                  ttttGCCACTACTAGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAAt                                                          <  2:834969/76‑1 (MQ=255)
                    ttGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAg                                                        <  2:615462/76‑1 (MQ=255)
                     tGCCACTACTAGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGt                                                       <  1:1629975/76‑1 (MQ=255)
                     tGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGt                                                       <  1:1717731/76‑1 (MQ=255)
                      gCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTg                                                      >  1:1005922/1‑76 (MQ=255)
                       ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgc                                                     >  2:850191/1‑76 (MQ=255)
                         actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc                                                   >  2:1287766/1‑76 (MQ=255)
                           tactaGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGccc                                                 <  1:98644/76‑1 (MQ=255)
                           tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc                                                 >  1:1003535/1‑76 (MQ=255)
                           tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc                                                 >  2:1302181/1‑76 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                                >  2:391206/1‑76 (MQ=255)
                              taGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGt                                              <  2:1403201/76‑1 (MQ=255)
                               aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa                                             >  1:1207295/1‑76 (MQ=255)
                                gCGTGCCACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg                                            <  1:688734/76‑1 (MQ=255)
                                gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg                                            <  1:592526/76‑1 (MQ=255)
                                gCGTGACCCCGCTTTTCCCCTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg                                            <  1:500294/76‑1 (MQ=255)
                                gCGTGACACCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg                                            <  1:12791/76‑1 (MQ=255)
                                 cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa                                           >  1:482435/1‑76 (MQ=255)
                                  gTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGAt                                          >  2:1529863/1‑76 (MQ=255)
                                    gACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa                                        <  2:1421281/76‑1 (MQ=255)
                                    gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa                                        >  1:1932930/1‑76 (MQ=255)
                                     gcacCGTTTTAACCCTTAATTTGAATCGCATATACTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa                                       <  1:1224397/75‑1 (MQ=255)
                                     acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa                                       >  2:377813/1‑76 (MQ=255)
                                      cacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa                                      >  1:997919/1‑76 (MQ=255)
                                           tttttCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa                                 <  1:1467914/75‑1 (MQ=255)
                                           cTTTTCCCCTTCATGTGCAGCGAATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa                                 <  1:1963226/76‑1 (MQ=255)
                                             tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                               >  1:1273461/1‑76 (MQ=255)
                                             tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                               >  1:1511092/1‑76 (MQ=255)
                                              ttccccTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTATATAAAATGTGAAtt                              <  1:1071269/73‑1 (MQ=255)
                                               tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                             <  2:1163211/74‑1 (MQ=255)
                                                ccccTTCTTTTGCAGCGCAAAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa                            <  1:1421371/75‑1 (MQ=255)
                                                ccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa                            <  1:836419/75‑1 (MQ=255)
                                                 cccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa                           >  1:381403/1‑76 (MQ=255)
                                                  ccTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt                          <  2:805335/76‑1 (MQ=255)
                                                  ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt                          >  2:1147492/1‑76 (MQ=255)
                                                  ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGGGCGCCCAGTAGATAAAATGTGAATTTAAt                          >  1:996776/1‑76 (MQ=255)
                                                  ccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt                          <  2:1561586/76‑1 (MQ=255)
                                                    ttCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa                        <  1:1108550/76‑1 (MQ=255)
                                                     tCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                       <  2:1115735/76‑1 (MQ=255)
                                                     tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                       >  1:1976751/1‑76 (MQ=255)
                                                      cATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt                      >  1:780689/1‑76 (MQ=255)
                                                      cATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCact                      <  2:75829/76‑3 (MQ=255)
                                                        tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg                    >  1:1622462/1‑76 (MQ=255)
                                                        tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg                    >  1:1713334/1‑76 (MQ=255)
                                                         ttGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg                   <  2:1718063/76‑1 (MQ=255)
                                                          tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa                  >  1:85400/1‑76 (MQ=255)
                                                           gCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt                 <  1:1747522/76‑1 (MQ=255)
                                                           gCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt                 <  2:1429592/76‑1 (MQ=255)
                                                           gCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt                 <  1:622839/76‑1 (MQ=255)
                                                           gCAGCCCATAAATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGAt                 <  2:686153/76‑1 (MQ=255)
                                                            cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa                >  2:1544702/1‑76 (MQ=255)
                                                            cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa                >  1:524718/1‑76 (MQ=255)
                                                                gcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc            >  1:482443/1‑76 (MQ=255)
                                                                gcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc            >  1:1266167/1‑76 (MQ=255)
                                                                gcATAGATTCCACCAGCCACTACCACTAATACTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc            >  1:1968760/1‑76 (MQ=255)
                                                                 cATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGTATACAGCg           <  1:389554/76‑1 (MQ=255)
                                                                 cATAGACCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg           <  1:1270434/76‑1 (MQ=255)
                                                                  aTAGATTCCACCAGGCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa          >  1:1967882/1‑76 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt         >  1:1108948/1‑76 (MQ=255)
                                                                    aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg        >  2:1724229/1‑76 (MQ=255)
                                                                    aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg        >  1:348461/1‑76 (MQ=255)
                                                                     gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa       <  2:727295/76‑1 (MQ=255)
                                                                      aTTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGaa      >  2:1167959/1‑76 (MQ=255)
                                                                         ccaccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGtt   >  2:1509361/1‑76 (MQ=255)
                                                                         ccaccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGtt   >  1:12312/1‑76 (MQ=255)
                                                                         ccaccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGtt   >  1:1178487/1‑76 (MQ=255)
                                                                          caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGttt  >  2:20408/1‑76 (MQ=255)
                                                                          caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGttt  >  1:1525789/1‑76 (MQ=255)
                                                                           |                                                                          
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTT  >  REL606/2296254‑2296403

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