Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,092	C→T	E93E (GAG→GAA)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,389,092	0	C	T	100.0%	137.7 / NA	45	E93E (GAG→GAA)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/25);  total (20/25)

TGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCC  >  NC_000913/4388856‑4389325                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          tgtgCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTACCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCAtgt                                                                                                                                                                                                                              <  1:97323/250‑1 (MQ=255)  ggtactgattttggggaggaTTGCCGGTAAGCGTGCCCAGGGTGCGCTTGATCTCATTTAACTGACGGCGGTCGTCGTTTTGCTGTTGCGGCAGTTGAGTCAGCGAATCGGCAAGCTCGCGGGCGCGCTCCGGCTCTTGCTGGGGCTGTCGGCTTGTATCCAGCAACTGGCTGCTGACCTGGATAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCTACACGCTGCTCGTTTCGTCACGCAtgtt                                                                                                                                                                                                                             <  2:358896/236‑1 (MQ=255)             aGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAACGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGCCGGCTTTTATCCCACAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCCACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATtgtgg                                                                                                                                                                                                                 >  1:377730/1‑250 (MQ=255)              gCGGAGTATGTCCGGTAAGCGGTCCCAGGCGGGGCTCGATCTCATTTAACGGACGGCGGGCGGCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAAGCGGCAATCTCGCGGGCGCGCTCCTGCGCTTGCTGGGCCTGACGGCTTTTATCCAGTAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTgcg                                                                                                                                                                                                                <  2:344814/251‑1 (MQ=255)                   gTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCag                                                                                                                                                                                                           <  2:72579/251‑1 (MQ=255)                       tGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGGTTTCGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTTTCCAGCAACTGGCTGCTGACCTGGAGAATTTTCTGATTCTGCGCGGCGGTAGACATTCCCGGCGGCACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTTTGCGCGCagaggg                                                                                                                                                                                                       >  2:207114/1‑249 (MQ=255)                        gCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAg                                                                                                                                                                                                      <  1:359762/251‑1 (MQ=255)                                             gcgcTCGATCTCATTTAACTGACGGCGGGCGGCGGTTTGCTGTCGCGGCAGGTGATTCAGCGAATCGGCAATCTCGCGGGCGCGTTCCTGCTCGTGCTGGGCCTGACGGCGTTTATCCGGCAACTGGCTGCTGGCCTTGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCTCGGCGACACTCTGCGCGGTTGGTCACGCATTGTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTAt                                                                                                                                                                                 <  2:243625/251‑1 (MQ=255)                                               ccTCGATCTCATTTAACTGACGGCGGGCTCCGGTTTGCTCTTGCCGCAGTTGATTCCGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCg                                                                                                                                                                               <  1:80790/250‑1 (MQ=255)                                                    ttctcATTTAACTGGCGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCGCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTAt                                                                                                                                                                                 <  2:242003/243‑1 (MQ=255)                                                    aTCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGGTTCAGCGCATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGGTGCTGACCTGGGGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCCGCGCGGTTCCGCACGCCTGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGGGTTTCTGGTAATTAt                                                                                                                                                                                 >  1:242003/1‑244 (MQ=255)                                                               cTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTTCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGGCACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCTGAGAGTTTCGGATAATTATCGATAACTTGCTGATAtt                                                                                                                                                               >  1:402823/1‑251 (MQ=255)                                                                 gactgctggcgTCTGTTTGCTGTTGCCGCAGTTGATTCATCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCCGAGTATCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGt                                                                                                                                                             <  2:409953/244‑1 (MQ=255)                                                                       gggcgGCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATg                                                                                                                                                       <  1:390799/251‑1 (MQ=255)                                                                           gTCGGTTTGCTGCTGCGGCAGTTGATTCAGCGAATCGGCAATCTAGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACAAGATGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGATAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATACGCt                                                                                                                                                   >  1:333881/1‑251 (MQ=255)                                                                              ggTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGTAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCa                                                                                                                                                 >  1:299533/1‑250 (MQ=255)                                                                                              aGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTc                                                                                                                                <  1:236855/251‑1 (MQ=255)                                                                                              aGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGATTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGCTTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTAGTCACGCATGTTGTTTAATTGTGAGCGCCGAGTAGCGGCGAGTTTCGGATACTTATCGATAACTTGCTGATATTGTTGGAGGCGCTCAAGGGAACCTtttctttc                                                                                                                                >  2:124003/1‑251 (MQ=255)                                                                                                     tCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAg                                                                                                                          <  1:352093/250‑1 (MQ=255)                                                                                                               ggCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCTCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAgg                                                                                                               <  1:346217/251‑1 (MQ=255)                                                                                                                       cgcgGTCGCGCTCCTGCTTTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCTGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTgg                                                                                                       <  2:402823/251‑1 (MQ=255)                                                                                                                       cgcgGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTgg                                                                                                       >  2:412742/1‑251 (MQ=255)                                                                                                                          gggcttgcTCCGGCTCTTGCTGGGCCTGTCGGCTTTTATCCAGCAACTGGCGGCTGACCTGGAGAATTGCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCGGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCAGTTAAGGCAGACTGGAgc                                                                                                    <  2:180327/245‑1 (MQ=255)                                                                                                                               cgcTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGCTAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCctc                                                                                               >  2:345628/1‑251 (MQ=255)                                                                                                                               cgcTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCctc                                                                                               <  1:207114/251‑1 (MQ=255)                                                                                                                                               gggCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCt                                                                                 <  1:412742/249‑1 (MQ=255)                                                                                                                                                   cTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGACTTTCCTGATTCAGCGCGTCGGTAGACATTACCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGGGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCACGGGAACCTTTTCGTTCCTCAAGCGCCTTTAAGGCAGACTGGAGAGCCTCTACGACTTCCGGCTGTGCgg                                                                           >  2:180912/1‑251 (MQ=255)                                                                                                                                                    tGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCggg                                                                          >  2:3925/1‑251 (MQ=255)                                                                                                                                                     gACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGCGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGt                                                                         >  1:388035/1‑251 (MQ=255)                                                                                                                                                            tttATCCAGCATCTGTCTGCTGACCTGGCAAATTTCCTGAGTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGTTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGGGTAGCGGTGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGTGTTTcgcc                                                                   <  2:188650/250‑1 (MQ=255)                                                                                                                                                                       aCTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCACGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCCTGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTTAGGCCGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGATTTTGCCTgt                                                       >  1:148003/1‑250 (MQ=255)                                                                                                                                                                        cTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCttt                                                                                                                                      >  2:21932/1‑171 (MQ=255)                                                                                                                                                                        cTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCttt                                                                                                                                      <  1:21932/171‑1 (MQ=255)                                                                                                                                                                         tGGCTGTTGACCTGGCGAATTTCCTGCTTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGCTAAATTGCTGATATTGTTTGATGCGCCCAAGGGAACCTTTTCGTTCCTCTAGCGCATTTAAGGCAGCCTGGAGCGCCTCTAAGACTTCCGGCTGTGGGGGTTTCGCCGCTTTTGCCGGCTc                                                     >  2:375246/1‑251 (MQ=255)                                                                                                                                                                         tGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGGGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCATCTACGCCTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTc                                                     >  2:375254/1‑251 (MQ=255)                                                                                                                                                                                  aCCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGGTTAGGATAATTATCGATAACTTGCTGATATTGTTTGCTGCGCTCAAGGGAACCTTTTCGTTCCTCACGCGCATTTAAGGCAGACTGGAGCGCCTCTACGGCTTCCGGCTGGGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTg                                            >  2:370221/1‑251 (MQ=255)                                                                                                                                                                                   ccTGGAGAATTTCCTGATTCAGCGCGGCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGTAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGa                                           <  2:299533/251‑1 (MQ=255)                                                                                                                                                                                              gccTGATTCAGCGCGTCGGTAGACATGCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAAGTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGAGAACTTGCTGATATTGGTTGAGGCGCTCAAGGGAACCTTTGCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGAgttgttt                                <  2:388035/250‑1 (MQ=255)                                                                                                                                                                                               ccTGATTCAGCGCGTCGGTAGACATGCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGAtttgttt                                <  1:375254/250‑1 (MQ=255)                                                                                                                                                                                               ccTGATACAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTGTTCGTGCCTCAAGCGCATTGAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGAtttgtgt                                <  1:375246/250‑3 (MQ=255)                                                                                                                                                                                                                gTAGACATTCCCGGCGACACGCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTcgc              >  1:240673/1‑251 (MQ=255)                                                                                                                                                                                                                        tCCCGGCGACACGCTGCGCGGTTCGTCACGCAGGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTAc       <  2:176764/250‑1 (MQ=255)                                                                                                                                                                                                                            ggCGACACGCTGCGAGGTTCGTCATGCATGTTGTTTAATTGTGCGCGCAGCGTCGCGGAGCGTTTCGGCTACTTATAGATAACTTGCTGATATTGTTTGATGCGCTCACGGGAACCTTTTCGTTCCTCACGCGCATTTAACGCAGACTGAAGCGCCACTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGCGTGATTTGTTTGCTATCGGGGGCCGTCCCGGCGTACGcccc  >  2:414582/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTc                                                                      >  1:113785/1‑175 (MQ=255)                                                                                                                                                                                                                                    gCTGCGCGGTTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTc                                                                      <  2:113785/175‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          TGTGCCTGATTGAGCGGAGTATTGCCGGTAAGCGTTCCCAGGCGGCGCTCGATCTCATTTAACTGACGGCGGGCGTCGGTTTGCTGTTGCGGCAGTTGATTCAGCGAATCGGCAATCTCGCGGGCGCGCTCCTGCTCTTGCTGGGCCTGACGGCTTTTATCCAGCAACTGGCTGCTGACCTGGAGAATTTCCTGATTCAGCGCGTCGGTAGACATTCCCGGCGACACGCTGCGCGGCTCGTCACGCATGTTGTTTAATTGTGCGCGCAGAGTAGCGGAGAGTTTCGGATAATTATCGATAACTTGCTGATATTGTTTGATGCGCTCAAGGGAACCTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCC  >  NC_000913/4388856‑4389325

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