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

GGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747830‑1747975
                                                                         |                                                                        
ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa                                                                        >  1:1411066/1‑76 (MQ=255)
ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCa                                                                        <  2:473766/76‑1 (MQ=255)
 gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:709969/1‑76 (MQ=255)
 gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:287129/1‑76 (MQ=255)
 gCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAc                                                                       >  1:1414815/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:1169987/1‑75 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcc                                                                      >  1:125864/1‑75 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                      >  1:283985/1‑75 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAca                                                                      >  1:751529/1‑75 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:973321/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:1195646/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:1217153/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:104759/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:482439/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:458827/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACt                                                                      >  1:438109/1‑76 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                      <  2:547438/76‑1 (MQ=255)
  cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACt                                                                      >  1:1512574/1‑76 (MQ=255)
   ttaTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtg                                                                     >  1:1226412/3‑75 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcct                                                                     >  1:1807214/1‑74 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcca                                                                     >  1:1409621/1‑74 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                     >  1:20934/1‑74 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcac                                                                     >  1:1518547/1‑74 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtg                                                                     >  1:277735/1‑75 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     <  1:1953719/76‑1 (MQ=255)
   tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTc                                                                     >  1:1573620/1‑76 (MQ=255)
   tGATAATGCTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcat                                                                     >  1:1180619/1‑74 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                    >  1:483068/1‑73 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcctt                                                                    >  1:1670991/1‑73 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcact                                                                    >  1:1655394/1‑73 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCAcaca                                                                    >  1:381100/1‑73 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtag                                                                    >  1:1605138/1‑74 (MQ=255)
    gATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGcc                                                                    >  1:1037536/1‑76 (MQ=255)
     aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACttag                                                                   >  1:589606/1‑73 (MQ=255)
     aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac                                                                   >  1:1865686/1‑73 (MQ=255)
     aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACtgac                                                                   >  1:1892/1‑73 (MQ=255)
     aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   >  2:501126/1‑76 (MQ=255)
     aTAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCg                                                                   >  1:748659/1‑76 (MQ=255)
       aaTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACGccaat                                                                 >  1:580178/1‑73 (MQ=255)
        aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCa                                                                >  1:1969852/1‑76 (MQ=255)
          gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                              >  1:1690155/1‑76 (MQ=255)
          gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCata                                                              >  1:1202030/1‑76 (MQ=255)
           ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             >  2:881966/1‑76 (MQ=255)
           ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             <  2:1062570/76‑1 (MQ=255)
           ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCatat                                                             <  1:1336898/76‑1 (MQ=255)
            tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            >  1:764005/1‑76 (MQ=255)
            tatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtt                                                            >  1:390391/1‑76 (MQ=255)
             atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGTCACTTAGGCatatcg                                                           >  1:144377/1‑74 (MQ=255)
              tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          >  2:684155/1‑76 (MQ=255)
              tttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAtttt                                                          >  2:48791/1‑76 (MQ=255)
               ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATAttttt                                                         >  2:120039/1‑76 (MQ=255)
                tataTTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTg                                                        >  1:1650327/1‑76 (MQ=255)
                 atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                       >  1:6602/1‑76 (MQ=255)
                  tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      <  2:1219608/76‑1 (MQ=255)
                  tatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGcc                                                      >  1:1836221/1‑76 (MQ=255)
                   tttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                     <  1:1428162/75‑1 (MQ=255)
                     tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                   <  2:1347281/76‑1 (MQ=255)
                     tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                   <  2:1829929/76‑1 (MQ=255)
                     tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                   <  2:30807/76‑1 (MQ=255)
                      ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGa                                                  <  1:1310554/76‑1 (MQ=255)
                       tATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                 >  1:628729/1‑76 (MQ=255)
                        aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                <  2:69969/76‑1 (MQ=255)
                        aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                >  1:1400946/1‑76 (MQ=255)
                        aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                >  1:781890/1‑76 (MQ=255)
                        aTTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACt                                                >  1:131162/1‑76 (MQ=255)
                           ttttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                             <  1:302398/75‑1 (MQ=255)
                             ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  1:495957/1‑76 (MQ=255)
                             ttATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                           >  2:658631/1‑76 (MQ=255)
                               aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatg                                         >  2:1479323/1‑76 (MQ=255)
                                tttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                        >  1:318903/1‑76 (MQ=255)
                                 ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                       >  2:695394/1‑76 (MQ=255)
                                 ttCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGaa                                       <  2:1512573/76‑1 (MQ=255)
                                   cATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATg                                     >  1:1791197/1‑76 (MQ=255)
                                     tGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   <  1:518810/76‑1 (MQ=255)
                                     cgCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                   <  1:1478527/75‑1 (MQ=255)
                                        aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                >  1:1027550/1‑76 (MQ=255)
                                        aGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGa                                >  2:647625/1‑76 (MQ=255)
                                         gCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGGGCGAc                               >  2:739608/1‑76 (MQ=255)
                                         gCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAc                               >  1:1213483/1‑76 (MQ=255)
                                          ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              <  2:672128/76‑1 (MQ=255)
                                          ccATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                              >  1:1491357/1‑76 (MQ=255)
                                            aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                            >  1:1787249/1‑76 (MQ=255)
                                            aTTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTg                            >  1:501841/1‑76 (MQ=255)
                                               acacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                         >  2:543232/1‑76 (MQ=255)
                                                cacGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                        <  1:535023/76‑1 (MQ=255)
                                                      gAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAg                  >  1:1194032/1‑76 (MQ=255)
                                                         cGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc               >  1:1053497/1‑76 (MQ=255)
                                                          gCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg              >  1:618478/1‑76 (MQ=255)
                                                             tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg           <  2:1737044/76‑1 (MQ=255)
                                                             tATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg           >  1:1274185/1‑76 (MQ=255)
                                                                gAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa        <  1:1970375/76‑1 (MQ=255)
                                                                   cTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTg     <  1:29760/76‑1 (MQ=255)
                                                                    cccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa    <  1:1424138/75‑1 (MQ=255)
                                                                     ccAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaa   >  1:1180043/1‑76 (MQ=255)
                                                                      cgggCCCTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa  <  1:1520611/74‑1 (MQ=255)
                                                                      cAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGaaa  <  1:555176/76‑1 (MQ=255)
                                                                         |                                                                        
GGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAA  >  REL606/1747830‑1747975

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: