Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,257,6980TG32.6% 54.0 / 10.3 43V195G (GTT→GGT) lexADNA‑binding transcriptional repressor LexA
Reads supporting (aligned to +/- strand):  ref base T (11/18);  new base G (13/1);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 8.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCTTG  >  NC_000913/4257453‑4257913
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           
gcgcGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATTGTGACTTGCTGGCAGTGCATAAAACTCATGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGTCGAAGGTACCGTTAAGCGCCTGAAAAAACATGGCAATAAAGTCGAACTGTTGCCATAAAATAGCGAGTTTAAACCAATTGTCGTTTACCTTCGTTAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTAt                                                                                                                                                                                                                          <  2:272667/249‑1 (MQ=255)
            gATGGCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGGCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAGTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTAtt                                                                                                                                                                                                                         <  1:225629/238‑1 (MQ=255)
                    tGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGCGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGGCTGGCTGtt                                                                                                                                                                                                    >  1:27006/1‑249 (MQ=255)
                      aaaGATGTCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGACGTTACCGTTCAGCGCCTGAAACAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGATTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGTGTT‑TTCG‑CAACGGCGAACGGCTGTaatc                                                                                                                                                                                                 >  2:353185/1‑249 (MQ=255)
                              cGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGATGGG‑GGTATTTG‑CAACGGCGACCGGCTGttaactatttact                                                                                                                                                                                         >  1:37012/1‑239 (MQ=255)
                                gCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGG‑GGTATTTG‑CAAAGGGGACTTGCTGTAACATATCTCTg                                                                                                                                                                                         >  1:276591/1‑249 (MQ=255)
                                 cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAAACAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTGATTTG‑CAACGGCGGCAGGCTGTAAAATATTTCTgtg                                                                                                                                                                                       >  1:104830/1‑248 (MQ=255)
                                         atggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGGCTGGCTGGAACATATTTCTGAGACcgcggtg                                                                                                                                                                              >  2:228615/1‑248 (MQ=255)
                                           ggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGACTGGCTGTAACATATATCTTGGACCGCGATgcc                                                                                                                                                                            >  2:9388/1‑251 (MQ=255)
                                                   gCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCg                                                                                                                                                                    <  1:275583/251‑1 (MQ=255)
                                                   gCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGGTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCT‑CAACGGCGACTGGCTGGAACATATCTCGGAGACCGCGATGCCGCCTGGc                                                                                                                                                                     <  2:27006/250‑1 (MQ=255)
                                                                              aCGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtca                                                                                                                                         <  2:91780/251‑1 (MQ=255)
                                                                                        cAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGt                                                                                                                              <  1:228615/251‑1 (MQ=255)
                                                                                                             gatgaCGAAGTTTCCGTTAAGCGCCTGAAAAAACAGGGCAATGAAGTCGAACTGTGGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGG‑GTTATGCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGTCCGCGATGCCGCCTGGCGTCGCGGTTTGTTGTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACtt                                                                                                          <  2:432122/250‑1 (MQ=255)
                                                                                                                   gAAGTTACCGTTAAGCGCCTGACAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGCCTGGGTGTAACCTATCTCTGGGGACGCGATGCCGCCATGCGGCGCGGGTTGGTTTTCTTCTCTCTTCATCAGG‑TCTGTCTGCATGGGATTCCTCACtttctctg                                                                                                    >  2:213950/1‑244 (MQ=255)
                                                                                                                        tACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAAAGGGGACTGGCTGTAACATATCTCTGGGAACGCGATACCGCCTGGCGTCGCGGGTTGTTTTTTAACTCTCTTCATCAGG‑CTTGTCTGCAAGGAATTCCTCACTTCATATGATaaa                                                                                               >  1:482627/1‑250 (MQ=255)
                                                                                                                          ccGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGTCAGAAAATAGCGAGTTTAAACCAATTGTCGGTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCGCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCAGTCCTCACTTCATCTGATAAAGCa                                                                                            <  2:95701/251‑1 (MQ=255)
                                                                                                                           cGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGG‑GTTATTCG‑CAACGGCGACTGGCTGTAAACTATCTCTGGGGCCGCGATGCCGCCTGGCGTCGCGGGTTGTTTTTCATTTCTCTTCATCAGG‑CTTTTCTTCCTGGCCTTCCTCACTTTCTCTGATAAAgac                                                                                            >  2:7746/1‑248 (MQ=255)
                                                                                                                                           aaaCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAGTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑tt                                                                          <  2:340759/251‑1 (MQ=255)
                                                                                                                                                gggCAATCAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCC‑TTAcccc                                                                     >  1:342768/1‑250 (MQ=255)
                                                                                                                                                  gCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATTTCTCTGTGACGGCGCGGACGCCTGGCGTCGCGGTTTGTTTTTTTTCTCTCTTCCTCAGGGCTT‑TCTGCATAGGATTCCTTACTTCATCTTATAAAGCACTCTGTCATCTCGCC‑TTTCCCATg                                                                   >  2:356032/1‑251 (MQ=255)
                                                                                                                                                       aaaGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGAtttt                                                              <  1:9388/251‑1 (MQ=255)
                                                                                                                                                        aaGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCCGCAGAGCTTCACCATTGAAGAGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGAATGGCAGTACAATATCTCTGAGACCGCGCTGCCGCCTGGCGTCGCGGGTTGTTTTTCCTCTCTCTTCACCAGG‑CTTGTCTGCAAGGCCTTCATCCCTTTATCTGATAAAGCACACTTGCATCTCGCC‑TTACCCCTGAttttt                                                             >  1:305185/1‑250 (MQ=255)
                                                                                                                                                              aaCTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGG‑GGTATTCG‑CAACGGCGCCTGGCTGTAACATATCTCTGAGGCCGCGGTTCCGCCCGGCGGCGCGGGTTTTTTTTTTTCTCTCTTCTTCAGG‑GTTTTCTGCATGGCATTTCTCACTTCATCTTTTAAAGCACTCTGGCAACTCCCC‑TTTCCCCTGCTTTTTTcccat                                                       >  2:157251/1‑248 (MQ=255)
                                                                                                                                                                  gTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑Ca                                                  <  1:100881/251‑1 (MQ=255)
                                                                                                                                                                   ttGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGATTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCGGCAA‑GGCGAATGGCTGCAACAAATCTCTGGGACCGCGAAGACGCCTGGAGGCGCGGGTTGGTTTTTATCTCTATTCATCAGG‑CTTGTCTGGCTGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CAc                                                 >  1:190281/1‑251 (MQ=255)
                                                                                                                                                                     gCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CAcc                                                <  1:258017/250‑1 (MQ=255)
                                                                                                                                                                       cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CAccgtt                                             <  1:157251/251‑1 (MQ=255)
                                                                                                                                                                       cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGACTG‑CTGTAACATATCTCTGAGACCGCGGTGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CAccgttc                                            >  1:140377/1‑251 (MQ=255)
                                                                                                                                                                               aGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCt                                     <  1:7746/251‑1 (MQ=255)
                                                                                                                                                                                    gTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGGAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTTAGACCGCGATTCCGCCTGGCGTCGCGGTTTGGTTTTCATCTCTCTTCATCAGT‑CTTGGCTGCGTGGCATTCCTCACTTCATCTTATAATGCACTCTGGCTTCGCGCC‑TTACCCATGATTTTCTTCAAGAT‑CTCCGTTCCTGTGCTGGGAc                                <  2:314232/251‑1 (MQ=255)
                                                                                                                                                                                       tAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGACTGGCTGTTACATATCTCTGGGACCGCGATGCCGCCCGGCGTCGCGGTTTGTTTTTCCTCTCTCTTCATCCGG‑CTTGTCTGCCTTGCATTCCTCCCTTCCTCTGATAAAGCACCCTGGCATCTCGCC‑TTTCCCCTGCTTTTCTCCAATAT‑CACCGTTCCGTTGCTgg                                   >  1:212729/1‑245 (MQ=255)
                                                                                                                                                                                       tAAACCAATTGTCGGTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTgg                                   <  2:212729/245‑1 (MQ=255)
                                                                                                                                                                                        aaaCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGGCTggg                            >  1:249554/1‑250 (MQ=255)
                                                                                                                                                                                           ccAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGGTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCCATAT‑CACCGTTCCGTTGCTGGGACTGGTCg                          >  1:261253/1‑250 (MQ=255)
                                                                                                                                                                                              aTTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCCACTGGCTGGAAAATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTTTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCC‑TTACCCATGATTTTCTTCAATAT‑CACCGTTCCGTTGCTGGGGCTGGTCGATAc                      >  1:406849/1‑251 (MQ=255)
                                                                                                                                                                                                      tGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGGTGCCGCCTGGCGTCGCGGGTTGTTTTTCATCTCTCTTCCTCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCCTGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGACTGGGCGGTAcggcgggg               >  1:418033/1‑248 (MQ=255)
                                                                                                                                                                                                            tCGTCAGCAGAGCGGCACCATTGATGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCTCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCGCACGGCATCTGATAAAGCACTCTGTCATCTCGGC‑TTACCCATGATTTTGTTCAATGT‑CGCCGGTCCGTGGCTGGGACTGGGCGATAGTGCGGTAATTGGt         <  2:231619/250‑1 (MQ=255)
                                                                                                                                                                                                              gTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CCACGGCCCCTCGCCGTAACATATCCCTGAGACCGCGATGCCCCCCCGCGGCGCGGGTTTTTTTTCTTTTCTCCTCATCCGG‑GCTGTCTGCAGTGCCTTCCCCCCCTCTTCTTATAAAGCCCTCCGGCCTCCCCCCCTTACC‑ACGCTTTTTTCCAA‑ATACACCGTTCCGTTGCTGGGGCCGGTCCATTcggcgggtatttgtgct      >  2:196582/1‑240 (MQ=255)
                                                                                                                                                                                                                 agcagAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTTG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGTTTGTTTTTCATCTTTCTTCATCCGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGGAATTAGTGATCtt   >  1:195864/1‑251 (MQ=255)
                                                                                                                                                                                                                 agcagAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCtt   <  2:261253/251‑1 (MQ=255)
                                                                                                                                                                                                                  gcagAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GGTATTCG‑CAAAGGCGACTGGCTGTAAAATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCAATCTGGCATATCGCC‑TTACCCCTGATTTTCTCCAATAT‑CAACGTTCCGTTGCTGGGACTGGGCGATACGGCGGGAAttttgtgtcttt  >  1:107015/1‑241 (MQ=255)
                                                                                                                                                                                                                    agagCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTtc                                                             <  1:446401/191‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           
GCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGG‑GTTATTCG‑CAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGG‑CTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCC‑TTACCCATGATTTTCTCCAATAT‑CACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCTTG  >  NC_000913/4257453‑4257913

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