Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,072,484	C→T	G127D (GGC→GAC)	rutB ←	ureidoacrylate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,072,484	0	C	T	98.0%	165.4 / NA	51	G127D (GGC→GAC)	rutB	ureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/23);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGAC  >  NC_000913/1072241‑1072719                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            aaaaaGGTTTCGATATTGAACAACGCGGCTTTCTGCGAAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGAGGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTg                                                                                                                                                                                                                                      >  1:230036/1‑251 (MQ=255)    aaGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCgg                                                                                                                                                                                                                                   >  1:511461/1‑251 (MQ=255)    aaGGTTTCGATATTGAACAACGAGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCgg                                                                                                                                                                                                                                   >  1:207505/1‑251 (MQ=255)         ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCAcc                                                                                                                                                                                                                               >  1:157938/1‑250 (MQ=255)         ttCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCAcc                                                                                                                                                                                                                               <  1:383999/250‑1 (MQ=255)         ttCGATATTGAACAACGCGGCTTTCTGCGAAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCAcc                                                                                                                                                                                                                               <  2:248666/250‑1 (MQ=255)            gATATTGAACCACGAGGCTTTCTGCGCAAATTTCGGCCGCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAACCCGCTGTAGAGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGtt                                                                                                                                                                                                                           >  2:288813/1‑251 (MQ=255)            gATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGtt                                                                                                                                                                                                                           >  1:394347/1‑251 (MQ=255)             atatTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTc                                                                                                                                                                                                                          <  2:180900/251‑1 (MQ=255)                                                      cTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAACGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAAtt                                                                                                                                                                                 >  2:95186/1‑251 (MQ=255)                                                        ggTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAAtt                                                                                                                                                                                 <  2:651870/249‑1 (MQ=255)                                                                        aaGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGgct                                                                                                                                                               >  1:197833/1‑251 (MQ=255)                                                                            cccacGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCtt                                                                                                                                                           <  1:142339/250‑1 (MQ=255)                                                                                  ccGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGCCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCAt                                                                                                                                                     <  2:215874/251‑1 (MQ=255)                                                                                     aaaTACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCCTGGt                                                                                                                                                  >  2:38184/1‑251 (MQ=255)                                                                                              aGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGc                                                                                                                                         <  1:409825/251‑1 (MQ=255)                                                                                              aGAAAAAAGCCGTCGCGAAGCGTCGATTCGAAGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGTTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAggg                                                                                                                                          >  1:419344/1‑250 (MQ=255)                                                                                                 aaaaaGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGtt                                                                                                                                      <  2:43574/251‑1 (MQ=255)                                                                                                                aGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAAtt                                                                                                                       >  2:131877/1‑251 (MQ=255)                                                                                                                   gTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTAATGAAAATTCgg                                                                                                                    >  1:420309/1‑251 (MQ=255)                                                                                                                     cGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTg                                                                                                                  >  2:56396/1‑251 (MQ=255)                                                                                                                           gACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCgg                                                                                                            >  1:577768/1‑251 (MQ=255)                                                                                                                           gACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCAC‑AGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCggg                                                                                                           <  1:589935/251‑1 (MQ=255)                                                                                                                              gCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTc                                                                                                         <  1:78541/251‑1 (MQ=255)                                                                                                                               cAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGt                                                                                                                   <  1:566695/240‑1 (MQ=255)                                                                                                                               cAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGt                                                                                                                   >  2:566695/1‑240 (MQ=255)                                                                                                                                    gTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGcc                                                                                                     >  2:575840/1‑249 (MQ=255)                                                                                                                                      tGGTAGCGATGCCGGTGAAAACCAGATGGCGTATGCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCAGGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGc                                                                                                  <  1:95186/250‑1 (MQ=255)                                                                                                                                        gTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAACCCGCTGTAGCGCGGCTGCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGTCGCCTTTCGCCAGCAATTTCCCCTACAGCTGCGGCTGCTTACGCATGGTTTTCAGTGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTc                                                                                               <  2:18087/251‑1 (MQ=255)                                                                                                                                           gCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGAc                                                                                            >  1:421683/1‑251 (MQ=255)                                                                                                                                                   ggTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGtt                                                                                    <  2:498375/251‑1 (MQ=255)                                                                                                                                                       aacaCCAGATGGCGTATTCCGCGGCTTGGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCCCGGCTTCGGCAGCACAATATCGGCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCGTTCGCCAGCAATTTCCCCGGCAGCTGCTGCTTCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGTGTCTACATCCTGTTCATc                                                                                <  2:419344/248‑1 (MQ=255)                                                                                                                                                            cAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGc                                                                           <  1:121267/251‑1 (MQ=255)                                                                                                                                                              gATGGCGTATTCAGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGATGTAGCGCGGATTCGGCAGAACAATATCGTAAGGCTGCGGGAACAGTTCATCAACCAGTTGATAATCCAAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGATTACGCATGGATTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGACCGCCAGCCTAGACATACTGGTCATACCAGcc                                                                          >  2:368528/1‑250 (MQ=255)                                                                                                                                                                                 gcgcAAAATGCTGTCCAGCGGCGTTTTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCGTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCa                                                      <  2:511461/251‑1 (MQ=255)                                                                                                                                                                                  cgcAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAg                                                     <  1:339478/251‑1 (MQ=255)                                                                                                                                                                                     aaaaTGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCAt                                                  >  1:508348/1‑251 (MQ=255)                                                                                                                                                                                         ggCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCt                                              <  2:71781/250‑1 (MQ=255)                                                                                                                                                                                                       gTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGc                                >  1:85441/1‑251 (MQ=255)                                                                                                                                                                                                       gTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGc                                >  2:468034/1‑251 (MQ=255)                                                                                                                                                                                                           tGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGt                             <  1:487453/250‑1 (MQ=255)                                                                                                                                                                                                                aaaCCGCGGTAGCGCGGCTTCGGCAGCACAAAATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTTAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCAcggc                       <  1:591934/251‑1 (MQ=255)                                                                                                                                                                                                                   ccGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGt                    >  1:69469/1‑251 (MQ=255)                                                                                                                                                                                                                    cGCTGTAGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGtt                   >  2:240220/1‑251 (MQ=255)                                                                                                                                                                                                                          aGCGCGGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATg             >  1:468851/1‑251 (MQ=255)                                                                                                                                                                                                                             gcgGCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTg          <  1:247055/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTATGTTGGc        >  1:571800/1‑251 (MQ=255)                                                                                                                                                                                                                                gCTTCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGACGATCAGCATCCCTGCCGCTCGCGCTGCGGTCCCGGCGGTTTGAATTTTgtg        >  2:284345/1‑248 (MQ=255)                                                                                                                                                                                                                                   tCGGCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCCCCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCaa      <  1:575840/249‑1 (MQ=255)                                                                                                                                                                                                                                   ccggcggCACAATCCCGCCAGCCTGCCTCACCTGTTCACCCACCAGTCGCTAATCCCAGGCGCCTTTCGCCAGCAATTTCCCCTTCTGCTGCGCCTGCTTACCCATGGTTTTCAGGCGGTTCGATTGCTGAAAATTCGCTCAGCCGGGTCCGCCAGTCTCTACATACTGTCCCTCCCAGCCATTTTGAACCCCGTTGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATg    <  2:303134/250‑1 (MQ=255)                                                                                                                                                                                                                                     ggCAGCACAATATCGTCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGAc  <  1:58709/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            AAAAAGGTTTCGATATTGAACAACGCGGCTTTCTGCGCAAATTTCGGCCCCGCCTGGTGAGTTGCGTCTTCAAGCACCACGCCGAAATACTCCAGAAAAAAGCCGTCGCGTAGCGTCGATTCGACGCAGACGTTGGTAGCGATGCCGGTGAAAACCAGATGGCGTATTCCGCGGCTGCGCAAAATGCTGTCCAGCGGCGTATTGAAGAAACCGCTGTAGCGCGGCTTCGGCAGCACAATATCGCCAGGCTGCGGCACCAGTTCATCCACCAGTTGATAATCCCAGGAGCCTTTCGCCAGCAATTTCCCCTGCAGCTGCGGCTGCTTACGCATGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGAC  >  NC_000913/1072241‑1072719

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