Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,900	C→T	L71L (CTA→TTA)	def →	peptide deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,900	0	C	T	95.3%	122.1 / ‑2.0	44	L71L (CTA→TTA)	def	peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (20/21);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTG  >  NC_000913/3433663‑3434139                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                aattCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTaatca                                                                                                                                                                                                                                     >  2:127265/4‑250 (MQ=255)   ttCCTGAACACATCTCTGGAGATTCAAGGCAGTTTTGCAAGTGTTCCATATTAAGGCCGAGCGGCTTCGCAACGTTGCTAAACCGGTAGAAGTAGTGAATGCAGAAATTCAGCGTATCGTCGATAATATGTTCGAGAAGATGTTCGAAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGATGATATCCATCAACGTATCATTGTTATTGATATTTCGGAAAACCGTGACGCACGGTTTGTGTTaatcaatc                                                                                                                                                                                                                                  >  2:475487/2‑251 (MQ=255)         aaCACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGc                                                                                                                                                                                                                            >  2:73630/1‑251 (MQ=255)             caTCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCATAAATTCAGCGTATCTTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCtttt                                                                                                                                                                                                                        <  1:19872/251‑1 (MQ=255)                                 aGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGAATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAAc                                                                                                                                                                                                      >  2:392034/1‑249 (MQ=255)                                      tGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGAAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAAAAGGCAtt                                                                                                                                                                                               >  1:273507/1‑251 (MQ=255)                                        cAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTga                                                                                                                                                                                             <  2:165621/251‑1 (MQ=255)                                          aGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGATATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTga                                                                                                                                                                                             >  1:298986/1‑249 (MQ=255)                                                       ccGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTc                                                                                                                                                                              >  1:91978/1‑251 (MQ=255)                                                              aGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCt                                                                                                                                                                       <  1:394273/251‑1 (MQ=255)                                                                   attCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGCTATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCATGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGGGACGAACGGTTAGTGTTAATCAATCCAAAGCTTTTAGAAAAAAGCGGCGCAACAGGCATTGAAGAAGGTTGACTGTCGATCCCTAaac                                                                                                                                                                   >  2:364268/2‑250 (MQ=255)                                                                      cGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACg                                                                                                                                                               >  2:46143/1‑251 (MQ=255)                                                                               gCTAACCCGTTAGAAGAAGTGAATGCAGAAATTCAGCGGATCGTCGATGATATGTGCGAGACGATGTACGCCGAAGGAGGTATTTGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTAGTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTGAATCTATCCAGAGCTGTGATAAAAAAGCGGCGAAACAGGCAGTGAAGAAGGTTGCTTGTCGATCCCTGAACAACGTGCTTTAGt                                                                                                                                                      <  2:96824/251‑1 (MQ=255)                                                                                       ggTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCcgcgcg                                                                                                                                              <  1:475487/251‑1 (MQ=255)                                                                                       ggTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCcgcgcg                                                                                                                                              <  1:197271/251‑1 (MQ=255)                                                                                                agTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGaa                                                                                                                                       >  1:270328/1‑249 (MQ=255)                                                                                                             aTTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTTATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCgcgc                                                                                                                        <  1:279167/251‑1 (MQ=255)                                                                                                                  gCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTg                                                                                                                   <  1:81716/251‑1 (MQ=255)                                                                                                                         gtcgATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACcgcg                                                                                                             >  1:264008/1‑250 (MQ=255)                                                                                                                              tgatATATTCGATACGATGTACGCCGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCAACAAAGAATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGGTAGTGTTAATCAATCCAGACCTTTTAGAACAAAGCGGCTAAACAGGCATTGAAGAAGGTTGCATGTCGATCCCTGAACAACGTGCTTTAGTACAGCGCGCAGAGCAAGTTAAAATTCGCGCCCTTGACCGCGACGGt                                                                                                        >  2:398124/1‑250 (MQ=255)                                                                                                                                  atGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTCTCATTGTTATTGATGTTTCGGAAAACCGTGCCGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAAcc                                                                                                   >  1:246465/1‑251 (MQ=255)                                                                                                                                       cGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTg                                                                                              >  2:272578/1‑251 (MQ=255)                                                                                                                                         agaCGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGa                                                                                             >  1:462370/1‑250 (MQ=255)                                                                                                                                         agaCGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGa                                                                                             >  2:31934/1‑250 (MQ=255)                                                                                                                                                     cAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAgcc                                                                                >  2:227031/1‑249 (MQ=255)                                                                                                                                                     cAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGAc                                                                                <  1:31934/251‑1 (MQ=255)                                                                                                                                                     cAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAATAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACAATTTGAACTGGAAGCAGAc                                                                                >  1:34772/1‑251 (MQ=255)                                                                                                                                                                aTTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGc                                                                     <  2:401374/251‑1 (MQ=255)                                                                                                                                                                        ggcggcAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTg                                                               <  1:355442/249‑1 (MQ=255)                                                                                                                                                                         gcggcAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAAGTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTAt                                                            <  2:172677/251‑1 (MQ=255)                                                                                                                                                                                   aGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATg                                                    >  2:162600/1‑249 (MQ=255)                                                                                                                                                                                       tGATAACCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATGCAGCATGAGATg                                               <  1:392034/250‑1 (MQ=255)                                                                                                                                                                                          tatCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATc                                           >  1:355047/1‑251 (MQ=255)                                                                                                                                                                                             ccATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCAc                                         <  2:462370/250‑1 (MQ=255)                                                                                                                                                                                                 cAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGt                                    <  1:272578/251‑1 (MQ=255)                                                                                                                                                                                                  aaCGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGACGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGCGCACCTGGTc                                   >  1:404703/1‑251 (MQ=255)                                                                                                                                                                                                   aCGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCg                                  <  2:305695/251‑1 (MQ=255)                                                                                                                                                                                                         taTTGTTACTGATGTTTCGGAAACCCGTGACGAACGGTTAGGGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACACCGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGCCGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCTCCTGGTCGGCaaa                             <  1:72264/249‑1 (MQ=255)                                                                                                                                                                                                               tATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGttt                       <  1:9037/250‑1 (MQ=255)                                                                                                                                                                                                                aTTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTAt                     <  1:53286/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaaCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCCGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTcc        >  2:139196/1‑250 (MQ=255)                                                                                                                                                                                                                             aaaaCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCa        <  2:148295/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaaCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCGTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCa        <  2:91978/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTg  >  2:46761/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                AACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTG  >  NC_000913/3433663‑3434139

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