Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT42.1% 81.4 / 23.5 58N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (20/13);  new base T (0/24);  total (20/38)
Fisher's exact test for biased strand distribution p-value = 5.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NC_000913/357725‑358006
                                                                                                                                                         |                                                                                                                                      
aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTa                                                                                                                                        <  1:739448/151‑1 (MQ=255)
     aCTCGTATGTCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCGGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACTCCATGTTCAAAAATGCCTT‑C‑GCGACGGCTGTCTGCCAGCTAttgtt                                                                                                                                   <  2:793562/151‑1 (MQ=14)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                   >  1:957946/1‑151 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                <  1:1116491/150‑1 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                <  2:177564/150‑1 (MQ=255)
        cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                 <  2:808279/149‑1 (MQ=255)
         gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgc                                                                                                                              >  1:1070631/1‑151 (MQ=255)
            ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                           <  2:172441/151‑1 (MQ=255)
             gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                           >  1:1172195/1‑150 (MQ=255)
                 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                        <  1:1232236/149‑1 (MQ=255)
                  aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaa                                                                                                                    <  1:1301364/151‑1 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaaa                                                                                                                   <  1:1205811/150‑1 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAc                                                                                                                  >  1:1303713/1‑151 (MQ=255)
                       cGTAGCGCATCAGGCAATTCCACCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAcc                                                                                                                 <  2:1219159/149‑1 (MQ=255)
                             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtc                                                                                                          >  2:554059/1‑150 (MQ=255)
                                tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑TCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtctctc                                                                                                      >  1:1057746/1‑151 (MQ=255)
                                  aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                     >  1:1185853/1‑150 (MQ=255)
                                        ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCCCTCCAGGAGTg                                                                                              >  2:674692/1‑151 (MQ=255)
                                          ccACCCGCAGCCTTGTGTCGGCGGCTCCCGTGTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCAACCGCCTGGCGAAAAAATCCCTC‑C‑GGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCCCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                            <  1:1324225/151‑1 (MQ=255)
                                           cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                          >  2:715539/1‑151 (MQ=255)
                                              ccGCAGCCTTGTCTCGGCGCCTCCTGTTATTTATTCCCGCAAACAC‑CCGCGCATTAATCCAACGCCTTGTG‑AAAAATGTCTTCC‑GCGGCGCCTGTCTGCAATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCagtag                                                                                                <  2:171586/143‑1 (MQ=255)
                                               cccagccatgtgCCGGCGGCTGCCGTGATCCATCCCCCCCAACAC‑CCGCGCTTTTACCCACCCCGATGTAAAAAAGTTCGTC‑G‑GCGGCGCCTGCCGCCCCGCTATTTTC‑CCCGCCAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                        <  2:1456042/141‑1 (MQ=255)
                                               cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                       >  1:547115/1‑150 (MQ=255)
                                                   gcccTGTGTCAGCGCCTACCGTTTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTT‑C‑GCGGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                   <  1:657477/148‑1 (MQ=255)
                                                      cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                       >  1:585558/1‑143 (MQ=255)
                                                         tttgggCGTCTACCGTGATTCTTCCCCCCCAACAC‑CCGCGCTTTTATCCAGCCCCATGGGAAAAAGTCCCTC‑C‑GCGGCGGCTGTGTCCCAGTTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaaa                                                                              <  1:664317/144‑1 (MQ=255)
                                                         tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggagggggaa                                                                             >  2:692297/1‑140 (MQ=255)
                                                              gcgcCTCCCGTGTTTCATTCCCCCCACCAC‑CCGCGCTTTTCCCCACCCCTTTGTAAAAAAAT‑CCTT‑CGGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg                                                                           <  2:98300/147‑1 (MQ=255)
                                                                ggCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGaaaaaaattt                                                                       >  1:858129/1‑142 (MQ=255)
                                                                ggCTACCGTAATTCTTTCCCGCCAACAC‑CCGGGCTTTTACCCACCCCCCTGTAAAAAAACGCCTG‑C‑GCGGGGGCTGCCTCCCAGTTTTTTTC‑CCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAt                                                                             <  1:15446/143‑1 (MQ=255)
                                                                gcctcccGTGTTTCTTTCCCGCCACCAC‑CCCCGCTTTTACCAACCGCCAGGT‑CAAAAAAGCCTTCC‑GCGGCGGTTGCCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaa                                                                               <  2:230456/136‑1 (MQ=255)
                                                                gcctcccGTGATTAATTCCCGCAACCAC‑CCGCGCTTTTACCCACCCCCATGTGAAAAAAGCCCTT‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCt                                                                         <  1:949156/142‑1 (MQ=255)
                                                                 gCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGgaaaaaattt                                                                       >  1:181921/1‑141 (MQ=255)
                                                                  cccccGTGTTTCTTTCCCGCCAACAC‑CCCCGATTTTACCCACCGCTATGTCAAAAAATCCTTC‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCACAAA‑AAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAg                                                                    <  2:1217260/147‑1 (MQ=255)
                                                                  cTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGAAGaaaaaaactt                                                                       >  2:1323704/1‑140 (MQ=25)
                                                                    ccggttaTTAATTCCCGCCAACAC‑CCCCCCATTTACCCACCCCATTGTGAAAAAATCCCTC‑G‑GCGGCGGCTTTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                 <  2:1346563/145‑1 (MQ=255)
                                                                    aCCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAAAAACTTTAAAccc                                                                 >  1:598663/1‑151 (MQ=255)
                                                                      cggtttttaTTTCCCCCAACAACCCGCGCATTTACCCACCCCCTTGTAAAAAATTCCCTC‑C‑GGGGCGGTTGTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCtt                                                                        <  2:1254825/135‑1 (MQ=255)
                                                                               ttccccccAAAAC‑CCGCCCTTTTACCCACCCCAATGTAAAAAAATCCCTG‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                          <  1:592878/141‑1 (MQ=255)
                                                                               tccccGCCACCAC‑CCGCTTTTTTACCCCACGCCATGTAAAAAAATCCTTC‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCACACA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGACCGGCAtgttg                                                        <  2:608206/146‑1 (MQ=255)
                                                                                  cccccAACAC‑CCGCGCTTTTACCCACCGCCTTGTGAAAAATTGCTTT‑G‑GGGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                     <  1:252997/146‑1 (MQ=255)
                                                                                   ccccaaccacccccgcTTTTACCCACCCCCATGGGAAAAATTCCCTG‑C‑GGGGCGGTTGTC‑GCCAATTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                          <  2:954952/131‑1 (MQ=255)
                                                                                       aacac‑cCGCGCATTTACCCACCCCCTTGTGAAAAAATTCCTT‑C‑GGGGCGGCTTTCTGCCATTTTTTTTC‑CCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                               <  1:1368108/150‑1 (MQ=255)
                                                                                             ccgcacttttccccACCGCCATGTAAAAAAGTGCATG‑G‑GCGGGGGCGGTCGGCGGTTTTTTTGC‑CCGCCACAAA‑AAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                              <  2:688549/134‑1 (MQ=255)
                                                                                               gcgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAACTTTAA‑CCCCTCCTGTTTTGTTGGAATccc                                            >  2:1309730/1‑144 (MQ=25)
                                                                                                       aTCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAACCCCTCCATTTGTGTTGCGAACGCTGCTggggaaaa                                >  2:774083/1‑144 (MQ=25)
                                                                                                         ccaccgccagggaaaaaaaTGCCTG‑C‑GCGGCGCCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                  <  2:1276166/132‑1 (MQ=255)
                                                                                                             ccccttgtaaaaaaaTGCTTC‑C‑GCGGCGGTTGCCGGCCAGTTTTTTTC‑CCCCCGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                              <  2:852803/136‑1 (MQ=255)
                                                                                                               ccagggaaaaaaaTGCCTC‑C‑GGGGCGGTTTTCGGCCAGTTTTTTTC‑CCCCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                          <  2:751548/140‑1 (MQ=255)
                                                                                                                   gtcaaaaaaTCCCTG‑C‑GGGGCGGTTGTCTCCCATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                      <  1:62683/144‑1 (MQ=255)
                                                                                                                   gtaaaaaaaGGCCTC‑C‑GCGGGGCTTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                      <  2:75230/144‑1 (MQ=255)
                                                                                                                    taaaaaaaTCCTTC‑C‑GCGGCGCTTGTCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                     <  2:476208/145‑1 (MQ=255)
                                                                                                                    taaaaaaaTCCCTC‑C‑GCGGCGGCTTTCTCCCATTTTTTTTT‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATAGCTGCTGGTAACAATGTGGAAAGGGa                    <  1:97284/146‑1 (MQ=255)
                                                                                                                    aacAAAAATGCCTC‑C‑GCGGCGGTTGTCTGCCGGTTATTTTC‑CCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                  <  2:1143154/149‑1 (MQ=255)
                                                                                                                         aaaTGCCTG‑C‑GGGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg             <  2:588244/151‑1 (MQ=255)
                                                                                                                          aaTGCCTT‑CGGGGGCGCCTTTCTCCCGTCT‑TTTTT‑CCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg             <  2:1513976/151‑1 (MQ=255)
                                                                                                                               cTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATTGCTGCTGGGAACAAAGGGGGAAGGGGAgggggggg           >  2:1343207/1‑140 (MQ=255)
                                                                                                                                       ggcggcTGTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc  <  1:437978/150‑1 (MQ=255)
                                                                                                                                       ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGACGAGGAAAAGAATATTACCCGACATATTGTGTGGGGCTCGGCGCTGGGGACAATATTGAAAAGGGAGGGCGGGGAGTCAGc    >  2:166503/1‑148 (MQ=255)
                                                                                                                                                         |                                                                                                                                      
ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTC  >  NC_000913/357725‑358006

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: ‑