Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,129,796	G→A	S12F (TCC→TTC)	flgN ←	flagellar biosynthesis protein FlgN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,129,796	0	G	A	97.7%	131.4 / NA	45	S12F (TCC→TTC)	flgN	flagellar biosynthesis protein FlgN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/19);  minor base C (1/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTT  >  NC_000913/1129552‑1130033                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑Ct                                                                                                                                                                                                                                          >  1:48486/1‑251 (MQ=255)                  ggCGTAGTTGCTTCGTTTTCCCAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAAGCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTTCCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGc                                                                                                                                                                                                                          <  2:397965/250‑1 (MQ=255)                  ggCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGc                                                                                                                                                                                                                          >  1:355744/1‑250 (MQ=255)                    cGTAGTTGCTGCGGTTTCACAGTAATTTCCCGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTGCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGGTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAg                                                                                                                                                                                                                       <  2:202797/251‑1 (MQ=255)                                                   gCCAGCGGTGACTAATATCAACGCTATTTGCTGGATTGGGTTCTTTCCCGCGTAACTGTGCGAGGTAATCCATCGTCGCCAGCATTGAGCTTTTATGTTCGGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTGTGAGATCGTTAAGCACAGCTAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGTTCAGGTACTCTGCAAGt                                                                                                                                                                                        <  2:493964/250‑1 (MQ=255)                                                    ctaGCGTTGACTAATATCTACGCTATTTGCTGTCTTGGGTTCTTTCCTGCGTAACTGATCGAGGTAATACAGCGTAGACAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTCCTCTTGATCCATTACCGTTTTGCGATCGTTAAGCACAGCG‑ATATTCTGGTCGAGGATCTCTGCATGAA‑GTGTCATACGAACAGTTACTCTGCAAGTCt                                                                                                                                                                                      >  2:189669/3‑251 (MQ=255)                                                    ccAGCGTTGACTAATATCAACGCTATTTTCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCACCATAGAGTGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCt                                                                                                                                                                                      <  1:432556/251‑1 (MQ=255)                                                    ccAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCt                                                                                                                                                                                      >  1:169241/1‑251 (MQ=255)                                                        cGTTGACTGATATCAACGCTGTTTGCTGTGTTGGGTTCTTTCATGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTGTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGAACCAGTACCGTTTTACGATAGTACAGCACCGCGAACAT‑CTGGTCGAGGATCTCTGCAA‑ACAGAGTCGTACGCTCAGTTAATCTACAAGTCTtgct                                                                                                                                                                                  >  2:167978/1‑251 (MQ=255)                                                             aCTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCt                                                                                                                                                                             >  1:213840/1‑251 (MQ=255)                                                                     cAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGACCAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGAt                                                                                                                                                                     >  2:191141/1‑251 (MQ=255)                                                                       aCGCTATTTGCTGTATTGGGTTCTTTCCTGCGTCACTGTTCGAGGACATCCCGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGAAATCCATAGCAACTGGCAGCCGTTGATCTGCCGCATAGAGCGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGCACAT‑CTGGTCGAGGGTCTCTGCAAGAC‑GGGTCATCCGATCAGTTACTCTGCAAGTCTTGCTGCGCTGCGTTGATc                                                                                                                                                                    >  2:504422/1‑250 (MQ=255)                                                                           tATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATAAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAgcgc                                                                                                                                                               >  2:81947/1‑251 (MQ=255)                                                                                 cTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGc                                                                                                                                                         >  1:400268/1‑251 (MQ=255)                                                                                  tGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTCTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATACGCGCATCGGc                                                                                                                                                         >  2:291793/1‑250 (MQ=255)                                                                                   gTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGc                                                                                                                                                         <  1:482065/249‑1 (MQ=255)                                                                                    tATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCCTTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTTTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCATCGCATCGGCAAt                                                                                                                                                      >  2:367395/1‑251 (MQ=255)                                                                                       tGGGTTCTTTCCTGAGTAACGTTTCGAGGTAATCCAGCGTCGCCAGCAGTGATCTTTTTTGTTCTGTAATCCATTGCCTCTGGCTGCCGTTGATCTGCCCAATAGAGAGATGTTGCTGCTCTTGATCCATTTCCGTTTTGAGATCGTTAAGCACAGCTAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTTTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTGGACCAGCGCATCGGCAAttt                                                                                                                                                    <  2:340242/250‑1 (MQ=255)                                                                                            tCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCGTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCg                                                                                                                                               <  2:263724/250‑1 (MQ=255)                                                                                                 ccTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTc                                                                                                                                          <  2:94624/250‑1 (MQ=255)                                                                                                    gCGTAACTGTTTGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATATAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCa                                                                                                                                        <  1:251032/249‑1 (MQ=255)                                                                                                               cGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCAcc                                                                                                                           >  1:342768/1‑251 (MQ=255)                                                                                                                    tAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGCGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTAc                                                                                                                      >  1:224422/1‑251 (MQ=255)                                                                                                                                gCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGtt                                                                                                          <  1:332556/251‑1 (MQ=255)                                                                                                                                 ccAGCAGTGAGCTGTTGTGTTCTGGAATCCATTTCAACTGGCTGTCGTTTATCTGCCCCATAGATAGAAGGTGCTGCTCTTGATCCCTTTCCGTTTTGAGATCGTGAAGCACAGCGTACAG‑CTGGTCGAGGATCTCTGCAAGAC‑GGGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGCTGTTCATCTTTAGTCCACCGTGACTACTCGCTAGttt                                                                                                         <  2:113746/251‑1 (MQ=255)                                                                                                                                  cagcagTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATAAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGtttt                                                                                                        >  1:4960/1‑251 (MQ=255)                                                                                                                                     cagTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAAc                                                                                                     <  1:15444/251‑1 (MQ=255)                                                                                                                                     cagTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAAc                                                                                                     >  1:169953/1‑251 (MQ=255)                                                                                                                                             ttttttGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGac                                                                                             >  2:421681/1‑251 (MQ=255)                                                                                                                                                    ttCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACAGTTAAGAATCGCCAGTTTTAACGCTTCTACACGTTc                                                                                       >  1:321571/1‑250 (MQ=255)                                                                                                                                                                   cAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACATCGACCAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCActg                                                                       <  2:1344/251‑1 (MQ=255)                                                                                                                                                                        ggCTGCCGTTGATCTGCACCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACCT‑CTGGTCGAGGATCTCTGCAAGCC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGAGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCCCCGTTCCGAATCGCCCGTTTTAACTCTTCGACACGTTCAAGATTGATATCACTGCTGcc                                                                  >  2:493727/1‑251 (MQ=255)                                                                                                                                                                             ccGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGtt                                                             <  1:499033/251‑1 (MQ=255)                                                                                                                                                                                                  agaTGTTGCTGCTCTTGCTCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGt                                        >  1:435592/1‑251 (MQ=255)                                                                                                                                                                                                   gaTGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTc                                       >  1:250814/1‑251 (MQ=255)                                                                                                                                                                                                         gctgctCTTGATCCATTACCGTTTTTAGATCGTTAAGCACAGCGAACTT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACTATCAGTTACTCTGCAAGTCTTGCTTCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCTACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTa                                  <  1:360157/250‑1 (MQ=255)                                                                                                                                                                                                          ctgctCTTGATCCATTACCGTTTTGAGCTCGTTAAGCACCGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATAAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGTCAATTTTGCCGGTGTCCATTTTTATTTCACCGTTACGAATAGCCAGTTTTAACGCTTAGTCACGTTCAATATTGATAGCACTGCTGCCGGGTTGGCTAAGTTTTTCTTGCGCGTCGCTTaa                                 >  2:460445/1‑250 (MQ=255)                                                                                                                                                                                                          ctgctCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAAc                                >  2:460459/1‑251 (MQ=255)                                                                                                                                                                                                           tgctCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGACCAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACg                               >  2:328229/1‑251 (MQ=255)                                                                                                                                                                                                                   attcATTACCTGTTTGAGCTCGTTAAGTACAGCTAACTT‑CTGTTCGATGATCTCTGAACGGC‑GTGTCATACGAGTAGTTACTCTGCACGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGGCCATTTTTAGTTCACCGTTCCGATTCGCCAGTTTTAACGCTTCGTCACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACAct                        <  2:408940/247‑1 (MQ=255)                                                                                                                                                                                                                         tACCGTTTTGAGATCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACActggtgct                  <  1:203669/250‑1 (MQ=255)                                                                                                                                                                                                                                 tgagagctgtAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAgg          <  2:169241/241‑1 (MQ=255)                                                                                                                                                                                                                                  gagaTCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGc         <  2:89887/250‑1 (MQ=255)                                                                                                                                                                                                                                  gagaTCGTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCg        >  1:238598/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTTAAGCACAGCGAACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGgttgtt  <  1:73777/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACAT‑CTGGTCGAGGATCTCTGCAAGAC‑GTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTT  >  NC_000913/1129552‑1130033

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