Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,296,3290AC33.9% 111.5 / 25.3 59intergenic (‑149/‑66)inaA/yfaHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/3);  new base C (0/20);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.34e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGT  >  REL606/2296254‑2296401
                                                                           |                                                                        
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                          >  1:1434458/1‑76 (MQ=255)
  cGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacc                                                                        <  1:1726776/76‑1 (MQ=255)
   gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca                                                                       <  2:201621/76‑1 (MQ=255)
   gTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTccacca                                                                       >  1:938856/1‑76 (MQ=255)
      ttAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcc                                                                    >  2:1314874/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                   >  2:17246/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                   <  1:235615/76‑1 (MQ=255)
         aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact                                                                 >  2:1730656/1‑76 (MQ=255)
           gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac                                                               >  1:1407134/1‑76 (MQ=255)
           gAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactac                                                               >  1:1407145/1‑76 (MQ=255)
               tttttttCCCACTCCTAGCGTGCCCCCGTTTTTCCCTTTCATTTGCAGCGCATAGATCCCCCCAGccactaccact                                                           <  1:842961/73‑1 (MQ=255)
               ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact                                                           >  1:1687556/1‑76 (MQ=255)
               ttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccact                                                           >  1:1400693/1‑76 (MQ=255)
                 cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa                                                         >  1:946350/1‑76 (MQ=255)
                   tttGCCACTACTAGCGTGCCCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATa                                                       <  2:1243679/76‑1 (MQ=255)
                   tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCCTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa                                                       >  1:994263/1‑76 (MQ=255)
                   tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa                                                       >  2:1357913/1‑76 (MQ=255)
                   tttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATa                                                       >  1:518263/1‑76 (MQ=255)
                    ttcccACTACTAGCGTGACCCCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAg                                                      <  2:1139672/73‑1 (MQ=255)
                       ccACTACTAGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTgc                                                   <  1:245770/76‑1 (MQ=255)
                          ctactaGGGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGcc                                                <  1:1264036/76‑1 (MQ=255)
                          ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc                                                >  1:116075/1‑76 (MQ=255)
                          ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc                                                >  2:525764/1‑76 (MQ=255)
                           tcctaGCGTCACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc                                               <  1:348658/74‑1 (MQ=255)
                           tactaGCGTGACCCCGCTTTTACCCTTCATTGGCAGCGCAAAAATCCCCCCAGCCACTACCACTAATAGTGCGccc                                               <  1:1288494/76‑1 (MQ=255)
                           tactaGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGccc                                               <  1:1606601/76‑1 (MQ=255)
                            ccaaGCGTGCCCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCa                                              <  2:1178754/73‑1 (MQ=255)
                            actgGCGTGACCCCGTTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCa                                              <  2:1449438/76‑1 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                              >  1:1360670/1‑76 (MQ=255)
                             ctaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAg                                             >  2:806673/1‑76 (MQ=255)
                              taGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGt                                            >  2:411228/1‑76 (MQ=255)
                               aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa                                           >  1:1439308/1‑76 (MQ=255)
                                gCGTGACCCCCTTTTTCCCCTTCATTTGCAGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAg                                          <  1:781543/76‑1 (MQ=255)
                                  gggACACCGTTTTTCCCCTTCATTTCCACCGCATATATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGAt                                        <  2:773493/74‑1 (MQ=255)
                                   tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa                                       >  1:1905518/1‑76 (MQ=255)
                                    gACCCCCCTTTTCCCCTTCATTGGCGGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa                                      <  2:1543639/76‑1 (MQ=255)
                                     acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa                                     >  2:795766/1‑76 (MQ=255)
                                       cccGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAt                                   <  2:675960/75‑1 (MQ=255)
                                       cccGCTTTTCCCCTTCATTGGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAt                                   <  2:1754331/75‑1 (MQ=255)
                                             tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                             >  1:1519127/1‑76 (MQ=255)
                                             tttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                             >  1:1495908/1‑76 (MQ=255)
                                              ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGAAGATAAAATGTGAAtt                            >  1:831251/1‑76 (MQ=255)
                                               tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                           <  2:1288893/74‑1 (MQ=255)
                                               tccccTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                           <  1:1317849/74‑1 (MQ=255)
                                               tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACAAATAGTGCGCCCAGTAGATAAAATGTGAAttt                           >  1:1735134/1‑76 (MQ=255)
                                                 cccTTCAGTTGCAGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTaa                         <  2:497845/76‑1 (MQ=255)
                                                    ttCATTTGCAGCGCATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa                      <  2:962728/76‑1 (MQ=255)
                                                    ttCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCa                      <  2:962795/76‑1 (MQ=255)
                                                     tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                     >  1:491514/1‑76 (MQ=255)
                                                     tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                     >  2:493370/1‑76 (MQ=255)
                                                      naTTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt                    >  1:1283133/2‑76 (MQ=255)
                                                      cATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGt                    >  1:85719/1‑76 (MQ=255)
                                                          tGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGa                >  1:158911/1‑76 (MQ=255)
                                                            cAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATa              >  1:1584951/1‑76 (MQ=255)
                                                              gcgcATAGATCCCCCAAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACa            <  1:958644/76‑1 (MQ=255)
                                                                 cATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCg         >  1:1220089/1‑76 (MQ=255)
                                                                  aTAGTTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa        <  2:521504/76‑1 (MQ=255)
                                                                  aTAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa        >  1:1693303/1‑76 (MQ=255)
                                                                  aTAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGa        <  1:1929141/76‑1 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt       >  2:930411/1‑76 (MQ=255)
                                                                    aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg      >  1:1085392/1‑76 (MQ=255)
                                                                    aGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg      >  2:1317185/1‑76 (MQ=255)
                                                                    aGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATg      <  2:956013/76‑1 (MQ=255)
                                                                     gATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa     <  1:1652658/76‑1 (MQ=255)
                                                                      aTTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGaa    >  2:1533494/1‑76 (MQ=255)
                                                                       ttCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAg   >  1:1264003/1‑76 (MQ=255)
                                                                        tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt  >  1:163489/1‑76 (MQ=255)
                                                                           |                                                                        
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGT  >  REL606/2296254‑2296401

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