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

GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA  >  REL606/2296254‑2296404
                                                                           |                                                                           
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                             >  2:818767/1‑76 (MQ=255)
gtCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTcca                                                                             <  1:1474078/76‑1 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                      >  1:374197/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                      >  1:374191/1‑76 (MQ=255)
       tAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGcca                                                                      >  1:1365327/1‑76 (MQ=255)
         aTGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccact                                                                    >  2:557285/1‑76 (MQ=255)
            aaTTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactacc                                                                 >  2:1187390/1‑76 (MQ=255)
              tttCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGccactaccac                                                               >  1:660572/1‑76 (MQ=255)
                 cTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTaa                                                            >  2:1248193/1‑76 (MQ=255)
                  ttttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAAt                                                           >  1:1402565/1‑76 (MQ=255)
                    ttGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAg                                                         >  2:253519/1‑76 (MQ=255)
                     tGCCACTACTAGCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGt                                                        <  2:187849/76‑1 (MQ=255)
                       ccACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACCACCACTAATAGTgc                                                      >  1:1150805/1‑76 (MQ=255)
                        cACTACTAGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTgcg                                                     <  2:1167830/76‑1 (MQ=255)
                        cACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTgcg                                                     >  1:1058757/1‑76 (MQ=255)
                          ctactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGcc                                                   >  2:312435/1‑76 (MQ=255)
                           tactaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGccc                                                  >  2:1224994/1‑76 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                                 >  1:848920/1‑76 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                                 >  1:1078048/1‑76 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                                 >  1:25214/1‑76 (MQ=255)
                            actaGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCa                                                 >  1:1158952/1‑76 (MQ=255)
                               aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTa                                              <  1:141235/76‑1 (MQ=255)
                               aGCGTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCACCAGCCACTACCACTAATAGTGCGCCCAGTa                                              <  2:268767/76‑1 (MQ=255)
                               aGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTa                                              >  1:65713/1‑76 (MQ=255)
                                gCGTGCCCCCGCTTTTACCCTTCATTTGCACCCAATGGATCCCACCAGCCACTGCCACTAATAGTGCGCCCAGTAg                                             <  2:1312752/76‑1 (MQ=255)
                                gCGTGACACCGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAg                                             <  2:945075/76‑1 (MQ=255)
                                 cGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGa                                            >  1:1169248/1‑76 (MQ=255)
                                  gTGACCCCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGAt                                           <  2:1127149/76‑1 (MQ=255)
                                   tGACACCGCTTTTACCCTTCATTTGCCGCGCATAGATTCCACCAGCCCCTACCACTAATAGTGCGCCCAGtatata                                          >  2:647961/1‑76 (MQ=255)
                                   tGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATa                                          >  2:677885/1‑76 (MQ=255)
                                    gACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaa                                         >  1:292141/1‑76 (MQ=255)
                                     acacCGCTTTTCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa                                        <  1:1117650/76‑1 (MQ=255)
                                     acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTCCCACTAATAGTGCGCCCAGTAGATaaa                                        >  1:799065/1‑76 (MQ=255)
                                     acacCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaa                                        >  1:799049/1‑76 (MQ=255)
                                      ccccGCTTTTCCCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATaaaa                                       <  1:1245720/74‑1 (MQ=255)
                                        ccGTTTTTCCCCTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg                                     <  1:955610/76‑1 (MQ=255)
                                        ccGCTTTTACCCTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg                                     <  2:754493/76‑1 (MQ=255)
                                        ccGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg                                     >  2:632883/1‑76 (MQ=255)
                                        ccGCTTTTACCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtg                                     <  1:350153/76‑1 (MQ=255)
                                         cGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgt                                    >  1:790/1‑76 (MQ=255)
                                          gCTTTTCCCCTTCATTTGCAGCGAATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAAtgtg                                   <  2:1319489/76‑1 (MQ=255)
                                            ttttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa                                 >  1:769725/1‑76 (MQ=255)
                                            ttttACCCCTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGaa                                 >  1:1181439/1‑76 (MQ=255)
                                             tttCCCCTTCATTTGCAGCGCATAGATCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAt                                <  2:223630/76‑1 (MQ=255)
                                              ttACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAtt                               >  1:1279346/1‑76 (MQ=255)
                                               tccccTTCATTGGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                              <  1:1052812/74‑1 (MQ=255)
                                               tACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAAttt                              >  2:807251/1‑76 (MQ=255)
                                                aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTa                             >  1:781710/1‑76 (MQ=255)
                                                aCCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGACCAGTAGATAAAATGTGAATTTa                             >  1:169070/1‑76 (MQ=255)
                                                  ccTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAAt                           >  1:1271044/1‑76 (MQ=255)
                                                   cTTCATTTGCAGCGCATAGATTCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATc                          <  1:327857/76‑1 (MQ=255)
                                                     tCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAg                        >  1:1503447/1‑76 (MQ=255)
                                                        tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg                     >  1:555059/1‑76 (MQ=255)
                                                        tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg                     >  2:939219/1‑76 (MQ=255)
                                                        tttGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCg                     >  2:939203/1‑76 (MQ=255)
                                                               cgcATAGCTCCCCCCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAg              <  1:138162/76‑1 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt          >  2:65791/1‑76 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt          >  2:65802/1‑76 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt          >  2:1025314/1‑76 (MQ=255)
                                                                   tAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGAt          >  1:1244307/1‑76 (MQ=255)
                                                                     gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa        >  1:1344988/1‑76 (MQ=255)
                                                                     gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa        >  1:1344993/1‑76 (MQ=255)
                                                                     gATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGa        >  1:1521724/1‑76 (MQ=255)
                                                                        tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt     >  1:751903/1‑76 (MQ=255)
                                                                        tCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGt     >  2:973108/1‑76 (MQ=255)
                                                                          caccaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGttt   >  1:454288/1‑76 (MQ=255)
                                                                           accaGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTa  >  2:1493749/1‑76 (MQ=255)
                                                                           |                                                                           
GTCGTATTAATGAATTTCTTTTGCCACTACTAGCGTGACACCGCTTTTACCCTTCATTTGCAGCGCATAGATTCCACCAGCCACTACCACTAATAGTGCGCCCAGTAGATAAAATGTGAATTTAATCAGTCGGATACAGCGATGAAGTTTA  >  REL606/2296254‑2296404

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