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	100.0%	131.5 / NA	40	S12F (TCC→TTC)	flgN	flagellar biosynthesis protein FlgN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/18);  total (22/18)

CCGTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCC  >  NC_000913/1129550‑1130039                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ccGTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACAt                                                                                                                                                                                                                                                 >  2:271562/1‑251 (MQ=255)   gTTAGGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCt                                                                                                                                                                                                                                               <  1:471040/251‑1 (MQ=255)           gATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAgg                                                                                                                                                                                                                                       >  1:133110/1‑251 (MQ=255)                    ggCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGc                                                                                                                                                                                                                               >  1:87051/1‑250 (MQ=255)                      cGTAGTTGCTGCGTTTTCACAGTAATTGTCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTGGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAg                                                                                                                                                                                                                            <  2:330306/251‑1 (MQ=255)                         aGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGCGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGACCATCTGGTCGAGGATCTCTGCAAGACg                                                                                                                                                                                                                         >  2:269486/1‑251 (MQ=255)                            tgctgcGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATAGCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACgtg                                                                                                                                                                                                                       <  1:453029/250‑1 (MQ=255)                                          aGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTa                                                                                                                                                                                                        <  1:106193/251‑1 (MQ=255)                                                 tCCTGCCAGCGTTGACTAATAGCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGGAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTGCTCTGc                                                                                                                                                                                                  <  2:383039/250‑1 (MQ=255)                                                   cggCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCGTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCTTTAAGCACAGCGAACATCTTGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAg                                                                                                                                                                                               <  2:471045/249‑1 (MQ=255)                                                                  aaTATCAACTCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCg                                                                                                                                                                                <  1:234941/251‑1 (MQ=255)                                                                   atatCAACGCTATGTGCTGTATTGGGTTCTGTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGt                                                                                                                                                                               <  2:124781/251‑1 (MQ=255)                                                                          cGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTGTGCAAGTCTTGCTGCGCTTCGTTGATCAg                                                                                                                                                                        <  1:417204/251‑1 (MQ=255)                                                                                   cTGTATTGGGTTCTTTCCTGCGTAACTGCTAGCGGTAACCCAGCGTCGCCCGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATCGAGAGATGTTGCTGCTCTTGATCCATTACCGCTTTGAGATCGTTAAGCACCGCGAACATCTGGTCGAGGATCTCTGCACGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCCTCGGc                                                                                                                                                               >  2:225830/1‑251 (MQ=255)                                                                                   cTGTATTGGGTTCTTTCATGCGTAACTGATCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATAAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGc                                                                                                                                                               >  1:87523/1‑251 (MQ=255)                                                                                        ttGGGTTCTTTCCTGCGTAACAGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGCCGAGGATCTCTGCAAGACGTGTCATACGATAAGTTACTCTGCCAGTCTTGCTGCGCTTCGTTTATCAGCGCATCGGCAAttt                                                                                                                                                          >  2:411043/1‑251 (MQ=255)                                                                                          gggTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTg                                                                                                                                                        <  1:155795/251‑1 (MQ=255)                                                                                                  tCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGgtgt                                                                                                                                                 <  1:376451/250‑1 (MQ=255)                                                                                                   ccTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTc                                                                                                                                                <  2:476718/250‑1 (MQ=255)                                                                                                          aaCTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAg                                                                                                                                                                                                           <  1:262943/184‑1 (MQ=255)                                                                                                          aaCTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAg                                                                                                                                                                                                           >  2:262943/1‑184 (MQ=255)                                                                                                                  gAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCg                                                                                                                                <  1:271562/251‑1 (MQ=255)                                                                                                                    ggTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGtt                                                                                                                              <  1:370027/251‑1 (MQ=255)                                                                                                                          ccAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCCGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAAt                                                                                                                        >  1:252275/1‑251 (MQ=255)                                                                                                                          ccAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAAt                                                                                                                        >  1:137963/1‑251 (MQ=255)                                                                                                                                            aGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTc                                                                                                      >  2:364831/1‑251 (MQ=255)                                                                                                                                               ttttttGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATc                                                                                                                                                                                                                                                <  2:59390/110‑1 (MQ=255)                                                                                                                                               ttttttGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATc                                                                                                                                                                                                                                                >  1:59390/1‑110 (MQ=255)                                                                                                                                               ttttttGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGac                                                                                                   >  1:260603/1‑251 (MQ=255)                                                                                                                                                     gTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGGTTCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTATCGCTTCGACACGtt                                                                                              <  1:364110/250‑1 (MQ=255)                                                                                                                                                                    gCAACTGGCTGCCGATGATCTGCCCCATAGAGCGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGCCAAGTTCAGGATTGATATCAct                                                                              >  2:232678/1‑251 (MQ=255)                                                                                                                                                                                    gATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGAGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATc                                                              >  2:200786/1‑251 (MQ=255)                                                                                                                                                                                           cccATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTttg                                                       >  2:312335/1‑251 (MQ=255)                                                                                                                                                                                                agagagATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACGCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTttgcttg                                                   <  2:260603/250‑1 (MQ=255)                                                                                                                                                                                                              gctCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGt                                    >  1:328184/1‑251 (MQ=255)                                                                                                                                                                                                                   tGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCccac                               >  2:251379/1‑251 (MQ=255)                                                                                                                                                                                                                                     agaTCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGCTCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCgg             >  2:186632/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGgttt          >  2:412604/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGttttt     >  1:56491/1‑251 (MQ=255)                                                                                                                                                                                                                                                cacaGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGttttttcc  >  1:46286/1‑248 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CCGTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCC  >  NC_000913/1129550‑1130039

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