Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,8470CT22.9% 78.8 / 11.6 48P300S (CCA→TCA) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (15/22);  new base T (11/0);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 3.73e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACGCGGCG  >  NC_000913/2438608‑2439092
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
cgcgTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCAATACAGATGACGCTAAAAAAACCTTTCGATGCCAGTGAACTGGCGGCGTTCTCGGCATATATCCCGCAACAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTACTTTATCCTTTTTCCCCGattatg                                                                                                                                                                                                                                                >  1:248472/1‑251 (MQ=255)
          ccgaaGCCAGCAAAATACGCCGACGCTGGAATCCCTACAGATGTCGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCCTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGAGCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGc                                                                                                                                                                                                                                      <  2:304901/248‑1 (MQ=255)
                   gCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCGTTTCCACCGATGATGGACACGTTGCAAAATATgg                                                                                                                                                                                                                             <  2:116202/251‑1 (MQ=255)
                                   cTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTCAATGGATTGATCCGCGGGTTTTTGC‑A‑CCCCTGCCTTATCCTTTTCCACCGCTGATGGACAAGTTGCAAAATAATGC‑CACGCGGCCCGCGCt                                                                                                                                                                                                            >  1:392711/1‑251 (MQ=255)
                                      gAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCtgtg                                                                                                                                                                                                         <  2:362446/251‑1 (MQ=255)
                                                      cGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTAACGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCCTCGCTTTCTTCTGCCCAGGCCGCGGATATTCTTAATCAGCTTTATCAACGGCGGGTTACAGGGATTTATCCGCGAGTTCT‑GC‑AAACGCTGCTTTATACTTTTTCTCCCCTGATTGACACGATTCAAAATATGGCGAAC‑CGTCCCGCGCTGTGGATACCGTTactagcc                                                                                                                                                                                         >  1:220267/1‑245 (MQ=255)
                                                       gCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGa                                                                                                                                                                                        <  2:48766/251‑1 (MQ=255)
                                                           aaaCAACCTTTAGATGCCAGTGAACTGGTGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTAt                                                                                                                                                                                    <  1:65051/251‑1 (MQ=255)
                                                            aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGAATATCTCCAGCAAAAATATGGTCTTTCTGCGAAATCATCGCTTTCTTCCGCCAAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑CACGCGTCCCGCGCTGTGGATACTGTTAGTTGCGATTATc                                                                                                                                                                                   >  1:319592/1‑251 (MQ=255)
                                                               aacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGATGATCCGCGTGGTAT‑GC‑AACCGCTTCTTTATTCTTTTTCCCCGATGATGGGAACGGTGCACAATATGGA‑AACGCGTCCCCCCCCGTGGGGACTGTTAGGCGCGGTTATCCtc                                                                                                                                                                                >  2:468382/1‑250 (MQ=255)
                                                                 ccTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt                                                                                                                                                                              <  1:319760/251‑1 (MQ=255)
                                                                 ccTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt                                                                                                                                                                              <  1:468382/251‑1 (MQ=255)
                                                                             aGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGTGTTTT‑GC‑AACCGCTGCTTAATCCTTTTCCCCCGATGATGGACCCGTTGCAAAATATTGC‑AAAGCGTCCCCCGCTGTGGATACTGTTTGTCGCGATTTTCCTGATGCTGGTctggct                                                                                                                                                                  >  1:47678/1‑251 (MQ=255)
                                                                                  aCTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGttt                                                                                                                                                             >  1:176216/1‑250 (MQ=255)
                                                                                      gcggagTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGATAGTCGCGGTTATCCTGATGCTGGTCTCTCTGGTTC‑GTT‑a                                                                                                                                                       <  2:46327/251‑1 (MQ=255)
                                                                                             tATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGTGTTAT‑TC‑AACCTCTGCTTTATCCTTTTTCACCGATGATGGACAAGTTGCAAAATATGGC‑AACGCGTCCCCCCCCGTGGGTACTGTTAGTCGCGATTTTCCTGATGCTGGTCTGGCGGGGTC‑GGT‑AAcccagc                                                                                                                                                >  2:15970/1‑248 (MQ=255)
                                                                                                  gCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑Gaa                                                                                                                                           >  1:183802/1‑250 (MQ=255)
                                                                                                   cATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAAt                                                                                                                                         <  2:176216/251‑1 (MQ=255)
                                                                                                    atatatCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATg                                                                                                                                        <  2:183802/251‑1 (MQ=255)
                                                                                                     tatatCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTTCAGGGGGTGATCCCCGTGTTTT‑TC‑CACCGCCGCTTAATCCTTTTTCACCCATGCTGGGAAAGGTGCAAAAAAAGGC‑AACGCGTCCCGCGCCGTGGATACTGTTTGTTGCGGTTATTCTGATGCGGGTCTGGCTGGGTC‑GTTTAAccccccgaaatg                                                                                                                                          >  2:169727/1‑242 (MQ=255)
                                                                                                                aaaaTATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTTATTCTTTTTCACCGATGATGGACACGTTGCCAAATATGGC‑CACCCCTCCCGCGTTGTGGGTACTGTTAGTCGCGGTTATCCTGATGCTGGGCTGGCTGGTTT‑GTT‑AAcccccgcgaaaatacaatattgtgg                                                                                                                             >  2:72727/1‑230 (MQ=255)
                                                                                                                 aaaTATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTTCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCCCGCTGTGGATACTGTTATTCGCGATTATCCTGATGCTGGTCTGGCCGGTTT‑GTT‑CACCCCGAC‑GAAATGACAGTATCGTGAc                                                                                                                           >  1:91697/1‑251 (MQ=255)
                                                                                                                     atGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTTAAGGGATTGATCCGCGTGTTTT‑GC‑AACCGCTGCTTCATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCCCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑TACCCCGAC‑GAAATGACAGTATCGTGGCaat                                                                                                                        >  1:155628/1‑250 (MQ=255)
                                                                                                                         tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGAGTGTTAT‑GC‑AAACGCTCCTTAATCCTTTTTCAACGATAATGCACACGTTGCAAAATATGGC‑AACGCGCGCCGCGCTGTGGGTACTGTTAGTCGCCATTTTCCTGATGCTGGACTAGCTGGGTC‑GTTTAACCCC‑AC‑GAAATGAAAGGATCGTGACAATAATAc                                                                                                                   >  1:214292/1‑251 (MQ=255)
                                                                                                                           tttCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTTATCCCCGTGTTAT‑TT‑AACACCTGCTTAATTCTTTTTCACCGATGCTGTACCAGTTTCAAAAAATGGT‑AACCAGGCCCCCCCTGTGGGGAATTTTAGTCGCGATTATTCTCATGCTGGTCTGGGTGGTTC‑GTT‑AACCCCCAC‑CGAATGACAGTATcctggtcgaaaaaccc                                                                                                                 >  2:379933/1‑236 (MQ=255)
                                                                                                                             tCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGCTTAAGGGATTGATCCGCGCGTTAT‑GC‑AACCGCTGCTTAATCCTTTTTCACAGATGATGGACACCTTGCAAAATATGGC‑AACGCCTCCCGAGCTGTGGAAACCGTTAGACGCGGTTATCCTGCAGGTGGTCTGGCTGGGTA‑GTT‑TACCCCGGC‑GGAATGACAGTATCGTggaaaacatacccac                                                                                                               >  1:377171/1‑237 (MQ=255)
                                                                                                                              cTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGGGGTTT‑GC‑CACCGCTGCTTCATCCTTTTTCACCCCTGATGGACACGTTGCCAAATATTGC‑AACGCGTCCCCCGCTTGGGAGACTGGTAGTCGCGAGTATCCTGATGCTGGGCTGGCTGGGTCTGTT‑ACCCCGGAC‑‑AAATGACAGTATCGTGAGAAAAATACCCAGc                                                                                                              >  2:363994/1‑251 (MQ=255)
                                                                                                                                ccgcAATCCTCGCTTTCTTCTGCCCAGGCCGTGGATATTCTTAGTCAGCTTTAGCAACGGCGGGTTAAAGTGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCACAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGACTGTCTGGTTC‑GTT‑AACTCCGAC‑GAAATGACCGTATCGTGATAATAATACCCAGCa                                                                                                             <  2:392711/249‑1 (MQ=255)
                                                                                                                                          cGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATc                                                                                                  <  1:72727/251‑1 (MQ=255)
                                                                                                                                              ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCCACGACGGGGTAAAGGGGTTCATCCCACTGTTAT‑GCCAACC‑CTGCTTAATCCTTTTCCACCGATGATGGACCCGCTGCAAAATATGGC‑CACGCGTCCCGCCCTGTGGCTACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AAACCCGAC‑GAAATGAACGTATCGTGAAAATAATACCCCGCACCGCACAGATCGcac                                                                                              >  1:40295/1‑249 (MQ=255)
                                                                                                                                              ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGc                                                                                                <  2:155628/249‑1 (MQ=255)
                                                                                                                                               tcttctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCTCGATTATCCTGATGCTGGTCTGGCGGGTTC‑GGG‑AACCCCGAC‑GATATGACAGGATCGTGACAATAATACCCATCACCGCAGAGATCGccc                                                                                              <  2:248472/250‑1 (MQ=255)
                                                                                                                                                    tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑TC‑AAACGCTGCTTAATCCTTTTCCACCGATGATTGACACGTTGCCAAATATTGC‑CACCCGTCCCGCGCTGTTGATACTGTTAGTCGCGATTATCCTGCTGCTGGGCTGGCGGGTTC‑GTT‑AACCCCGACCGAATTG‑CAGATTCGTGGCAATAATTCCCAGCCCCCCAGGGATCGCCCCCGCa                                                                                         >  1:99415/1‑250 (MQ=255)
                                                                                                                                                     gCCCAGGCCGAGGATATTCTTAATCAGCTTGATCAACGGCGGGTTAAAGGGATTGATCCGCTTGTTAT‑GC‑AACCGCGGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTGCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCGTGTCTGGCTGGGTC‑GTT‑AACCCCGAC‑GAAATGGCAGTATCGTGGCAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAg                                                                                       <  1:503694/251‑1 (MQ=255)
                                                                                                                                                                ggATAATCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGAGTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAAGATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCGGGTCTGGCTGGTTC‑GTT‑AACCCCTAC‑GAAATGACAGTAGCGTGAGAAGAATTCCCAGCGCCGCAGAGATCGCCCCGGCAAGAAATACCGAag                                                                            <  2:530218/251‑1 (MQ=255)
                                                                                                                                                                             cAGCTTTATCAACGGCGGGTGAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCTCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGTTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGATATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACg                                                                <  1:510237/250‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGAAGATGGACAAGTAACACAATATGGC‑CACCCGTACCGCGCCATGGATACTGTCAGTAGCGATTATCCTAATGCTGGCCTGGCTGGTTC‑GAT‑AACCCCCAC‑GAAATGACAGTATCGTGACAATAATACCCAGCCCCGCAGAGATCGCCCAGGCAAGAAATACCGACGAGTAACCAAACGTGCCCGCCAGCATTCCCGCCAGCGGCCCGGAGACGcct                          >  1:31550/1‑250 (MQ=255)
                                                                                                                                                                                                                  gtgtTAG‑GC‑AACCGGTGCTTAACCCGTGTCCACCGATGATGGACACGTTGCAAAATATGGG‑AACGCGGCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCGGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAAGACCGAAGAGTAACCAAACGGGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCg                          <  2:282101/251‑1 (MQ=255)
                                                                                                                                                                                                                  gggtTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGccc                          >  2:416489/3‑250 (MQ=255)
                                                                                                                                                                                                                       aT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGCGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGc                     >  2:5135/1‑251 (MQ=255)
                                                                                                                                                                                                                       aT‑GC‑AACCGCTGCTTAATCCTTTTCAACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GCAATGACAGAATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCACGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGCgggg                     >  2:535601/1‑250 (MQ=255)
                                                                                                                                                                                                                            cacccGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGGTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAAGACCGAAGAGTAACCAAACGGGTTCGCGAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGat                  <  1:363994/248‑1 (MQ=255)
                                                                                                                                                                                                                                cGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCCATATcc              >  2:397809/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGa           >  1:60496/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aaTCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑ACCCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAAcgc      <  1:431901/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCCTTTTCCCCCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGCCTGGTTC‑GTT‑ACCCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGACCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCACGCATTCCCTCCAGCGGCCCGGACACGCCGAGGGCGATATCCTGAAAcgcg     >  2:203333/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccTTTGCCCCCGAGGAGGGACACGTGGCAAAATATGGC‑AACGCGTCCCGCGCTTTGGATCCTGTTAGTCGCGATTACCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATCCCGAAGAGTAACCAAACGTGGTCGCCAGCATCCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACgcggc   <  2:21272/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cgtttCCACCGATGATGGACTCGTTGCAAAATAGGTC‑ACCGCGTCCCGCGCTGTGGATACTGTTATGCGCGATTATCCTGATGCTGGTCTGGCTGTTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATCCCCAGCACCGCAGAGATCTCCCCGGCAAGAAATACCGAAGATTACCCAAACGTTGTCGCCATCATTCCCGCCAGCGTCCCGGAGACGCCGAGGGCGATATCCTGAAACgcggcg  <  1:52168/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
CGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTAT‑GC‑AACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGC‑AACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTC‑GTT‑AACCCCGAC‑GAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACGCGGCG  >  NC_000913/2438608‑2439092

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