Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,051,826	G→T	G2297W (GGG→TGG)	yeeJ →	inverse autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,051,826	0	G	T	97.8%	168.9 / NA	45	G2297W (GGG→TGG)	yeeJ	inverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (21/23);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGTTGATGCTACCGGTAAAGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCG‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTT  >  NC_000913/2051590‑2052064                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                aTGTTGATGCTACCGGTAAAGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGt                                                                                                                                                                                                                                  >  1:362980/1‑251 (MQ=255)   gTTGATGCTACCGGTAAAGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGa                                                                                                                                                                                                                                <  2:18759/251‑1 (MQ=255)                gTAAAGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATg                                                                                                                                                                                                                   <  2:416019/251‑1 (MQ=255)                 aaaaGTGACATTTAAAAATGGCGGCAGCAATTCGGAAAGGATTGCGGCGACGCCAAAATCAGGAGACGCGAGCGATGGATACGAAATCCGTGTGAATAGTTGGTGGGTGAACCCCGGCGAGGCTTTCATGAGATACAGCCTTGCAGAAAATTTATGCATCATCAATGGCTACACTCTCCCCAGAGCAAACTATTTAAACCACTGTAGGTTCCTAGGCA‑CGTGCTCGCTGTACAGTGAATGGGGTGATATgg                                                                                                                                                                                                                  <  2:380011/250‑1 (MQ=255)                  aaaGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATggg                                                                                                                                                                                                                 >  1:101976/1‑251 (MQ=255)                                      cGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCt                                                                                                                                                                                              <  1:36621/250‑1 (MQ=255)                                         cagcaATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGtt                                                                                                                                                                                          <  1:288130/251‑1 (MQ=255)                                              aTTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCGACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTtca                                                                                                                                                                                      >  1:345893/1‑250 (MQ=255)                                                  ggAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTtcaatca                                                                                                                                                                                  >  2:394343/1‑250 (MQ=255)                                                  ggAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTtcaatca                                                                                                                                                                                  >  2:319376/1‑250 (MQ=255)                                                   gAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTtcaatcaa                                                                                                                                                                                 >  2:57684/1‑250 (MQ=255)                                                   gAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTtcaatcaa                                                                                                                                                                                 >  1:41703/1‑250 (MQ=255)                                                    aaaGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATCCGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAat                                                                                                                                                                               <  2:344953/251‑1 (MQ=255)                                                              ggCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCTCGAGCCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGt                                                                                                                                                                     <  1:148955/251‑1 (MQ=255)                                                              ggCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGt                                                                                                                                                                     >  1:173561/1‑251 (MQ=255)                                                                   cGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCt                                                                                                                                                                >  2:144565/1‑251 (MQ=255)                                                                    gCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTa                                                                                                                                                               <  2:32173/251‑1 (MQ=255)                                                                         aaTCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTccc                                                                                                                                                          >  1:407582/1‑251 (MQ=255)                                                                              ggaggCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCaaa                                                                                                                                                     <  1:394343/251‑1 (MQ=255)                                                                                   cccTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCa                                                                                                                                                 >  2:52985/1‑250 (MQ=255)                                                                                    ccTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAg                                                                                                                                               <  1:140280/251‑1 (MQ=255)                                                                                           aTGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAAc                                                                                                                                          <  2:271357/249‑1 (MQ=255)                                                                                                   aaaaTCCGGGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTGGCTGAAAATATGTGCAGCAGCCATGGCTACACGCTCCCCAGAGCAAGGTATTTAAACCACTGTAGTTCGCGAGGCACCT‑GGTCAGTGTTGAGTGAATGGGGAGATATGGGGCATTACTCGATTGACGCTGGTGTTCAGTCGAAGATGTATTGTTCCTCTAGTCCCGCATACTGAAGCCAAGAATAGGTAGt                                                                                                                                <  2:329022/250‑1 (MQ=255)                                                                                                      aTCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTc                                                                                                                             >  2:81667/1‑251 (MQ=255)                                                                                                              gtagagTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCa                                                                                                                      >  2:416048/3‑250 (MQ=255)                                                                                                                     tggtggGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTg                                                                                                               <  2:359927/250‑1 (MQ=255)                                                                                                                     tggtggGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTg                                                                                                               <  2:362980/250‑1 (MQ=255)                                                                                                                       gtgggtgAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATc                                                                                                            >  2:272714/1‑251 (MQ=255)                                                                                                                       gtgggtgAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATc                                                                                                            <  2:41703/251‑1 (MQ=255)                                                                                                                                cGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTAt                                                                                                   <  2:352724/251‑1 (MQ=255)                                                                                                                                  ccGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTAttt                                                                                                 <  1:212515/251‑1 (MQ=255)                                                                                                                                                  gATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTGCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGttt                                                                                  <  1:166689/250‑1 (MQ=255)                                                                                                                                                     ataCAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCtt                                                                              >  2:366392/1‑251 (MQ=255)                                                                                                                                                                  gAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTta                                                                 <  1:50166/251‑1 (MQ=255)                                                                                                                                                                       tttttGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATaaa                                                             >  2:231878/1‑250 (MQ=255)                                                                                                                                                                         tttGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATaaaaa                                                           <  1:231878/250‑1 (MQ=255)                                                                                                                                                                          ttGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAAc                                                          <  1:307955/250‑1 (MQ=255)                                                                                                                                                                            gcagcagcaATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCt                                                        <  1:4564/250‑1 (MQ=255)                                                                                                                                                                                 agcaATGCATACACGCCCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGCGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGCACAATACGTAGTTTCCCTGGCAACAGGTGATCAACGCGTATTTGAAAATCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGAtt                                                  >  2:155581/1‑251 (MQ=255)                                                                                                                                                                                          tACACGCTCCCGAGAGCAAACTATTTAAACCACTGTAGTTCCCGCGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTCTTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTtat                                         >  2:362462/1‑251 (MQ=255)                                                                                                                                                                                           acacGCGCCCCATAGCAAACTATTTAAACCACGGTAGTGCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAGCAATACGTAGTTTCCCGGGCATCAGGTGATCAAAGCTTATGTGAAAAGCTTGAGTTGGCTTATGCGACATGGTATAAAAACCGGGGATTTTCCGTta                                          <  2:340869/249‑1 (MQ=255)                                                                                                                                                                                                           aaaCTATTTAAACCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAgtgt                        >  2:437723/1‑251 (MQ=255)                                                                                                                                                                                                                    aaaCCACTGTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCgg               >  1:379192/1‑251 (MQ=255)                                                                                                                                                                                                                            gTAGTTCCCGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGTGGCActct       <  2:131399/251‑1 (MQ=255)                                                                                                                                                                                                                                  cccGAGGCATCT‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCtttttt  >  1:265339/1‑250 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ATGTTGATGCTACCGGTAAAGTGACATTTAAAAATGTCGGCAGCAATTCGGAAAGGATTACGGCGACGCCAAAATCAGGAGGCCCTAGCTATGTATACGAAATCCGTGTGAAGAGTTGGTGGGTGAACGCCGGCGAGGCTTTCATGATATACAGCCTTGCTGAAAATTTTTGCAGCAGCAATGGCTACACGCTCCCCAGAGCAAACTATTTAAACCACTGTAGTTCCCGAGGCATCG‑GGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTT  >  NC_000913/2051590‑2052064

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