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.

CGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACGCGGCG  >  NC_000913/2438608‑2439092
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
cgcgTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCAATACAGATGACGCTAAAAAAACCTTTCGATGCCAGTGAACTGGCGGCGTTCTCGGCATATATCCCGCAACAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTACTTTATCCTTTTTCCCCGattatg                                                                                                                                                                                                                                                >  1:248472/1‑251 (MQ=255)
          ccgaaGCCAGCAAAATACGCCGACGCTGGAATCCCTACAGATGTCGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCCTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGAGCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGc                                                                                                                                                                                                                                      <  2:304901/248‑1 (MQ=255)
                   gCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCGTTTCCACCGATGATGGACACGTTGCAAAATATgg                                                                                                                                                                                                                             <  2:116202/251‑1 (MQ=255)
                                   cTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTCAATGGATTGATCCGCGGGTTTTTGCACCCCTGCCTTATCCTTTTCCACCGCTGATGGACAAGTTGCAAAATAATGCCACGCGGCCCGCGCt                                                                                                                                                                                                            >  1:392711/1‑251 (MQ=255)
                                      gAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCtgtg                                                                                                                                                                                                         <  2:362446/251‑1 (MQ=255)
                                                      cGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTAACGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCCTCGCTTTCTTCTGCCCAGGCCGCGGATATTCTTAATCAGCTTTATCAACGGCGGGTTACAGGGATTTATCCGCGAGTTCTGCAAACGCTGCTTTATACTTTTTCTCCCCTGATTGACACGATTCAAAATATGGCGAACCGTCCCGCGCTGTGGATACCGTTactagcc                                                                                                                                                                                         >  1:220267/1‑245 (MQ=255)
                                                       gCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGa                                                                                                                                                                                        <  2:48766/251‑1 (MQ=255)
                                                           aaaCAACCTTTAGATGCCAGTGAACTGGTGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTAt                                                                                                                                                                                    <  1:65051/251‑1 (MQ=255)
                                                            aacaacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGAATATCTCCAGCAAAAATATGGTCTTTCTGCGAAATCATCGCTTTCTTCCGCCAAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCCACGCGTCCCGCGCTGTGGATACTGTTAGTTGCGATTATc                                                                                                                                                                                   >  1:319592/1‑251 (MQ=255)
                                                               aacCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGATGATCCGCGTGGTATGCAACCGCTTCTTTATTCTTTTTCCCCGATGATGGGAACGGTGCACAATATGGAAACGCGTCCCCCCCCGTGGGGACTGTTAGGCGCGGTTATCCtc                                                                                                                                                                                >  2:468382/1‑250 (MQ=255)
                                                                 ccTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt                                                                                                                                                                              <  1:319760/251‑1 (MQ=255)
                                                                 ccTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGAt                                                                                                                                                                              <  1:468382/251‑1 (MQ=255)
                                                                             aGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGTGTTTTGCAACCGCTGCTTAATCCTTTTCCCCCGATGATGGACCCGTTGCAAAATATTGCAAAGCGTCCCCCGCTGTGGATACTGTTTGTCGCGATTTTCCTGATGCTGGTctggct                                                                                                                                                                  >  1:47678/1‑251 (MQ=255)
                                                                                  aCTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGttt                                                                                                                                                             >  1:176216/1‑250 (MQ=255)
                                                                                      gcggagTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGATAGTCGCGGTTATCCTGATGCTGGTCTCTCTGGTTCGTTa                                                                                                                                                       <  2:46327/251‑1 (MQ=255)
                                                                                             tATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGTGTTATTCAACCTCTGCTTTATCCTTTTTCACCGATGATGGACAAGTTGCAAAATATGGCAACGCGTCCCCCCCCGTGGGTACTGTTAGTCGCGATTTTCCTGATGCTGGTCTGGCGGGGTCGGTAAcccagc                                                                                                                                                >  2:15970/1‑248 (MQ=255)
                                                                                                  gCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGaa                                                                                                                                           >  1:183802/1‑250 (MQ=255)
                                                                                                   cATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAAt                                                                                                                                         <  2:176216/251‑1 (MQ=255)
                                                                                                    atatatCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATg                                                                                                                                        <  2:183802/251‑1 (MQ=255)
                                                                                                     tatatCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTTCAGGGGGTGATCCCCGTGTTTTTCCACCGCCGCTTAATCCTTTTTCACCCATGCTGGGAAAGGTGCAAAAAAAGGCAACGCGTCCCGCGCCGTGGATACTGTTTGTTGCGGTTATTCTGATGCGGGTCTGGCTGGGTCGTTTAAccccccgaaatg                                                                                                                                          >  2:169727/1‑242 (MQ=255)
                                                                                                                aaaaTATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTTATTCTTTTTCACCGATGATGGACACGTTGCCAAATATGGCCACCCCTCCCGCGTTGTGGGTACTGTTAGTCGCGGTTATCCTGATGCTGGGCTGGCTGGTTTGTTAAcccccgcgaaaatacaatattgtgg                                                                                                                             >  2:72727/1‑230 (MQ=255)
                                                                                                                 aaaTATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTTCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCCCGCTGTGGATACTGTTATTCGCGATTATCCTGATGCTGGTCTGGCCGGTTTGTTCACCCCGACGAAATGACAGTATCGTGAc                                                                                                                           >  1:91697/1‑251 (MQ=255)
                                                                                                                     atGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTTAAGGGATTGATCCGCGTGTTTTGCAACCGCTGCTTCATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCCCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTTACCCCGACGAAATGACAGTATCGTGGCaat                                                                                                                        >  1:155628/1‑250 (MQ=255)
                                                                                                                         tCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGAGTGTTATGCAAACGCTCCTTAATCCTTTTTCAACGATAATGCACACGTTGCAAAATATGGCAACGCGCGCCGCGCTGTGGGTACTGTTAGTCGCCATTTTCCTGATGCTGGACTAGCTGGGTCGTTTAACCCCACGAAATGAAAGGATCGTGACAATAATAc                                                                                                                   >  1:214292/1‑251 (MQ=255)
                                                                                                                           tttCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTTATCCCCGTGTTATTTAACACCTGCTTAATTCTTTTTCACCGATGCTGTACCAGTTTCAAAAAATGGTAACCAGGCCCCCCCTGTGGGGAATTTTAGTCGCGATTATTCTCATGCTGGTCTGGGTGGTTCGTTAACCCCCACCGAATGACAGTATcctggtcgaaaaaccc                                                                                                                 >  2:379933/1‑236 (MQ=255)
                                                                                                                             tCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGCTTAAGGGATTGATCCGCGCGTTATGCAACCGCTGCTTAATCCTTTTTCACAGATGATGGACACCTTGCAAAATATGGCAACGCCTCCCGAGCTGTGGAAACCGTTAGACGCGGTTATCCTGCAGGTGGTCTGGCTGGGTAGTTTACCCCGGCGGAATGACAGTATCGTggaaaacatacccac                                                                                                               >  1:377171/1‑237 (MQ=255)
                                                                                                                              cTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGGTTGATCCGCGGGGTTTGCCACCGCTGCTTCATCCTTTTTCACCCCTGATGGACACGTTGCCAAATATTGCAACGCGTCCCCCGCTTGGGAGACTGGTAGTCGCGAGTATCCTGATGCTGGGCTGGCTGGGTCTGTTACCCCGGACAAATGACAGTATCGTGAGAAAAATACCCAGc                                                                                                              >  2:363994/1‑251 (MQ=255)
                                                                                                                                ccgcAATCCTCGCTTTCTTCTGCCCAGGCCGTGGATATTCTTAGTCAGCTTTAGCAACGGCGGGTTAAAGTGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCACAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGACTGTCTGGTTCGTTAACTCCGACGAAATGACCGTATCGTGATAATAATACCCAGCa                                                                                                             <  2:392711/249‑1 (MQ=255)
                                                                                                                                          cGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATc                                                                                                  <  1:72727/251‑1 (MQ=255)
                                                                                                                                              ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCCACGACGGGGTAAAGGGGTTCATCCCACTGTTATGCCAACCCTGCTTAATCCTTTTCCACCGATGATGGACCCGCTGCAAAATATGGCCACGCGTCCCGCCCTGTGGCTACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAAACCCGACGAAATGAACGTATCGTGAAAATAATACCCCGCACCGCACAGATCGcac                                                                                              >  1:40295/1‑249 (MQ=255)
                                                                                                                                              ttcttcTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGc                                                                                                <  2:155628/249‑1 (MQ=255)
                                                                                                                                               tcttctGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCTCGATTATCCTGATGCTGGTCTGGCGGGTTCGGGAACCCCGACGATATGACAGGATCGTGACAATAATACCCATCACCGCAGAGATCGccc                                                                                              <  2:248472/250‑1 (MQ=255)
                                                                                                                                                    tGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATTCAAACGCTGCTTAATCCTTTTCCACCGATGATTGACACGTTGCCAAATATTGCCACCCGTCCCGCGCTGTTGATACTGTTAGTCGCGATTATCCTGCTGCTGGGCTGGCGGGTTCGTTAACCCCGACCGAATTGCAGATTCGTGGCAATAATTCCCAGCCCCCCAGGGATCGCCCCCGCa                                                                                         >  1:99415/1‑250 (MQ=255)
                                                                                                                                                     gCCCAGGCCGAGGATATTCTTAATCAGCTTGATCAACGGCGGGTTAAAGGGATTGATCCGCTTGTTATGCAACCGCGGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTGCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCGTGTCTGGCTGGGTCGTTAACCCCGACGAAATGGCAGTATCGTGGCAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAg                                                                                       <  1:503694/251‑1 (MQ=255)
                                                                                                                                                                ggATAATCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGAGTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAAGATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCGGGTCTGGCTGGTTCGTTAACCCCTACGAAATGACAGTAGCGTGAGAAGAATTCCCAGCGCCGCAGAGATCGCCCCGGCAAGAAATACCGAag                                                                            <  2:530218/251‑1 (MQ=255)
                                                                                                                                                                             cAGCTTTATCAACGGCGGGTGAAAGGGATTGATCCGCGTGTTATGCAACCTCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGTTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGATATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACg                                                                <  1:510237/250‑1 (MQ=255)
                                                                                                                                                                                                                  gtgtTATGCAACCGCTGCTTAATCCTTTTCCACCGAAGATGGACAAGTAACACAATATGGCCACCCGTACCGCGCCATGGATACTGTCAGTAGCGATTATCCTAATGCTGGCCTGGCTGGTTCGATAACCCCCACGAAATGACAGTATCGTGACAATAATACCCAGCCCCGCAGAGATCGCCCAGGCAAGAAATACCGACGAGTAACCAAACGTGCCCGCCAGCATTCCCGCCAGCGGCCCGGAGACGcct                          >  1:31550/1‑250 (MQ=255)
                                                                                                                                                                                                                  gtgtTAGGCAACCGGTGCTTAACCCGTGTCCACCGATGATGGACACGTTGCAAAATATGGGAACGCGGCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCGGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAAGACCGAAGAGTAACCAAACGGGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCg                          <  2:282101/251‑1 (MQ=255)
                                                                                                                                                                                                                  gggtTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGccc                          >  2:416489/3‑250 (MQ=255)
                                                                                                                                                                                                                       aTGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGCGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGc                     >  2:5135/1‑251 (MQ=255)
                                                                                                                                                                                                                       aTGCAACCGCTGCTTAATCCTTTTCAACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGCAATGACAGAATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCACGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGCgggg                     >  2:535601/1‑250 (MQ=255)
                                                                                                                                                                                                                            cacccGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGGTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAAGACCGAAGAGTAACCAAACGGGTTCGCGAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGat                  <  1:363994/248‑1 (MQ=255)
                                                                                                                                                                                                                                cGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCCATATcc              >  2:397809/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gctTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGa           >  1:60496/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aaTCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTACCCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAAcgc      <  1:431901/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCCTTTTCCCCCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGCCTGGTTCGTTACCCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGACCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCACGCATTCCCTCCAGCGGCCCGGACACGCCGAGGGCGATATCCTGAAAcgcg     >  2:203333/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccTTTGCCCCCGAGGAGGGACACGTGGCAAAATATGGCAACGCGTCCCGCGCTTTGGATCCTGTTAGTCGCGATTACCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATCCCGAAGAGTAACCAAACGTGGTCGCCAGCATCCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACgcggc   <  2:21272/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cgtttCCACCGATGATGGACTCGTTGCAAAATAGGTCACCGCGTCCCGCGCTGTGGATACTGTTATGCGCGATTATCCTGATGCTGGTCTGGCTGTTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATCCCCAGCACCGCAGAGATCTCCCCGGCAAGAAATACCGAAGATTACCCAAACGTTGTCGCCATCATTCCCGCCAGCGTCCCGGAGACGCCGAGGGCGATATCCTGAAACgcggcg  <  1:52168/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
CGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACAACCTTTAGATGCCAGTGAACTGGCGGCGTTATCGGCATATATCCAGCAAAAATATGGTCTTTCTGCGCAATCATCGCTTTCTTCTGCCCAGGCCGAGGATATTCTTAATCAGCTTTATCAACGGCGGGTTAAAGGGATTGATCCGCGTGTTATGCAACCGCTGCTTAATCCTTTTCCACCGATGATGGACACGTTGCAAAATATGGCAACGCGTCCCGCGCTGTGGATACTGTTAGTCGCGATTATCCTGATGCTGGTCTGGCTGGTTCGTTAACCCCGACGAAATGACAGTATCGTGACAATAATACCCAGCACCGCAGAGATCGCCCCGGCAAGAAATACCGAAGAGTAACCAAACGTGGTCGCCAGCATTCCCGCCAGCGGCCCGGAGACGCCGAGGGCGATATCCTGAAACGCGGCG  >  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: