Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,383,7220GA30.8% 87.9 / 18.7 52F80F (TTC→TTT) ycjWDNA‑binding transcriptional repressor YcjW
Reads supporting (aligned to +/- strand):  ref base G (27/9);  new base A (0/16);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.97e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACAT  >  NC_000913/1383472‑1383961
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
cGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CGTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCagag                                                                                                                                                                                                                                                 <  1:521/251‑1 (MQ=255)
  aCAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCGTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACTTCGAAACCGT‑CTTCGTCATCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGaa                                                                                                                                                                                                                                               <  2:894716/251‑1 (MQ=255)
             aTGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTAc                                                                                                                                                                                                                                    <  2:130461/251‑1 (MQ=255)
                  cATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTCCccgcc                                                                                                                                                                                                                               <  2:448040/251‑1 (MQ=255)
                                   ttGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGGAAAAAGGATTACCCGCCGTAGTTTTGGTCGAACg                                                                                                                                                                                                              >  1:352140/1‑251 (MQ=255)
                                      ccAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACgggg                                                                                                                                                                                                           >  2:306898/1‑250 (MQ=255)
                                       caatccccccccccGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCACCC‑GCTGGCGAACTGAGCATAATAATGCTTTTAATCTTTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATACCCCGTCAAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAAAAAAACGTATTACCCGCCGTAGTTTTGGTCGAACGGGTa                                                                                                                                                                                                          <  1:721163/239‑1 (MQ=255)
                                             ccccccccgaaaTATCATATTTATCGAGTTGGGAAAAAAATGA‑TCCTCCCGCTGGCGAACTGACCATAAAAATCCTTTTAATTTTTTTTTGCTTAATTTTGTTTTCCCTTTTTTGCAGATCTTCCGCCGTA‑TTGTGCGCCTTTTTTAATATCACGTCGAACCCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCAGGAAAAACTTCAGAAAAAACCGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATc                                                                                                                                                                                                    <  2:617442/241‑1 (MQ=255)
                                                  cccGGAATACTATATTTATCGGGTTGAGAAAAAAATGACTCCTCC‑GCGGGCAATCTGAACATAATAATCCTTTTAATCATTTTCTGTTTAATTTTGCTTTCGATTTTTTGCAGATCTTCTGCCGGT‑TTGTGCGCCTTTTGTAATATCACGTCAAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGT‑TGGAATGCAAAACTTAAGAAAAAAACGTTTTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCAt                                                                                                                                                                                               <  1:146316/250‑1 (MQ=255)
                                                            cATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATc                                                                                                                                                                                      >  1:539012/1‑250 (MQ=255)
                                                             atatTTTTCAAGTTGGGCAAAAACTGCCTCTACC‑GCTGGCGAACTGAGCATAATAATCCCTTTAATCATTTTCTGTTAAATTTTGCTTTCGAATTTTTGCAGATCTGCTGCCGGT‑TTGTGCGGCGTTGTTAATATCACGACGAACCCTT‑CTCCTACAGCTTTGG‑GGGTGA‑TGGCGAAAAAAACTACAGAAAAAAACGTATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCa                                                                                                                                                                                     <  1:1037568/250‑1 (MQ=255)
                                                                   aTCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACggg                                                                                                                                                                                                            >  1:606465/1‑221 (MQ=255)
                                                                           gggCAAAAAATGACTCACCC‑GCTGGCGAACTGAGCATAATAATGCTTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCGGGA‑TTGTGCGACGTTTGAAATATCACGTCGAAACTTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGAAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTc                                                                                                                                                                      <  1:1098033/251‑1 (MQ=255)
                                                                            gcaaaaaaaTGCCTCATCC‑GCGGGCGAACTGAGCATAATAATCCTTTTAATCTTTTTCTCCTTAATTTTTTTTGCGCTTTTTTGCATATCTTCTGCCGGG‑TTGTGCGGCGTTTGAAATATCACGTCGAAACCTT‑CTTCTTCATCTTTGG‑GTGTGA‑TGGCATGCAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCa                                                                                                                                                                     <  2:944237/248‑1 (MQ=255)
                                                                                 aaattgccTCATCC‑GCTGGGGAACGGAGCAAAATAATGCTTTTAATATTTTTCGTCTTAATTTGTCTTTATCTTTTTTGCAGTTCTTGTGCTGGT‑TGGGGCGACTTTTGTAAAATCACGACAAAACCTT‑CTTCTTCAGTTTGGG‑CGGTGA‑TGGATTGAAAATCTTAAAAAAAAAGCGGAATCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAgcgcg                                                                                                                                                                <  2:979739/244‑1 (MQ=255)
                                                                                   aaTGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGGCAGcgcgcg                                                                                                                                                              >  2:594015/1‑251 (MQ=255)
                                                                                             cccgctggcaaacgggaaaaaaaattGCCTTTAATCTTTTTTTGTTTATTTTTGTTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGGGCGACTTTTGAAATACCACGTCGAACCCTT‑CTTCTTCAGTTTGGG‑GGGTGA‑GGGAATGAAAAACTACAAAAAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCtt                                                                                                                                                     <  2:539012/226‑1 (MQ=255)
                                                                                                        aCTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                                          >  2:1085499/1‑251 (MQ=255)
                                                                                                         cTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGAt                                                                                                                                          >  1:669184/1‑250 (MQ=255)
                                                                                                              cATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCTCCGTAGTTTTGGTCGAACGGGTAGaaa                                                                                                                                                                                                      <  2:925266/185‑1 (MQ=255)
                                                                                                              cATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCTCCGTAGTTTTGGTCGAACGGGTAGaaa                                                                                                                                                                                                      >  1:925266/1‑185 (MQ=255)
                                                                                                                       gCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCa                                                                                                                            >  1:1145235/1‑250 (MQ=255)
                                                                                                                          ccgtttttttttCTTCTTATTTTTGCTTCCTCTTTTTTGAAAATCTTCTGCGGTT‑TTGTAGGAGGTTTATAAAACCACGACAAACCCTT‑CTTCTTCAGTTTTGG‑GGGTGA‑TGACATACAAAACTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATgg                                                                                                                        <  2:617029/243‑1 (MQ=255)
                                                                                                                            aaatatttttCTGCTTATTTTTGTTTCCGAATTTTTGAAGATCTTCTGCGGGA‑TTGTGCGACGTTTGTAATACCCCGTCGAACCTTT‑CTTTTTCAGCTTTGG‑CGGTGA‑GGGCATGAAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGc                                                                                                                       <  1:547614/244‑1 (MQ=255)
                                                                                                                                 aTTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTACGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATgg                                                                                                                        >  1:617029/1‑244 (MQ=255)
                                                                                                                                     tttgttaatttttTCTTTCGTATTTTTGCAGTTCTTCTGCGGGT‑TTTGGCGTCGTTTGAAATACCAC‑TTAAAACCTTTCTTCTTCA‑TTTTGG‑GGGGGAAGGGAATGCAAACCTACAGAGAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAacac                                                                                                              <  2:787501/241‑1 (MQ=255)
                                                                                                                                          ttatttttGCTTTCGCTTTTTTGCAGATCTTCTGCCGGT‑TTGTGGTACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGAAAAACTTAAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACttttt                                                                                                         <  1:534904/246‑1 (MQ=255)
                                                                                                                                            aaTTTTTCTTCCGCTTTTTTGCGGATCTTCTGC‑GGAGTTGTGCGTCGTTTGTAATATCACGTCAAAACCTC‑CTTCTCCACTTTTGG‑CGGTGA‑TGCCATAAAAAACTCCAAAGAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGTTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcg                                                                                                       <  1:730516/250‑1 (MQ=255)
                                                                                                                                            aaTTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgc                                                                                                      >  2:1087191/1‑251 (MQ=255)
                                                                                                                                              ttttGCTTTCGCATTTTTGCAGATTTTCTGCCGGT‑TTGTGCGCCGTTTGAAATATCACGTCAAACCCTT‑CTTCTCCGGTTTTGG‑GGGTGA‑TGGCATACAAAATTAAAGAAAAAAACGTATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTcgcg                                                                                                     <  2:1145235/250‑1 (MQ=255)
                                                                                                                                                  tctttcctttttttgacgATCTTTTGCCGGT‑TTGGGCGGCGTTTAAAATACCCCGTCAAACCCTT‑CTTCTTCAGCTTGGG‑GGTTGA‑GGGCATGAAAAACTAAAAAAAAAAACGTTTTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCt                                                                                                 <  2:1075354/233‑1 (MQ=255)
                                                                                                                                                   cTTTCGCCTTATTGCAGATCTTCTGCCGGA‑TCGTGCGACGTTTGTACTATCACGTGGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTGGATTTGGTCGAAAGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGAGGGCCAGCCTGTTTGGCTGATCCTGTAATTCTTCCAAGGCCCTTCACACCTTTTCGCGCGCttt                                                                                               >  2:425184/1‑250 (MQ=255)
                                                                                                                                                    tttCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTcc                                                                                              >  2:71301/1‑251 (MQ=255)
                                                                                                                                                     ttcttttttttGCAGATCTGCTCCCGGT‑TTGTGCGCCTTTTGAAATACCACGCCGAACCCTT‑CTTCTCCA‑CTTTGGGCGGTGA‑GGACATACAAACCCTCAAAGAAAACCGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTcc                                                                                              <  1:411667/244‑1 (MQ=255)
                                                                                                                                                     ttCGCATTTTTGAAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTAAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTTGCTGCTACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCg                                                                                             >  2:1069143/1‑251 (MQ=255)
                                                                                                                                                      ggtgttctcttccgATCTTCTGCGGGA‑TTGTGCGACGTTTGTAATACCACGTCGAACCCTT‑CTTCTTCAGCTTTGG‑GGGTGA‑GGGCATGAAAACCTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCgg                                                                                            <  1:365824/238‑1 (MQ=255)
                                                                                                                                                             tttGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGat                                                                                      >  1:573224/1‑250 (MQ=255)
                                                                                                                                                                  aGATCTTCTGACGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGt                                                                                >  1:509060/1‑251 (MQ=255)
                                                                                                                                                                        tctGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCgg                                                                                            >  2:365824/1‑233 (MQ=255)
                                                                                                                                                                             cGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATAACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGATTGCTTATTc                                                                      >  2:1116994/1‑250 (MQ=255)
                                                                                                                                                                              ggA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAAAACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTc                                                                      >  2:1076465/1‑249 (MQ=255)
                                                                                                                                                                                     ttcgacgaTTGTAATATCCCGTCGAAACCTT‑CTTCTGCAGCTTTGG‑CGGAAC‑TGGCATGCAAAAATTCAAAGAAAAACGGATTACCCGCCGTAGGTTTGGTCGAACGGGTAGAAATCACCATAATGGCAACAAAACCCGAAGAGGTCCGCGCGCGGGCCAGCCTGTTTGGCTGATACTATAATTATTCAATGGCCCGTCACCCTTTTTCGCGGGCTTCCGGGGAGATATTCGTTTGCTTATTCATCAcgcg                                                              >  2:975829/3‑251 (MQ=255)
                                                                                                                                                                                         aCGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGAt                                                          >  2:639095/1‑251 (MQ=255)
                                                                                                                                                                                             ttGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGCGC‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCCTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACgg                                                      >  1:425243/1‑251 (MQ=255)
                                                                                                                                                                                               gTAAGATCACGACTAACCCTT‑CTCCTTCATCTTTGG‑CGGTGT‑GGGCAGGCAACCTTTCAAAAAAAAACGAATCCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTAGTCAGCACGCGTGATACGGt                                                     <  2:422546/250‑1 (MQ=255)
                                                                                                                                                                                                    atCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGAttt                                               <  1:594015/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCt                                     >  1:1062080/1‑251 (MQ=255)
                                                                                                                                                                                                                      ttcttcAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAcccc                              >  1:1125143/1‑251 (MQ=255)
                                                                                                                                                                                                                            cagtttggggggggA‑TGGCATGAAAAACTTCAAAAAAAAACGTATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCa                         <  1:290915/238‑1 (MQ=255)
                                                                                                                                                                                                                               tttGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAatat                     >  1:730784/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggA‑TGGCATGAAAAATTTCAGAAAAAAACGGATTCCCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAAta             <  1:1076465/250‑1 (MQ=255)
                                                                                                                                                                                                                                               cATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAggg       >  1:122804/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     aaaCTTCAGAGAAAAACGGATTACCCGCAGTAGTTTTGGTCGAACGGGTAGAAATAAACATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATAAGCCTGCAACCCTGGCCATATAATAAATAGCAGGGGAcgt  >  2:234959/1‑248 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
CGACAGAATAAACATGGGCATATTGACCTTCAACTTTGCCAATCACCACCACCGGAATATCATATTTATCGAGTTGGGCAAAAAATGACTCATCC‑GCTGGCGAACTGAGCATAATAATGCCTTTAATCATTTTCTGCTTAATTTTGCTTTCGCATTTTTGCAGATCTTCTGCCGGA‑TTGTGCGACGTTTGTAATATCACGTCGAAACCTT‑CTTCTTCAGCTTTGG‑CGGTGA‑TGGCATGCAAAACTTCAGAGAAAAACGGATTACCCGCCGTAGTTTTGGTCGAACGGGTAGAAATCACCATAATGGCATCAAAACCCGAAGAGGTCAGCGCGCGGGCCAGCTTGTTTGGCTGATACTGTAATTCTTCAATGGCCCGTAACACTTTTTCGCGCGCTTCCGGGGAGATATTAGTTTGCTTATTCAGCACGCGTGATACGGTGGATTTTGATACGCCTGCAACCCTGGCAATATCATAAATAGTAGGGGACAT  >  NC_000913/1383472‑1383961

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