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.

GCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCTTG  >  NC_000913/4257453‑4257913
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           
gcgcGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATTGTGACTTGCTGGCAGTGCATAAAACTCATGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGTCGAAGGTACCGTTAAGCGCCTGAAAAAACATGGCAATAAAGTCGAACTGTTGCCATAAAATAGCGAGTTTAAACCAATTGTCGTTTACCTTCGTTAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAt                                                                                                                                                                                                                          <  2:272667/249‑1 (MQ=255)
            gATGGCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGGCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAGTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAtt                                                                                                                                                                                                                         <  1:225629/238‑1 (MQ=255)
                    tGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGCGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGGCTGGCTGtt                                                                                                                                                                                                    >  1:27006/1‑249 (MQ=255)
                      aaaGATGTCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGACGTTACCGTTCAGCGCCTGAAACAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGATTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGTGTTTTCGCAACGGCGAACGGCTGTaatc                                                                                                                                                                                                 >  2:353185/1‑249 (MQ=255)
                              cGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGATGGGGGTATTTGCAACGGCGACCGGCTGttaactatttact                                                                                                                                                                                         >  1:37012/1‑239 (MQ=255)
                                gCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTATTTGCAAAGGGGACTTGCTGTAACATATCTCTg                                                                                                                                                                                         >  1:276591/1‑249 (MQ=255)
                                 cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAAACAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTGATTTGCAACGGCGGCAGGCTGTAAAATATTTCTgtg                                                                                                                                                                                       >  1:104830/1‑248 (MQ=255)
                                         atggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGGCTGGCTGGAACATATTTCTGAGACcgcggtg                                                                                                                                                                              >  2:228615/1‑248 (MQ=255)
                                           ggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTAACATATATCTTGGACCGCGATgcc                                                                                                                                                                            >  2:9388/1‑251 (MQ=255)
                                                   gCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCg                                                                                                                                                                    <  1:275583/251‑1 (MQ=255)
                                                   gCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGGTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCTCAACGGCGACTGGCTGGAACATATCTCGGAGACCGCGATGCCGCCTGGc                                                                                                                                                                     <  2:27006/250‑1 (MQ=255)
                                                                              aCGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTtca                                                                                                                                         <  2:91780/251‑1 (MQ=255)
                                                                                        cAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGt                                                                                                                              <  1:228615/251‑1 (MQ=255)
                                                                                                             gatgaCGAAGTTTCCGTTAAGCGCCTGAAAAAACAGGGCAATGAAGTCGAACTGTGGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGGTTATGCGCAACGGCGACTGGCTGTAACATATCTCTGAGTCCGCGATGCCGCCTGGCGTCGCGGTTTGTTGTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACtt                                                                                                          <  2:432122/250‑1 (MQ=255)
                                                                                                                   gAAGTTACCGTTAAGCGCCTGACAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGCCTGGGTGTAACCTATCTCTGGGGACGCGATGCCGCCATGCGGCGCGGGTTGGTTTTCTTCTCTCTTCATCAGGTCTGTCTGCATGGGATTCCTCACtttctctg                                                                                                    >  2:213950/1‑244 (MQ=255)
                                                                                                                        tACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAAAGGGGACTGGCTGTAACATATCTCTGGGAACGCGATACCGCCTGGCGTCGCGGGTTGTTTTTTAACTCTCTTCATCAGGCTTGTCTGCAAGGAATTCCTCACTTCATATGATaaa                                                                                               >  1:482627/1‑250 (MQ=255)
                                                                                                                          ccGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGTCAGAAAATAGCGAGTTTAAACCAATTGTCGGTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCGCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCAGTCCTCACTTCATCTGATAAAGCa                                                                                            <  2:95701/251‑1 (MQ=255)
                                                                                                                           cGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGTTATTCGCAACGGCGACTGGCTGTAAACTATCTCTGGGGCCGCGATGCCGCCTGGCGTCGCGGGTTGTTTTTCATTTCTCTTCATCAGGCTTTTCTTCCTGGCCTTCCTCACTTTCTCTGATAAAgac                                                                                            >  2:7746/1‑248 (MQ=255)
                                                                                                                                           aaaCAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAGTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCtt                                                                          <  2:340759/251‑1 (MQ=255)
                                                                                                                                                gggCAATCAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCCTTAcccc                                                                     >  1:342768/1‑250 (MQ=255)
                                                                                                                                                  gCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATTTCTCTGTGACGGCGCGGACGCCTGGCGTCGCGGTTTGTTTTTTTTCTCTCTTCCTCAGGGCTTTCTGCATAGGATTCCTTACTTCATCTTATAAAGCACTCTGTCATCTCGCCTTTCCCATg                                                                   >  2:356032/1‑251 (MQ=255)
                                                                                                                                                       aaaGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGAtttt                                                              <  1:9388/251‑1 (MQ=255)
                                                                                                                                                        aaGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCCGCAGAGCTTCACCATTGAAGAGCTGGCGGTTGGGGTTATTCGCAACGGCGAATGGCAGTACAATATCTCTGAGACCGCGCTGCCGCCTGGCGTCGCGGGTTGTTTTTCCTCTCTCTTCACCAGGCTTGTCTGCAAGGCCTTCATCCCTTTATCTGATAAAGCACACTTGCATCTCGCCTTACCCCTGAttttt                                                             >  1:305185/1‑250 (MQ=255)
                                                                                                                                                              aaCTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTATTCGCAACGGCGCCTGGCTGTAACATATCTCTGAGGCCGCGGTTCCGCCCGGCGGCGCGGGTTTTTTTTTTTCTCTCTTCTTCAGGGTTTTCTGCATGGCATTTCTCACTTCATCTTTTAAAGCACTCTGGCAACTCCCCTTTCCCCTGCTTTTTTcccat                                                       >  2:157251/1‑248 (MQ=255)
                                                                                                                                                                  gTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCa                                                  <  1:100881/251‑1 (MQ=255)
                                                                                                                                                                   ttGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGATTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGGCAAGGCGAATGGCTGCAACAAATCTCTGGGACCGCGAAGACGCCTGGAGGCGCGGGTTGGTTTTTATCTCTATTCATCAGGCTTGTCTGGCTGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAc                                                 >  1:190281/1‑251 (MQ=255)
                                                                                                                                                                     gCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAcc                                                <  1:258017/250‑1 (MQ=255)
                                                                                                                                                                       cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAccgtt                                             <  1:157251/251‑1 (MQ=255)
                                                                                                                                                                       cAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGCTGTAACATATCTCTGAGACCGCGGTGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCAccgttc                                            >  1:140377/1‑251 (MQ=255)
                                                                                                                                                                               aGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCt                                     <  1:7746/251‑1 (MQ=255)
                                                                                                                                                                                    gTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGGAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTTAGACCGCGATTCCGCCTGGCGTCGCGGTTTGGTTTTCATCTCTCTTCATCAGTCTTGGCTGCGTGGCATTCCTCACTTCATCTTATAATGCACTCTGGCTTCGCGCCTTACCCATGATTTTCTTCAAGATCTCCGTTCCTGTGCTGGGAc                                <  2:314232/251‑1 (MQ=255)
                                                                                                                                                                                       tAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTTACATATCTCTGGGACCGCGATGCCGCCCGGCGTCGCGGTTTGTTTTTCCTCTCTCTTCATCCGGCTTGTCTGCCTTGCATTCCTCCCTTCCTCTGATAAAGCACCCTGGCATCTCGCCTTTCCCCTGCTTTTCTCCAATATCACCGTTCCGTTGCTgg                                   >  1:212729/1‑245 (MQ=255)
                                                                                                                                                                                       tAAACCAATTGTCGGTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTgg                                   <  2:212729/245‑1 (MQ=255)
                                                                                                                                                                                        aaaCCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGGCTggg                            >  1:249554/1‑250 (MQ=255)
                                                                                                                                                                                           ccAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGGTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCCATATCACCGTTCCGTTGCTGGGACTGGTCg                          >  1:261253/1‑250 (MQ=255)
                                                                                                                                                                                              aTTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCCACTGGCTGGAAAATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTTTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCCTTACCCATGATTTTCTTCAATATCACCGTTCCGTTGCTGGGGCTGGTCGATAc                      >  1:406849/1‑251 (MQ=255)
                                                                                                                                                                                                      tGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGGTGCCGCCTGGCGTCGCGGGTTGTTTTTCATCTCTCTTCCTCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCCTGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGGCGGTAcggcgggg               >  1:418033/1‑248 (MQ=255)
                                                                                                                                                                                                            tCGTCAGCAGAGCGGCACCATTGATGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCTCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCGCACGGCATCTGATAAAGCACTCTGTCATCTCGGCTTACCCATGATTTTGTTCAATGTCGCCGGTCCGTGGCTGGGACTGGGCGATAGTGCGGTAATTGGt         <  2:231619/250‑1 (MQ=255)
                                                                                                                                                                                                              gTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCCACGGCCCCTCGCCGTAACATATCCCTGAGACCGCGATGCCCCCCCGCGGCGCGGGTTTTTTTTCTTTTCTCCTCATCCGGGCTGTCTGCAGTGCCTTCCCCCCCTCTTCTTATAAAGCCCTCCGGCCTCCCCCCCTTACCACGCTTTTTTCCAAATACACCGTTCCGTTGCTGGGGCCGGTCCATTcggcgggtatttgtgct      >  2:196582/1‑240 (MQ=255)
                                                                                                                                                                                                                 agcagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTTGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGGCGCGGTTTGTTTTTCATCTTTCTTCATCCGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAACACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGGAATTAGTGATCtt   >  1:195864/1‑251 (MQ=255)
                                                                                                                                                                                                                 agcagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCtt   <  2:261253/251‑1 (MQ=255)
                                                                                                                                                                                                                  gcagAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAAAGGCGACTGGCTGTAAAATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCAATCTGGCATATCGCCTTACCCCTGATTTTCTCCAATATCAACGTTCCGTTGCTGGGACTGGGCGATACGGCGGGAAttttgtgtcttt  >  1:107015/1‑241 (MQ=255)
                                                                                                                                                                                                                    agagCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTtc                                                             <  1:446401/191‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                           
GCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCCTGGCGTCGCGGTTTGTTTTTCATCTCTCTTCATCAGGCTTGTCTGCATGGCATTCCTCACTTCATCTGATAAAGCACTCTGGCATCTCGCCTTACCCATGATTTTCTCCAATATCACCGTTCCGTTGCTGGGACTGGTCGATACGGCGGTAATTGGTCATCTTG  >  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: