Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA46.2% 58.5 / 52.3 78L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (32/10);  new base A (0/36);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.90e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039
                                                                                                                                                       |                                                                                                                                                        
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                         >  2:2026022/1‑150 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  1:1243860/151‑1 (MQ=255)
 aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       <  1:1848082/151‑1 (MQ=255)
  tAAGCGTAGCGCTTCAGGCAATTCCAGCCGCGGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAC‑CCGCGCATTCATCCAACGCCATGTGCAAAAAT‑GCATT‑CGCGTCGGCTGTCTGCCAGCTA‑TTTTCCGCCGcaacaa                                                                                                                                                        <  2:1237731/148‑1 (MQ=14)
  tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                       >  2:768295/1‑150 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                                      <  2:1258916/150‑1 (MQ=255)
   aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                     >  1:1792281/1‑151 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                    >  2:352727/1‑151 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                                 >  1:727587/1‑151 (MQ=255)
       gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgac                                                                                                                                                 >  1:1616099/1‑149 (MQ=255)
          gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                              >  1:1635819/1‑151 (MQ=255)
            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                                            >  1:668219/1‑151 (MQ=255)
             cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                           >  1:1114480/1‑151 (MQ=255)
             cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                           >  1:1683744/1‑151 (MQ=255)
                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                        >  2:1068227/1‑151 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                                  >  1:281082/1‑151 (MQ=255)
                      aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                                   >  2:1875817/1‑150 (MQ=255)
                       tccccgCCGCAGCCCTGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCAACAC‑CCGCGCTTTTACCCAACGCCTTGTGAAAAAAT‑CCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                                    <  1:1319676/143‑1 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGAGGCTAACGTGATTCATTTCCGAAAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAACT‑GCCTT‑CGCGGACGCTGTCTGCCAGCTATTTTTCAGCAGCAACAACACCGCCGTTCACTCCAGTAGTGGc                                                                                                                               >  2:2139920/1‑150 (MQ=255)
                                     cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggg                                                                                                                         >  2:1700009/1‑143 (MQ=255)
                                      tgtgCCGGCGCCTCCCGTGATTATTTCCCCCCAACAC‑CCGCGCTTTTACCCACCGCCTTGTGAAAAAAT‑GCCTC‑CGCGGCGGTGGTCGGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCaa                                                                                                                       <  2:1704582/145‑1 (MQ=255)
                                           cggcgcCTCCCGTCTTTCATTCCCCCCAACAC‑CCCCGCATTTACCCACCCCGATGGGAAAAAGG‑GCCTC‑GGCGCCGGCTCTCTCCCGTTTATTTCCCGGCCGCACAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAg                                                                                                            <  1:1490716/147‑1 (MQ=255)
                                            aGCGGCTACCGTGAATCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGAGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggggaagca                                                                                                                 >  2:1358711/1‑136 (MQ=25)
                                             gcgcCTCCCGTGTTTCATTCCCCCCAACAC‑CCGCCCTTTTACCCACCGCCATGTGAAAAAAG‑GCCTC‑CGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCttt                                                                                                        <  1:36197/151‑1 (MQ=255)
                                             gCGGCTCCCGTGTTTTTTTCCCCCAAACAC‑CCGCGCTTTTACCCAACGCCATGTAAAAAAAT‑CCCTT‑CGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCt                                                                                                          <  2:2131483/149‑1 (MQ=255)
                                               ggCTCCCTTGATTCATTCCCGCCAACAC‑CCGCGTTTTTACCCACCCCATTGTGAAAAAAT‑GCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCt                                                                                                          <  1:420232/147‑1 (MQ=255)
                                               ggCTCCCGTATTTCTTTCCCGCCAACAC‑CCGCGCTTTTATCCACCCCCATGTGAAAAAAT‑GCCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTa                                                                                                       <  1:640800/150‑1 (MQ=255)
                                                gCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggtgccggggaata                                                                                                               >  2:1180879/1‑131 (MQ=255)
                                                 atACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGTAc‑gagggaaagatct                                                                                                         >  2:1299634/2‑137 (MQ=25)
                                                    cggTGATTCATCCCCCCCAACAC‑CCGCGAATTAATCCACCGCATTTAGAAAAATT‑GCCTC‑CGCGGCGGCTGTCCCCCAGCTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGccc                                                                                                  <  1:1910652/148‑1 (MQ=255)
                                                    ccGTGTTTCATTCCCGCCAACAA‑CCGCCATTTTACCCAACCCCTTGTAAAAATTT‑GCCTC‑CGGGGCGCCTGTTTCCCAGCT‑TTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGccc                                                                                                  <  1:2177338/149‑1 (MQ=255)
                                                    ccGTGTTTC‑TTTCCCCCAAAAACCCGCCCTTTTACCCACCGCCATGTCAAAAAAT‑CCCTT‑CGCGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAg                                                                                                     <  1:972508/147‑1 (MQ=255)
                                                      gtttttCTTTCCCCCCACCAC‑CCGCGCTTTTATCCACCCCCTTGTGAAAAAAG‑GCCTT‑GGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGt                                                                                                <  1:482270/146‑1 (MQ=255)
                                                       tgttcctttccccccAACAC‑CCGCCCTTTCCCCCACCGCCTTGTGAAAAAGC‑GCCTG‑CGCGGCGGTTTTCGCCCAGTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGt                                                                                                <  2:1097151/136‑1 (MQ=255)
                                                        gATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAG‑AAAAAACTTTaatacc                                                                                                  >  1:146954/1‑142 (MQ=25)
                                                              tttCCCCCAACAC‑CCGCCCTTTTCACCACCCCCATGTGAAAAAAT‑CCCTT‑CGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                                         <  2:560526/149‑1 (MQ=255)
                                                              ttccccccAACAC‑CCCCCCTTTTACCCACCGCCATGCGAAAAAAT‑CCCTT‑GGCGGCGGCTGTCGGCCAGTTATTTTGCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGt                                                                                                <  2:1653482/136‑1 (MQ=255)
                                                               ttCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAAAAACTTTTAACCCCCCCTGTtttgt                                                                                      >  1:267943/1‑151 (MQ=25)
                                                                 cccccaccaC‑CCGCGCATTTATCCACCGCCTTGTGAAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTtgtgt                                                                                      <  1:464253/146‑1 (MQ=255)
                                                                      aacac‑cCGCCCTTTTACCCACCGCCTTGTGAAAAAAG‑CCCTT‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                                  <  1:1032545/148‑1 (MQ=255)
                                                                         aa‑cCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAAAAACTTTAAACCCGTCATGTTTTGTGGCaa                                                                                 >  1:382756/1‑146 (MQ=255)
                                                                             cgcgCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAA‑GGG‑AAAAGACTTTAAACCC‑TCATGGTGTGTggggttag                                                                              >  2:333337/1‑140 (MQ=25)
                                                                              gccaaTTTACCCAACCCCTTGT‑CAAAAATCCCCTT‑CGGGGCGGTTTTCCGCCGTCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                         <  1:1508304/146‑1 (MQ=255)
                                                                                   tttATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑AAG‑AAAAAAATTTTAAACCCCCCTGTTTTGTGGGGAACGCTGCTgg                                                                       >  1:787669/1‑147 (MQ=21)
                                                                                   tttATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGA‑GAG‑AAAAAAATTTTAAACCCCCATGTTTTGTGGGGAAAGCTGCTgg                                                                       >  1:465384/1‑147 (MQ=25)
                                                                                   tttATCCAACCCCATGTAAAAAAAT‑CCCTT‑CGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                                       <  1:342028/147‑1 (MQ=255)
                                                                                   tttACCCACCGCCATGTCAAAAAAT‑CCCTC‑CGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                        <  1:904170/146‑1 (MQ=255)
                                                                                   tttACCCACCCCCATGTGAAAAAAT‑GCTTC‑CGCGGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                     <  1:1969844/149‑1 (MQ=255)
                                                                                   tttACCCAACGCCATGTGAAAAATT‑CCCTC‑GGGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                         <  1:962394/145‑1 (MQ=255)
                                                                                        ccAACCCCATGTGCAAAAAT‑CCCTC‑CGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                              <  1:1429417/151‑1 (MQ=255)
                                                                                         caccccaagggcaaaaaaT‑CCCTC‑CGGGGGGGCTTTCTCCCAGGTTTTTTCCCCCACCAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                                  <  2:721553/133‑1 (MQ=255)
                                                                                         cAACCCCATGTGAAAAAAG‑TCCTT‑CCGGGCGTGTTGCCCCCATGTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                                  <  2:2177442/146‑1 (MQ=255)
                                                                                           aCGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAG‑AAAAAAATTTAAACCCCGCCTGTTTTGTGGGGAACGCTGCTGGGAaaaaaaggg                                                            >  1:887493/1‑145 (MQ=25)
                                                                                            cccctggtaaaaaaaT‑GCCTC‑CGGGGGGGTTGCCTGCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                                <  2:468096/135‑1 (MQ=255)
                                                                                            cGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAA‑AAG‑AAAAAAAATTATAACCCCCCTCTTGTGTGGGGAACGCTGCTGGTGACaaaatgg                                                            >  1:1242600/1‑144 (MQ=21)
                                                                                              ccttgtaaaaaatT‑GCCTG‑CGGGCGGGCTGCCGGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                             <  2:626541/138‑1 (MQ=255)
                                                                                               cATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAG‑AAAAAAATTTTATCCCCTCCTGTTGTGTGGGGAACGCTGCTGGTGACAAtgtg                                                             >  1:1613105/1‑145 (MQ=37)
                                                                                                    gaaaaaaa‑GCCTTCCGGGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCGTGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                                   <  2:2076678/149‑1 (MQ=255)
                                                                                                       aaaaT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAGAAAAAACCTTTAAGCCC‑CCCTGTTGTTTTGGGATCCCTGCTGGTGAAAAATTGGGAAGGGGAggggg                                               >  1:810743/1‑147 (MQ=255)
                                                                                                            tcccTT‑CGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                                            <  1:1498079/148‑1 (MQ=255)
                                                                                                              ccTC‑CGCGGCGCCTTTCGCCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                         <  2:434068/151‑1 (MQ=255)
                                                                                                               cTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAGAAAAAAGTTTTAACCCCCCCA‑GTGTTGTGGGGATAGCTGCTTGGGAAAAAGTGGGGAGGGAAGAgcgg                                              >  1:2049355/1‑145 (MQ=255)
                                                                                                                   cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAG‑AAAAAAATTTTAAACCCCCCTTTTTTGTGGGGATACCTGCTGGGTACAAAGTGGGAAggggaaggggggggatc                                        >  1:307131/1‑134 (MQ=255)
                                                                                                                      ggcggcTGTCCCCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                                         <  2:1993011/144‑1 (MQ=255)
                                                                                                                       gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACTGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  2:89198/151‑1 (MQ=255)
                                                                                                                            tttcgccgagtttttttcCCGCCAAAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  2:1370889/133‑1 (MQ=255)
                                                                                                                               ctcccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGT‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                              <  2:506444/143‑1 (MQ=255)
                                                                                                                                  ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                          <  1:1550277/147‑1 (MQ=255)
                                                                                                                                  ccAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt                       <  1:932212/150‑1 (MQ=255)
                                                                                                                                  ccAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                                 <  1:1800706/140‑1 (MQ=255)
                                                                                                                                  ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑GGG‑AAAAAAACTTAAACCCCCCCCGTTCTGTGGCGATCGCTCCGGGGGACCAAGTGGGAAGGGGAGGGCGGCGGAACAGCTccccaaacccgc                        >  2:814227/1‑139 (MQ=255)
                                                                                                                                         tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGA‑GAG‑AAAAAACTTTTAACCCCCCCTTTTTTGTGGGGATCGCTGCTGGGAAAAAAGTGGAAAAGGAAGGGGGGGGGAACAGCCCCCCAACCGCGCCAATTTAg                >  1:135170/1‑150 (MQ=255)
                                                                                                                                         ttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                            <  2:1308814/138‑1 (MQ=255)
                                                                                                                                         ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa                 <  1:395976/149‑1 (MQ=255)
                                                                                                                                          tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt               <  1:1516533/150‑1 (MQ=255)
                                                                                                                                             cccgccaAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc             <  2:748908/148‑1 (MQ=255)
                                                                                                                                              ccgccgAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGc           <  2:281047/150‑1 (MQ=255)
                                                                                                                                              cccccacaaaaaACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc             <  1:1950260/142‑1 (MQ=255)
                                                                                                                                              ccacggcAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                          <  2:1432384/130‑1 (MQ=255)
                                                                                                                                                  cgcAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAA‑GAGAAAAAAATTTTAAACCCG‑CCTGTTGTGTGGGAATTGCTGCTGGTAAAAAATTGGGAAGGGGGGGGGGGGGGATTCACTCCCGAACCCCGCTAATTTAGTTTGCCt         >  1:1069530/1‑148 (MQ=255)
                                                                                                                                                  cgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt         <  1:1426215/148‑1 (MQ=255)
                                                                                                                                                    caacaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGGAA‑AAAAAAATTTAACCCCTCCTTGTT‑TTTGGCGCACGCCGCTGGTGAAAATGTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGGAACCCGCCAAATGAGGTCGCCTccc      >  2:1498248/1‑147 (MQ=255)
                                                                                                                                                      acaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac  <  1:1056026/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                        
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGAC  >  NC_000913/357742‑358039

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