Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	15,174	G→A	R336K (AGG→AAG)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,174	0	G	A	97.6%	121.6 / ‑6.4	41	R336K (AGG→AAG)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (18/22);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGT  >  NC_000913/14936‑15410                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             gAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAgct                                                                                                                                                                                                                                  >  2:188439/1‑251 (MQ=255)    cGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTCTGCCGGCGCTGGGTGGCGAAATCCAACTACAGACCCTTGCTGGCCGCGTCCAACTGACAGTGACTGGCGGAACCCCGCCCGGTAAGCTTTTCCGTATGCGCGGTACAGGCGTCGGGTCTGTCCGCGGCGGCGCACCGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAgctgct                                                                                                                                                                                                                               >  1:133428/1‑251 (MQ=255)     gTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGCAATACCGACCCTTGATGGTCGCGTCAAACTGAAAGCGCCTGGCGAAACCCAGACCGGGAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTCGGCCTGAACGCGAAGCAGAAACAgctggtg                                                                                                                                                                                                                              >  1:163263/1‑248 (MQ=255)     gTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGctgctg                                                                                                                                                                                                                              >  1:103773/1‑251 (MQ=255)              acaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctg                                                                                                                                                                                                                     >  1:216412/1‑251 (MQ=255)                aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgca                                                                                                                                                                                                                   <  1:121119/251‑1 (MQ=255)                aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgca                                                                                                                                                                                                                   <  1:134156/251‑1 (MQ=255)                aacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgc                                                                                                                                                                                                                    >  1:63590/1‑250 (MQ=255)                     ctATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaagagctgcaaga                                                                                                                                                                                                                <  2:165362/248‑1 (MQ=255)                            gAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCgg                                                                                                                                                                                                       <  1:145904/251‑1 (MQ=255)                             aaGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGt                                                                                                                                                                                                      <  1:247557/251‑1 (MQ=255)                                   cGATCAGCTTCGCTATGGCGGCGCTCGGTGGCGAAATAGAAGAACCGACCCTTGATGGTCGCGCCCACCTGAAAGTGCCTGGCGAAACCCCGACCGGGAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGGATGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGACGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAAACGCTGCTGCAAAAGCTGCCAGAAAGcgttgggggcccc                                                                                                                                                                                                >  2:38220/1‑239 (MQ=255)                                              gCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGGGGCCCAACCGGGGAGc                                                                                                                                                                                      >  2:69687/1‑250 (MQ=255)                                                 aTGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGGATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCTGTGAGcac                                                                                                                                                                                    <  2:125152/249‑1 (MQ=255)                                                 aTGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGGCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGcaca                                                                                                                                                                                   <  2:71884/250‑1 (MQ=255)                                                       gcgcTGGGTGGCGAAAACGAAGTACCGACCCTTGATTGTCGCGGCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGGCCGCGGTGGCGCACAGGGTGAGTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGccc                                                                                                                                                                            <  2:214015/251‑1 (MQ=255)                                                                 gCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCGTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGTTTCGGTGGTCCAACCGGTGAGCACAACAGCCCGCGCTCAAag                                                                                                                                                                  <  2:124613/251‑1 (MQ=255)                                                                          aaGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAACGCTTCGGTGGCCCAACCGGTGAGCACACCAGCCCGCGCTCAAAGAGCTTCttt                                                                                                                                                         >  1:55629/1‑251 (MQ=255)                                                                               ccGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATgg                                                                                                                                                    <  1:188439/251‑1 (MQ=255)                                                                                       tGATGGTCGCGTCAAACTGAAAGTGCCTGGCTAAACCCAGACCGGTAAGCTGTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACGCCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaa                                                                                                                                              <  2:216412/249‑1 (MQ=255)                                                                                        gATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAATCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTgaagaa                                                                                                                                           <  1:202968/251‑1 (MQ=255)                                                                                               gcgTCAAACTGAAAGTGGCTGGGGAAACCCATACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCATGTCTGTCCGCGGTGGCGCACAGGGTGAGGTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGGTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCGAACCGGTGAGCACAACAGCCTGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGtttttt                                                                                                                                    <  2:16560/251‑1 (MQ=255)                                                                                                     aaCTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTGCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGGCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTgacgac                                                                                                                              <  2:49177/251‑1 (MQ=255)                                                                                                         gAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCATCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTg                                                                                                                           >  1:56761/1‑250 (MQ=255)                                                                                                            aGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGAccc                                                                                                                       <  1:212512/251‑1 (MQ=255)                                                                                                                      gAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGCGATTTGCTGTGCCGCGTTGTCGTCGGAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTcc                                                                                                             >  2:230625/1‑251 (MQ=255)                                                                                                                             aGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCaaa                                                                                                        >  1:98942/1‑249 (MQ=255)                                                                                                                              gACCGGTAAGCTTTTCCTTATTCGCGGTAGAGGCTTCAAGTCTGTCCGCGGTGGCGCACACGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACGCCGGTAGGTCTGAACGAGAAGCAGAAACAGCTGCTGCAGGTGCTGCGAGAAAGCTTCGGTGGCCTAACCGGTGATCACAACAGCCCGCTCGCAGAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGc                                                                                                     <  1:120085/251‑1 (MQ=255)                                                                                                                                     aagcggtgcagTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGGGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGAACTGACCCGCTATCCTCCCCAAAAGCCTGccc                                                                                               <  2:55629/240‑1 (MQ=255)                                                                                                                                          tttCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTGCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTggg                                                                                          <  2:103773/249‑1 (MQ=255)                                                                                                                                                     gcgGTAAAGGCGTCAAGTCTGGCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCGGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTa                                                                              <  1:267052/251‑1 (MQ=255)                                                                                                                                                          aaaGGCGTCAAGTCTGTCCGCGGGGGCGCACAGGGTGATTTGCTGTGGCGCGTTGGCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTTAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGGTTTTTGACGAGCTTACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTaaaat                                                                          <  2:193871/250‑2 (MQ=255)                                                                                                                                                                gTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTg                                                                   >  2:32043/1‑251 (MQ=255)                                                                                                                                                                  caatgcTGTCCGCGTTGGCGCACAGGGTGATTTGCGGTGCCGCGTTGTCGTCGAAACACCGGTAGTCCTGAACGAGAAGCAGAAACAGCTGCGGCAAGAGCTGCAAGAAAGCTTCGGTTGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTGTATGGTGGTAAGAAGGTTTTGGACGACCTGACCCGCTAACCGCCCCAAAAGCCGGCCCGTGGGCAGGCTTGGTTACAATTAGGGTGc                                                                  <  2:243439/245‑1 (MQ=255)                                                                                                                                                                   aaGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGt                                                                >  1:262714/1‑251 (MQ=255)                                                                                                                                                                     gTCTGTCCGAGGTACCGCACAGCGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAGGCAGAAAAAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGGGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCCAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAgggggggttg                                                              >  2:258272/1‑244 (MQ=255)                                                                                                                                                                         gTCCGCGGGGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGATGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGACGTTTTTTGACGACCTGACCCGCTCACCTCCCCCAAAGACAACCCGTGGGCAGGCTTGGGTCCAATTAGGGTGAGTTGa                                                             >  2:149708/1‑248 (MQ=255)                                                                                                                                                                                       acaGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGTAGATATGCGAGCACCTg                                            >  1:62228/1‑251 (MQ=255)                                                                                                                                                                                               aTTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTACAATTAGGGTGCGTTGCAGATATGCGAGCACCTGAAAAGTgg                                    >  2:113751/1‑251 (MQ=255)                                                                                                                                                                                                                             cGGTAGGCCTGAACGAGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa      <  1:69687/251‑1 (MQ=255)                                                                                                                                                                                                                                 aGGCCTGAACGAGAAGCAGCAACAGCTGCTGCCAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGACGTTTTTTGACGACATGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTccccaaaacgctaaattcagggg  >  1:222714/1‑231 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             GAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGT  >  NC_000913/14936‑15410

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