Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,086,634	G→A	A143A (GCC→GCT)	plaP ←	putrescine:H(+) symporter PlaP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,086,634	0	G	A	97.8%	142.6 / NA	45	A143A (GCC→GCT)	plaP	putrescine:H(+) symporter PlaP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/24);  minor base C (1/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAAACAGG  >  NC_000913/2086389‑2086877                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    gccgtcaggtggcgtaTGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGACGCCCATGTAAGCATCGTGATAACATGCTCATCGCCAGACCATAATGGACGAGTGCTCGCCAACGTACCAGCGCCGTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCGGCACTACCTGCAACACGACGATTACAGTGTTGAATTTCGCTACGGATGTCAGACTACGCAGGTTAAAAGCgtt                                                                                                                                                                                                                                                <  2:549813/236‑3 (MQ=255)       aGGTGGCTGATGCCGTCAAAGACGGTAAAGGAGAAGCACATGAGTGTCGCCCCGGTAAGCAGCGGGATAACATGCGCATCGCCAGACCAGACTAGACGAGTGCTCGCCCGCGTACAAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCGGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATaaa                                                                                                                                                                                                                                          <  2:448317/251‑1 (MQ=255)              tGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTACAAGCGGTCCTAAAGGccacc                                                                                                                                                                                                                                   >  1:104020/1‑251 (MQ=255)                       cAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGccacc                                                                                                                                                                                                                          <  1:254017/251‑1 (MQ=255)                            ccGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACaa                                                                                                                                                                                                                      >  1:232181/1‑250 (MQ=255)                               gTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACAt                                                                                                                                                                                                                  <  1:170450/251‑1 (MQ=255)                                    ggcgaAGCACAGGATTGTCGCCCCGGTAATCAACTAGAAAACATGCGCATCGCCAGCCCAGAATGAACGTGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCAATAAGCAACCATGCCCCGAATCACCGCGCTCAGCACTACCTGCAACACGACGATTACGGTGTTGAAGATCGCTACGGATTTCAGACTACGCAGGTTAACCGCGGTCATAAAGCCCACCAGCCCCACCACAAACATccccg                                                                                                                                                                                                             >  2:44493/4‑248 (MQ=255)                                    ggAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACg                                                                                                                                                                                                             >  1:285125/1‑251 (MQ=255)                                                 aTTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCATCACAAACATCCACGATGGGATGGAAgg                                                                                                                                                                                                >  2:379739/1‑251 (MQ=255)                                                       gccacGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCag                                                                                                                                                                                          <  1:452480/247‑1 (MQ=255)                                                        ccccGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCaga                                                                                                                                                                                         <  1:440813/251‑1 (MQ=255)                                                               aaTCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCaa                                                                                                                                                                                  <  1:466131/251‑1 (MQ=255)                                                                atcatcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCaaa                                                                                                                                                                                 <  1:152466/251‑1 (MQ=255)                                                                  caacGGGATAAAATGCGCAGCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTACAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTGCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGATCTTCaaaa                                                                                                                                                                                <  2:95515/247‑1 (MQ=255)                                                                   atcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCATACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAAt                                                                                                                                                                               >  1:131546/1‑250 (MQ=255)                                                                     cGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGa                                                                                                                                                                            <  1:379739/251‑1 (MQ=255)                                                                                     gcgcaagaccaaaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTTCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGTAAGGCACCAGATCTTCAAAATAGATTTTCGCCAGTAGa                                                                                                                                                              <  2:373064/244‑1 (MQ=255)                                                                                            accagaATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGatc                                                                                                                                                     >  1:385339/1‑251 (MQ=255)                                                                                                  aTGGCCGAGTGCTCGCCAGCTTACCAGCGCCTGCAACTTCAAAAACGCCATAAACAACCATGCCCAGAATCTCCGCGATCAGCACTACCTGCAACACGACGATTACAGTGGTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATTGAAGGCACCAGAGCTGCAAAATAGTTTTTCGCCAGTATAATGTTGATGATCGgc                                                                                                                                               <  2:104020/251‑1 (MQ=255)                                                                                                                                cttcaccttcaACTACGCAATAAACAACCATGCCCAGAATCACCTCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGACGTTCGCTACGGATTTCAGACAACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACAACAAACATCCACGATGGGATGGAAGCAACCAGAGCTTCAAAACAGATTTTCGCCAGTATAATGTTGGACATCGGCGCGACAAGATAGTCAAGCAGAGAATACCa                                                                                                                  >  2:272355/1‑250 (MQ=255)                                                                                                                                    accttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGCAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGCAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTAGAGCAGCGAAGACCTAcccc                                                                                                             >  1:564742/1‑250 (MQ=255)                                                                                                                                      cttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccacc                                                                                                           >  1:222756/1‑251 (MQ=255)                                                                                                                                               cGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCg                                                                                                  >  1:178795/1‑251 (MQ=255)                                                                                                                                                    tAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcgc                                                                                                                                             <  1:144415/203‑1 (MQ=255)                                                                                                                                                     aaaCAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTc                                                                                            <  1:533947/251‑1 (MQ=255)                                                                                                                                                         aacCATACCTAGAATCACCGCGATCAGCACTACCTGCAACACGACGCTTACAGTGTAGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGGTGGAAGGCACCAGAGCATCACCATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGGCCAACCCACCATAAAGCCGCCCGTCGGGc                                                                                        >  2:214583/1‑251 (MQ=255)                                                                                                                                                                cccAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCATCATAAAg                                                                                                     >  2:393826/1‑231 (MQ=255)                                                                                                                                                                cccAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAg                                                                                                     <  1:393826/231‑1 (MQ=255)                                                                                                                                                                        cACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAGAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCATAGAAGACCAACCCACCATAAAGCCGACAGTCGGGc                                                                                        <  2:522240/236‑1 (MQ=255)                                                                                                                                                                        cACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAGAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGc                                                                                        >  1:522240/1‑236 (MQ=255)                                                                                                                                                                        cACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTggg                                                                         >  1:388413/1‑251 (MQ=255)                                                                                                                                                                              gATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTa                                                                     >  1:326719/1‑249 (MQ=255)                                                                                                                                                                                     ccTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGc                                                            <  1:351986/250‑1 (MQ=255)                                                                                                                                                                                        aCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCa                                                           <  1:485988/249‑1 (MQ=255)                                                                                                                                                                                         ccTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCagag                                                        <  1:176268/251‑1 (MQ=255)                                                                                                                                                                                            gCAACACGACGATTACAGTGGTGAAGTTCGCTACGGAATTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGcc                                                      <  1:226203/250‑1 (MQ=255)                                                                                                                                                                                                   gacgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGa                                                >  1:26175/1‑249 (MQ=255)                                                                                                                                                                                                    acgaTTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAgg                                             <  2:232181/251‑1 (MQ=255)                                                                                                                                                                                                                        cGCTCCGGATGTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTcc                         <  2:159711/251‑1 (MQ=255)                                                                                                                                                                                                                         gCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGCTGGAAGGCACCAGAGCTTAAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGACGCGAACAGATAATCGAGCAGAGAAGACCAACCCACAATAAAGCCGACAGACGGGCTAATGGATTTCTGGGCGTAGGGGTATACAGAGCCAACCGAAGGGTAGCGGAGAACCAGCTTccc                        >  2:159459/1‑251 (MQ=255)                                                                                                                                                                                                                                    cAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATAATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCag               >  1:196782/1‑249 (MQ=255)                                                                                                                                                                                                                                      gACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGcc           <  1:323277/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTaa       >  1:239426/1‑250 (MQ=255)                                                                                                                                                                                                                                              gtggtAAAAGCGGTCATAAAGTCCACCAGCGCAACCCCAAAGATCCACGTTGGGATGGAAGTCACCAGAGCTTCAAAATAGATGTTCGCCAGTAGAATGTTGAGCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCTACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGATCCGTAAACAg   <  2:285125/247‑1 (MQ=255)                                                                                                                                                                                                                                                gTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAAACAgg  >  1:164778/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    TCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAAACAGG  >  NC_000913/2086389‑2086877

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