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

GTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAG  >  REL606/2296257‑2296400
                                                                        |                                                                       
gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca                                                                      <  2:327146/76‑1 (MQ=255)
 tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg                                                                     >  2:426515/1‑76 (MQ=255)
 tATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAg                                                                     >  2:1473934/1‑76 (MQ=255)
  aTTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGc                                                                    >  1:1493744/1‑76 (MQ=255)
   ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc                                                                   >  2:720408/1‑76 (MQ=255)
    tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                  >  2:1170485/1‑76 (MQ=255)
      aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact                                                                >  2:496562/1‑76 (MQ=255)
       tGAATTTGTTTTGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGccacta                                                               <  1:1632419/76‑1 (MQ=255)
         aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacc                                                             >  2:443602/1‑76 (MQ=255)
         aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacc                                                             >  2:1308996/1‑76 (MQ=255)
            ttCTTTTGCCACTACTAGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGccactaccact                                                          <  2:871852/76‑1 (MQ=255)
             tCTTTTGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGccactaccacta                                                         <  2:503793/76‑1 (MQ=255)
                tttGCCACTACTAGCGTGACCCCGTTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATa                                                      <  2:1579265/76‑1 (MQ=255)
                tttGCCACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATa                                                      <  2:382435/76‑1 (MQ=255)
                 ttGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAg                                                     <  1:1170140/76‑1 (MQ=255)
                  tGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGt                                                    >  2:932986/1‑76 (MQ=255)
                   gCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTg                                                   >  1:85445/1‑76 (MQ=255)
                     cACTACTAGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTgcg                                                 <  1:1179496/76‑1 (MQ=255)
                     cACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgcg                                                 <  1:1307499/76‑1 (MQ=255)
                      actactAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcgc                                                >  1:117201/1‑76 (MQ=255)
                        tactaGCGTGCCACCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc                                              <  1:1401061/76‑1 (MQ=255)
                           taGCGTGCCCCCGTTTTTCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGt                                           <  1:332001/76‑1 (MQ=255)
                           taGCGTGACACCGCTTTTATCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt                                           >  1:1403806/1‑76 (MQ=255)
                            aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa                                          >  1:812665/1‑76 (MQ=255)
                             gCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTan                                         <  1:1472527/76‑2 (MQ=255)
                              cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa                                        >  1:94384/1‑76 (MQ=255)
                              cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa                                        >  1:7984/1‑76 (MQ=255)
                               gTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt                                       <  2:781222/76‑1 (MQ=255)
                               gTGACACCGCTTTTACCCTACATTTGCGGCGCATATATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt                                       <  1:1580790/76‑1 (MQ=255)
                                tGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATa                                      <  1:1266230/76‑1 (MQ=255)
                                 gACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa                                     <  2:455931/76‑1 (MQ=255)
                                      cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt                                >  2:1360162/1‑76 (MQ=255)
                                        tttttCCCCTTGATTGGCAGCGCATAGATCCCGCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa                              <  1:1388459/75‑1 (MQ=255)
                                        cTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa                              >  2:1537222/1‑76 (MQ=255)
                                        cTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGa                              >  1:1169146/1‑76 (MQ=255)
                                         ttttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa                             >  1:1724296/1‑76 (MQ=255)
                                          tttCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                            <  2:1677188/76‑1 (MQ=255)
                                            tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                          <  1:870234/74‑1 (MQ=255)
                                              cccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa                        >  1:1679046/1‑76 (MQ=255)
                                               ccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt                       <  2:1505135/76‑1 (MQ=255)
                                                cTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc                      >  2:1161031/1‑76 (MQ=255)
                                                cTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc                      <  1:1327355/76‑1 (MQ=255)
                                                cTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc                      <  1:180735/76‑1 (MQ=255)
                                                  tcttttGCAGCGCAAAAATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                    <  2:386170/73‑1 (MQ=255)
                                                  tCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                    <  1:28517/76‑1 (MQ=255)
                                                    aTTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTc                  <  2:1592675/76‑1 (MQ=255)
                                                      ttGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg                <  1:405894/76‑1 (MQ=255)
                                                      tggCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCgg                <  2:223687/74‑1 (MQ=255)
                                                         cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa             >  2:1523730/1‑76 (MQ=255)
                                                         cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa             >  1:1035596/1‑76 (MQ=255)
                                                           gcgcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa           >  2:405105/1‑76 (MQ=255)
                                                           gcgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa           <  2:1001907/76‑1 (MQ=255)
                                                           gcgaATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa           <  1:1521641/76‑1 (MQ=255)
                                                           gcccATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa           <  1:1460798/73‑1 (MQ=255)
                                                            cgcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg          <  2:307257/76‑1 (MQ=255)
                                                             gcATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc         >  1:1429483/1‑76 (MQ=255)
                                                             gcATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc         <  1:1399202/76‑1 (MQ=255)
                                                             ccATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGc         <  1:406387/75‑1 (MQ=255)
                                                              cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTGGGATACAGCg        >  1:138770/1‑76 (MQ=255)
                                                               aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa       >  2:743261/1‑76 (MQ=255)
                                                               aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa       >  2:794482/1‑76 (MQ=255)
                                                                tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt      >  2:1254002/1‑76 (MQ=255)
                                                                tAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt      <  1:1194536/76‑1 (MQ=255)
                                                                aggATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt      <  1:1194549/74‑1 (MQ=255)
                                                                 aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg     >  1:635942/1‑76 (MQ=255)
                                                                 aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCACTCGGATACAGCGATg     >  1:1483383/1‑76 (MQ=255)
                                                                  gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa    >  1:826447/1‑76 (MQ=255)
                                                                  gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa    <  1:785559/76‑1 (MQ=255)
                                                                  gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa    <  1:1135584/76‑1 (MQ=255)
                                                                    ttCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAg  >  2:924438/1‑76 (MQ=255)
                                                                        |                                                                       
GTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAG  >  REL606/2296257‑2296400

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: ‑