Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,693,1400CG38.7% 48.7 / 11.9 32P256P (CCG→CCC) unknownputative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (15/4);  new base G (0/12);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.61e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCAGCAGTTCTGGCTGCAGCATTTGCCATGGTAGGATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTACCTCGGTAAGTAG  >  1/1692896‑1693382
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
cGGCAGCAGTTCTGGCTGCAGCATTTGCCATGGTCGGATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATaa                                                                                                                                                                                                                                              >  1:228641/1‑251 (MQ=255)
      cagTTCTGGCTGCAGCATTTGCCATGGTCGGATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCg                                                                                                                                                                                                                                           <  1:222425/248‑1 (MQ=255)
                         tGCCATGGTCGGATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTc                                                                                                                                                                                                                                                    >  1:707980/1‑220 (MQ=255)
                           ccgtggtggGATCAGGTTTGCCCCCGTAACTGGGATAATTTTTTCTTTTGTATACATGGTTTTTTGCG‑GGGGGGAGTAGTTTTCTAGGTGGGGTTTCATTTTTATTACCGGCCGATTAAAAAAAGCCACTGCTTTGG‑GCGAAAA‑GT‑TTTTGTCGCATATCGGCCAAAAA‑A‑CGTGGTCTTAACCCGGCTTCAGCATTATCCCCCCACTGACCACAGGGAGAATAACCGCCATTACCGTTGTTTGGTAAccc                                                                                                                                                                                                                    <  1:239495/247‑1 (MQ=255)
                                    gATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCg                                                                                                                                                                                                                                  >  1:402411/1‑227 (MQ=255)
                                                              taagtttttCTTTGGGATACGTCAATATTTCCC‑GGGGGGGGAAGTTTTCTAAGGAGGGTTTCATTTAAATTATCGACCGAGTTAAAAAAGCCTCTGACTTGT‑GGAGAAT‑GT‑TTTTGTCGCGTTTCGGCAAAAAA‑A‑CGTGTTCTTAACCCGGCTTCAGTATTATTCCCCCACAAACCACAAGGGGGATAGCCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAAt                                                                                                                                                                                <  1:128850/246‑1 (MQ=255)
                                                                    tttCTTTTGGATACATCG‑TTTTTTCCGGGGGGGAGTGGTTTTCTCAGTGGGGTTTAATTTTGTTTACCGACCGAGTAAAAAAAGTCACTG‑TTTGG‑GGGGAAA‑GGTTTTTGTCGCATACCGGCAAAAAA‑G‑CGTGGTTTTAATCCGGCTTCAGCTTTTTCCCCCCCCTGATCAAAAGTGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCtg                                                                                                                                                                          <  1:742603/251‑1 (MQ=255)
                                                                               tACATCGATATTTCCG‑GGGCGGAGTAGTGTTCTCAGGGGGGTTTCTTTTTAATTATCGACCGATTAAAAAAAGTCACTGATTTG‑‑GGGGAAAGGG‑TTTTGTCGCATACCGCCAAAAAA‑A‑CGTGGTTTAAACCCGGCTTCAGCTTTATCCCCCCCCTGATCACAAGTGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGAt                                                                                                                                                               <  1:79316/251‑1 (MQ=255)
                                                                                                  gCGGAGTAGTGTTCTCAGTAGGGCCTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCa                                                                                                                                              >  1:307162/1‑250 (MQ=255)
                                                                                                     gAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTAt                                                                                                                                          >  1:229942/1‑251 (MQ=255)
                                                                                                       gtagtGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTAttt                                                                                                                                        >  1:324279/1‑251 (MQ=255)
                                                                                                          tttttctcggtgggttttctttttGATTATCGACGGATTAAAAAAAGTCACTGTCTTGT‑GGGGAAA‑GT‑TTTTGTCGCATACCGGCCAATAA‑A‑CGTGGTTTTACTCCGGTTTCAGTTTTTTCCCCCCACTGATCACAGGGGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTAttttt                                                                                                                                      <  1:513443/228‑1 (MQ=255)
                                                                                                           ttttctcagtgggttttctttttGATTATCGACCGATTAAAAAAATTCACTGTTTTGT‑GGGGAAA‑GT‑TTTTGTCGCATAGCGGCAAATAA‑A‑GGGGTTTTTACCCCGTCTTCAGCTTTTTCCCCCCACTGACCAAAAGTGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTtaa                                                                                                                                    <  1:895013/230‑1 (MQ=255)
                                                                                                            tttttcagtgggtcttcttttttaTTATCGCCGGAGTAAAAAAAGTCCCTGTTTTG‑‑GGGGAAA‑GGTTTTTGTCGCATA‑CGGGCAAAAAGG‑CGTGGTTTTAACCCGGTTTCAGTTTTATCCCCCCCCTAATCCCAGGTGGGAAAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTtaat                                                                                                                                   <  1:252157/228‑1 (MQ=255)
                                                                                                                                              ggtaaaaaaaGTCAC‑GGCTTGGGGCGGAAG‑GT‑TTTTGTCGCATATC‑GCCAAAAA‑AGGGTGTTTTTAACCCGGCTTCAGCTTTATCCCCCCCCGGATCACAAGTGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACg                                                                                                  <  1:503887/246‑1 (MQ=255)
                                                                                                                                              aataaaaaaaGTCCCTGTTTTGT‑GGGAAAA‑GT‑TTTTGGCCCATATCGCCAAATAA‑A‑GGTGTTCTTAATCCGTCTTCAGCTTTCTCCCCCCCCTAATCAAAAGTGGAAAAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACg                                                                                                  <  1:730858/246‑1 (MQ=255)
                                                                                                                                              aGTTAAAAAAGTCACTGACTTGA‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAgg                                                                                                                                                                                            >  1:58757/1‑160 (MQ=255)
                                                                                                                                                                 ttGT‑GCGGAAT‑GT‑TTTTGCCGCATAGCGGCAAATAA‑G‑CGTGTTCTTACCCCGGTTTCAGCTTTATCCCCCCACGGACCACAAGTGGAATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAAc                                                                              <  1:274150/251‑1 (MQ=255)
                                                                                                                                                                  ttgtgggaaaa‑GT‑TTTTGCCACATAGCGGCAAAAAA‑A‑GGTGTTCTTAACCCGGCTTCAGCATTTTCCCCCCCCTGACCACAAGTGGAAAAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAAc                                                                              <  1:632083/243‑1 (MQ=255)
                                                                                                                                                                               gtttttGTCGCATACCGCCAAATAA‑C‑CGTGTTCTTAATCCGTCTTCAGTTTTATCCCCCCACCAACCACAATCCGAATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCt                                                                    <  1:798209/248‑1 (MQ=255)
                                                                                                                                                                                          ggagatgtgtaaaagagaagggTCTTATCCCGGCTTCAGCATTATCCCCCCACGGATCAAAGGGCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCAc                                                           <  1:591327/230‑1 (MQ=255)
                                                                                                                                                                                                                       ccGGCTTCAGCTTTATCCCCCCACTGATCACAAGTGGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCAttgttg                             <  1:153715/251‑1 (MQ=255)
                                                                                                                                                                                                                            aaagagacattccccccattGATCAAAAGTGGAATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTaa                        <  1:59664/232‑1 (MQ=255)
                                                                                                                                                                                                                                gCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTag                                                                                                                                                    >  1:531375/1‑123 (MQ=255)
                                                                                                                                                                                                                                     attCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTa               >  1:591468/1‑251 (MQ=255)
                                                                                                                                                                                                                                     attCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTa               >  1:761021/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ttCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTtaat                                                                                                                                   >  1:77434/1‑134 (MQ=255)
                                                                                                                                                                                                                                      ttCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTa               >  1:714540/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ttCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTa               >  1:378178/1‑250 (MQ=255)
                                                                                                                                                                                                                                             acTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTgccgc                                                                                                                     >  1:888651/1‑141 (MQ=255)
                                                                                                                                                                                                                                              cTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTACCTCGGTaa      >  1:497179/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTCACAAGTCGAATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTACCTCGGTAagt    <  1:202308/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTACCTCGGTAagtag  >  1:551233/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
CGGCAGCAGTTCTGGCTGCAGCATTTGCCATGGTAGGATCAAGGTTGCCACCGTAACTCGGATAAACTTTTCCTTTGGATACATCGATATTTTCC‑GGGCGGAGTAGTGTTCTCAGTAGGGCTTCATCTTGATTATCGACCGAGTTAAAAAAGTCACTGACTTGT‑GCGGAAA‑GG‑TTTTGTCGCATATCGGCCAATAA‑A‑CGTGGTCTTAATCCGGCTTCAGCATTATTCCACCACTGATCACAAGTCGGATAACCGCCATTACCGTTGTTTGGTAACCCTGCGTCACGTGCGTTATCGTAAGGAAATAGTTTAATCGGCTGTTGTGCCCGATAGCGAGGGTAGAGCTCACTATTTTTAATAAGCAGTTTGCCGCCTAACCAACTGATTTCACGAGCTTGCGCTTCATTTAAACTGACCTGGCTTTCTTGCACTTTCCGTCGTGCCGGAACGTAGCATTGTTGCGTAAATTGCTGTACCTCGGTAAGTAG  >  1/1692896‑1693382

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑