Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,383,7220GA27.3% 75.2 / 10.1 44F80F (TTC→TTT) ycjWDNA‑binding transcriptional repressor YcjW
Reads supporting (aligned to +/- strand):  ref base G (28/4);  new base A (1/11);  total (29/15)
Fisher's exact test for biased strand distribution p-value = 1.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGT  >  NC_000913/1383475‑1383965
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
cAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGaa                                                                                                                                                                                                                                                     >  2:381999/1‑250 (MQ=255)
          aTGGACATATTTACCTTCAATCTTGCCAATCACCACCACAGGAATATCATATTTATCGAGTTGGGCAAAAAATGTCTCATCAGCTGGCGAACCGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCGTTCTCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGGCGAAACCGA‑CTTCTTCAGCGTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACCGATTA‑c                                                                                                                                                                                                                                         <  2:232020/251‑1 (MQ=255)
             ggCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑Ccc‑g                                                                                                                                                                                                                                     >  1:584197/1‑251 (MQ=255)
               cATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑Ccc‑gcc                                                                                                                                                                                                                                   <  1:218474/251‑1 (MQ=255)
                                      aTCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGGAGTTTTGGTCGAACGGGtcg                                                                                                                                                                                                             >  2:304728/1‑248 (MQ=255)
                                                  ggAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAAt                                                                                                                                                                                                >  1:525177/1‑251 (MQ=255)
                                                     atatCATATTTATCGATTTGGGAAAAAAATGACTCACCCGCTGGCGAACGGAGCATAATAATGCCTTAAAT‑CATTTTCTGTTTAATTTTGCTTTCGAATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACTTTTGTAATATCCCGTAAAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGAAAAACTT‑CAGAAAAAAACGAATTC‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGc                                                                                                                                                                                             <  1:381999/251‑1 (MQ=255)
                                                       atCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGTGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCAt                                                                                                                                                                                           >  1:619028/1‑251 (MQ=255)
                                                          atatTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCa                                                                                                                                                                                         >  1:324138/1‑250 (MQ=255)
                                                                        gggcacacaaTGACTAATCCGCGGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑CGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGGCGAACGGGTAGAAATCACCATAATGGAAACATAACCCGGAGAGGTc                                                                                                                                                                          >  2:255113/9‑251 (MQ=255)
                                                                           cAAAAAATGACACATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgc                                                                                                                                                                       >  2:514528/1‑251 (MQ=255)
                                                                           aaaaaaaTGACTCTTCCTTGGGCGATCTGAGAAAAATATTGC‑TTTAATTCTTTTTCTGTTTATTTTTGTTTTTGCATTTTT‑GCATATCTTCTGCCGATTT‑GTGCGACTTTTGTAATATCACGTAAAAACCTT‑CTTCTTCA‑CTTTGGGGGGTGTTGGCAGGAAAAATTA‑GAGAAAAAAAGGTATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgc                                                                                                                                                                       <  2:558096/250‑1 (MQ=255)
                                                                           aaaaaaaTGACTCACCCGCTGGCGAATAGAGAATAATAATGCTTTTAAT‑AATTTTCTTTTTAATTTTTTTTTCGCTTTTTT‑GCAAATCTTCTGCCGGTTT‑GTGCGATGTTTGTAATCTCCCTTCAAAACCTT‑CTTCTCCAGCTTTGGCTG‑TGATGGCATG‑AAAACTTTCAGAAAAAAACGTATCC‑ACC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgc                                                                                                                                                                       <  2:558090/250‑1 (MQ=255)
                                                                                  tGACTCATCAGCTGGCGAACTGAGCCTAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTTGTCGAACGGGTAGAAATCACCATAATGGCATCAAACCCCGAAGAGGTCAGCGCGCggg                                                                                                                                                                >  2:186748/1‑251 (MQ=255)
                                                                                   gccccATCCGCGGGCGAACTGAGCATAATAATGCTTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCCCATTTTT‑GCAGATCTTCTGCGGGTTT‑GTGCGACGTTTGTAATATCACTTCGAACCCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGAATGAAAAACTT‑CAGAAAAAAACGTATTC‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGc                                                                                                                                                               <  1:514528/247‑1 (MQ=255)
                                                                                   gACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGc                                                                                                                                                               >  1:24011/1‑251 (MQ=255)
                                                                                       cacccGCTGGCGAACTGAGCATAATAATTCTTTTAAT‑CATTTCCTGCTTATTTTTGTTTTCGAATTTTT‑GAAGATCTTCTGCCGGTTT‑GTGCGACGTTTGTAATATCACGTCGAACCCTT‑CTTCTTCAGCTTTGGGGG‑TGATGGCATGCAAAACTA‑CAAAAAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGc                                                                                                                                                           <  1:635249/248‑1 (MQ=255)
                                                                                               ggCGAATTGAGCATAAAAATGCTTTAAAT‑AATTTTCTGCTTAATTTTGTTTCCGCTTTTTT‑GAAGATCTTCGCCCGGATT‑GTGGGCCTTTTGTAATATCCCGGCGAAACCTT‑CTTTTTCAGTTTGGGCGG‑GGAGGGCATGAAAAACTT‑CAGAAAAAAACGATTTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTg                                                                                                                                                    <  1:381040/250‑1 (MQ=255)
                                                                                                              ataatGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCCGCTTGTTTGGCTGATACTGTAAttc                                                                                                                                    >  2:224943/1‑251 (MQ=255)
                                                                                                                  tGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAACACCCGAAGAGGTCAGCGAGCGGGACAGCTTGTTTGGCTGATACTGTAAttcttc                                                                                                                                 >  2:362667/1‑250 (MQ=255)
                                                                                                                    ccTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCaac                                                                                                                              >  2:189349/1‑250 (MQ=255)
                                                                                                                        aaaT‑CATTTTCTTCTTATTTTTTCTTTCGCATTTTT‑GGAGTTCTTCTTCCGGATT‑GTGTGCCGTTTGTAATATCACTTCTAAACCTT‑CTTTTTCAGGTTTGGCGG‑TGGGTGCTTGAAAACCTT‑AAGAAAAAAACGGATTG‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGc                                                                                                                           <  1:510854/249‑1 (MQ=255)
                                                                                                                         taatcttttttgtgttattttttttttCGCCTTTTT‑GTAGATCTTCTGC‑GGATTTGTGCGACGTTTGTAATATCACAGCGAAACCTT‑CTCCTCCGGTTTTGGCGG‑TGGTGGCAGGAAAAACTG‑CAAAAAAAAACG‑TTTA‑CCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGc                                                                                                                           <  2:127984/226‑1 (MQ=255)
                                                                                                                              aTTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTa                                                                                                                      >  1:85034/1‑250 (MQ=255)
                                                                                                                                tttCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGGCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGGAac                                                                                                                    >  2:610022/1‑250 (MQ=255)
                                                                                                                                  tcgtcTTAATTTTTCTTTC‑CATTTTTTGCAGTTCTTCGGCCGTGTG‑GGGCGATGTTGGTAATATCTCGTAAAAACCTT‑TTTCTTCGGCTTTGGCGG‑TGAGGGAATGAAAAACTT‑CAGAAAAAAACG‑ATTAACCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAacac                                                                                                                  <  2:130871/246‑1 (MQ=255)
                                                                                                                                          aTTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcg                                                                                                         >  1:20902/1‑251 (MQ=255)
                                                                                                                                                    tcgcttttttgaaagtcttctGCCGATTT‑GTGCTTCTTTTGTAATACCGCGTC‑AAACCTTTCTCCTCCACTTTTGGCGG‑TGTTGGTATGCAAATCTT‑CAGAGAACAACGAATTC‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAATCCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCTCGCTTCCg                                                                                                 <  1:155802/235‑1 (MQ=255)
                                                                                                                                                                        tGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGACATGCAAAACTT‑CAGAAAAAAACGGGTTA‑CCC‑TCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTAt                                                                            >  1:158025/1‑251 (MQ=255)
                                                                                                                                                                                       aCGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGAt                                                              >  2:218444/1‑251 (MQ=255)
                                                                                                                                                                                              tAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTgg                                                       >  1:275542/1‑251 (MQ=255)
                                                                                                                                                                                                atatCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCACTGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGCGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAt                                                     >  2:303506/1‑251 (MQ=255)
                                                                                                                                                                                                atatCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAt                                                     >  1:529753/1‑251 (MQ=255)
                                                                                                                                                                                                 tatCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑CGATGGCATGCACAACTT‑CAGAGAAAAGCGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCTCGTAAACCTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAtt                                                    >  1:157193/1‑251 (MQ=255)
                                                                                                                                                                                                  atCACGTCGAAACCTT‑CTTCTTCCGCTTTGGCGG‑TGATGGCCTGCAAAACTT‑CCGAGAAAACCGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCACTGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAttt                                                   >  2:158342/1‑251 (MQ=255)
                                                                                                                                                                                                  atCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAttt                                                   >  1:344500/1‑251 (MQ=255)
                                                                                                                                                                                                   tCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAACTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAtttt                                                  >  1:422916/1‑251 (MQ=255)
                                                                                                                                                                                                          gAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGc                                           >  2:469941/1‑251 (MQ=255)
                                                                                                                                                                                                              ccTT‑CTTCGTCAGTTTGGGCGG‑TGAGGGAATGAAAAACTA‑CAGAAAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGc                                       <  2:529753/251‑1 (MQ=255)
                                                                                                                                                                                                                               tGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGTTACTGTAATTCTTCAATGGCCCGTAAAAATTTTTCGCGCGCTTCCGGGGAGGTATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGGCTGCAACCCTGGCCATATcc                       >  2:382195/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gg‑TGATGGAACGCAAAACAT‑AAAAGAAAAACGGTTTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATaaa                   <  2:584197/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaaaCTT‑CAGAAAAAAACGATTTC‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACAt      <  1:469941/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaCTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGGTTTGGGCGAACGGGTAGAAATAACCATAATGGCATCAAAACCCGTAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGCGATATTAGTTTGCTTCTTCAGCACTCGTGATACGGTGGATTTTGATACGCCTGCAACCTTCGCACTATCATAAATAGTAgggtacct      >  1:121315/1‑244 (MQ=255)
                                                                                                                                                                                                                                                     cTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGt  >  1:550145/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
CAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCCGCTGGCGAACTGAGCATAATAATGCCTTTAAT‑CATTTTCTGCTTAATTTTGCTTTCGCATTTTT‑GCAGATCTTCTGCCGGATT‑GTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGGCGG‑TGATGGCATGCAAAACTT‑CAGAGAAAAACGGATTA‑CCC‑GCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACATAGGT  >  NC_000913/1383475‑1383965

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