Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,981	G→C	T23S (ACC→AGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,981	0	G	C	98.2%	163.9 / ‑6.1	56	T23S (ACC→AGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (21/34);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTA  >  NC_000913/69741‑70221                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gcCAGCACGTTACCGTCGGCAGCAATCGGTGCGGGCGTCGAGCCGTGACGGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTc                                                                                                                                                                                                                                        <  2:400173/251‑1 (MQ=255)  ccAGAACGTTACCGTCGGCATCAATCGGTTCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCAGGGATTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAACATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATATCTTCACCGCTAGCGAAGtac                                                                                                                                                                                                                                       >  2:354895/1‑249 (MQ=255)            accGTAGGCATCAATAGGTGAGGGCGTCGAGCCGGTCCTGTCAAAGCCAATCCCGACCACAGCGGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGAGTCATCACAAAATTGCACTTTTTGCCAACGGGGATACCACTCGACGCTTGTGGCGATCTCTTCACCGCTTGCGCCGTCCACCGCCaa                                                                                                                                                                                                                               >  2:44446/2‑249 (MQ=255)                        aaTCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAAt                                                                                                                                                                                                                 <  1:334486/251‑1 (MQ=255)                           gggTGCGGGCGACGAGCCGGTCCTGTCAACGCCAATCCCGACCACAGCGGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCCCAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCAc                                                                                                                                                                                                              <  2:169799/250‑1 (MQ=255)                                 gggCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGCACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCa                                                                                                                                                                                                         >  2:199331/1‑250 (MQ=255)                                     gccGAGCCGGTACTGTCAACGCCAAGCTCGACTGCAGCGGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGGTTGCAGGGCCGCTTCAATTGACTCAATGGAGTCACGAGGATGATGACGTAACTGGTTAGTCGGGGCGTCACAAAATTGCCCTTTTTTCAAACGGGGATACCACTCGACGCAGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATc                                                                                                                                                                                                    <  2:167012/249‑1 (MQ=255)                                        gAGCCGGTACTGTCCACGCCACTCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATAGCCCTTTTTGACAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCCGTCCACCGCAAAAGCTCGCACAGAATCACtacccaactcg                                                                                                                                                                                                   >  2:66993/1‑239 (MQ=255)                                         aGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAgg                                                                                                                                                                                                <  1:106201/251‑1 (MQ=255)                                              gtccTGTCAACGCCAATCCCGACCACAGCTTCGCGCTGGTCGACGCTAAGCCCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATTACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATCCCACTCTACGCTGGTGGCTATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCCCAGAATCACGGCCAAAATCGAGGCCAAt                                                                                                                                                                                           <  1:115777/248‑1 (MQ=255)                                                acggTCAACTCCAATCCCGTCCTCATCTGCGCGCTGTTCGACGTTAAGCTCTGCAAGCGCGGTTTTCAGTGCCGCTTCCATTGGCTCAATTGAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAAGTGCCCTTTTTGCCAACGGGGATACCACTCTACGTTGGTGGCGATCTCTTCACCGCTGGCGCAGGCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCTAGGCCaagt                                                                                                                                                                                          <  2:154773/247‑3 (MQ=255)                                                     cAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAAt                                                                                                                                                                                     <  1:697414/250‑1 (MQ=255)                                                         gCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGCCTCAATGTAGTCACGCGGATGATGACGGACCTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAACGCTCGCACAGAATCACTGCCAAAATCGAGGCCACTTGCCATCGcc                                                                                                                                                                                 >  2:176431/1‑250 (MQ=255)                                                           cAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATGCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATc                                                                                                                                                                              <  1:420084/251‑1 (MQ=255)                                                            aaTCCCGACCACAGCTGCGTGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCg                                                                                                                                                                             <  1:121711/251‑1 (MQ=255)                                                               cccGACCACATCTGCGCGCTGTCCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCACCGGGGATACCACGCTACGCGGGTGGCGATCTCTTCACCGCTGGCGTAGTCCACCGTCAAAGCGCGCACAGAATCACTGCCAAAAGCGAGGCCAATTGCAATCGCCATCGttt                                                                                                                                                                          <  2:134483/251‑1 (MQ=255)                                                                ccGACCTCAGCGGCGCGCTGTCCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGGGGCGATCTCTTCACCGCTGGCGCATTCCACCGCCAAAGCTCGCACAGAAGCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTc                                                                                                                                                                         <  2:151853/251‑1 (MQ=255)                                                                  gACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCAc                                                                                                                                                                       >  1:424685/1‑251 (MQ=255)                                                                             ggcgcTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCaaaa                                                                                                                                                            <  1:240288/250‑1 (MQ=255)                                                                              gcgcTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCaaaa                                                                                                                                                            >  1:313525/1‑250 (MQ=255)                                                                               cgcTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCaaaaaa                                                                                                                                                          <  1:706730/251‑1 (MQ=255)                                                                                   gTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGGGCCGCGTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACgg                                                                                                                                                      <  2:249457/251‑1 (MQ=255)                                                                                      cGACGGTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTGTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGGCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTa                                                                                                                                                   <  1:646417/251‑1 (MQ=255)                                                                                      cGACGCTAAGCACTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGt                                                                                                                                                    <  1:486958/250‑1 (MQ=255)                                                                                               gCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGCGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAAc                                                                                                                                          >  2:653138/1‑251 (MQ=255)                                                                                                    gCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTag                                                                                                                                     <  1:489844/251‑1 (MQ=255)                                                                                                    gCAAGCACGGTTTTCAGGGCCGCTTCCATTGTCTCAATGTAGTCACGCGGATGATTACGGAACTGGTTATTCGGGGCAGCACAAAAGTGCCCTTTTTGCCAACGGGGATACCACTCGACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAAACGAGGCCAATTGCAAACGCCATCGTTTCACTCCATCCAAAAAAACGGTTATGGAGAAACAGt                                                                                                                                       <  2:100664/249‑1 (MQ=255)                                                                                                      aaGCACGTTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTagag                                                                                                                                   <  1:471955/251‑1 (MQ=255)                                                                                                                    gTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATaaaa                                                                                                                      >  1:78352/1‑250 (MQ=255)                                                                                                                      tccGCTTCCATTGACTCAAATTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAg                                                                                                                    <  2:552361/249‑1 (MQ=255)                                                                                                                       ccGCTTCCATTGACTCAATGTCGTCACGCGGATGATGACGGAACTGGTTATTCGAGGCATCACAACATTGCCCTTTTTGACACCGGGGATACCACTCTCCGCTGGTGGCGCTCTCATCCCCGCTGGCGCAGTCAACCGCAAACGCTCGCACAGCATCACTGCCAAAATCGAGGCCAATTGCACTCGCCACCGTTTCACTCCATACAAAAAAACGGGTATGGAGCAACAGTAGAGAGTTGCGATAAAAAGCg                                                                                                                  >  1:648022/1‑251 (MQ=255)                                                                                                                                   aCTCAATGTATTCACGCGGATGATGGCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTGTTTGCCAACGGGTATACCACTCTTCGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                                                                                      <  2:30256/251‑1 (MQ=255)                                                                                                                                      cAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCg                                                                                                    <  1:40920/250‑1 (MQ=255)                                                                                                                                                  gcgGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGCTAAAAAGCGTCAGGTAGGATCCGCTAATCTTATgg                                                                                        >  1:386902/1‑250 (MQ=255)                                                                                                                                                    ggATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAt                                                                                      >  1:118025/1‑250 (MQ=255)                                                                                                                                                        gatgaCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTa                                                                                                                                      <  2:277633/198‑1 (MQ=255)                                                                                                                                                        gatgaCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTa                                                                                                                                      >  1:277633/1‑198 (MQ=255)                                                                                                                                                           gaCGGAACTGGTTATTCGGGGCATCACAAAAGTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATg                                                                               <  1:599957/250‑1 (MQ=255)                                                                                                                                                            aCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCt                                                                             >  1:303769/1‑251 (MQ=255)                                                                                                                                                                 aCTGGTTATTCGGGGCATCACAAAATGGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGGCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGc                                                                        <  2:313525/251‑1 (MQ=255)                                                                                                                                                                      tactaaggggCATCACAAAATTGCCCTTTTTGCAAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCCGTCCACCGCCAACGCTCGCACAGAATAACTGCCCAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCTAACAAACGTGTATGGAGAAACAGTAGAGAGTTGCGCTAAAAAGCGTCAGGTCGGTTCAGCTAATCTTATGGCTAAAACTGCTATGGCATAGc                                                                   >  2:576406/7‑251 (MQ=255)                                                                                                                                                                          tCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa                                                                <  1:353139/250‑1 (MQ=255)                                                                                                                                                                          tCGGGGCATCACAAAATTGCCCTTGTTGCCAACGGGGAAACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAAGTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGGTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa                                                                <  2:55103/250‑1 (MQ=255)                                                                                                                                                                           cGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgt                                                              >  1:383132/1‑251 (MQ=255)                                                                                                                                                                             gggCGTCACAAAAGTGCCCTTTTTGCCAACTGGGATACCCCGCTACGCTGGTGGTGATTTCTTCACCGCTGGCGCGGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAAGCGAGGCCAATTGCAAGCGCCCTCGTTTCACTCCATCCAAAAAAACGGGTATGTAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCGAATCGTATGGATAAAAATTCTATGCCATAGGAAAgtgt                                                            <  2:196321/251‑1 (MQ=255)                                                                                                                                                                                cATCTCAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGAACTCGCCACCGCTGGCGCAGTCCACCGCCAAAGTTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGAc                                                         <  1:29195/251‑1 (MQ=255)                                                                                                                                                                                         aTTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCa                                                 >  1:81429/1‑250 (MQ=255)                                                                                                                                                                                                tttCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGCAACAGTAGAGAGTTGCGATAAAAAGCATCAGGTAGGAGACGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGATCAAATAATc                                          >  2:11023/1‑250 (MQ=255)                                                                                                                                                                                                     gCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                                                                                      <  1:596125/185‑1 (MQ=255)                                                                                                                                                                                                     gCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTAGCACAGAATCACTGCCAAAATCGAGGCAAATTGCAATCGCCATCGTTTCACTCAATCCAAAAATACGGGAATGGAGAAACAGTAGAGGGTTGCGATAAAAAGCGTAAGGTAGGATc                                                                                                      >  2:596125/1‑185 (MQ=255)                                                                                                                                                                                                      ccAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTgg                                   >  1:406910/1‑251 (MQ=255)                                                                                                                                                                                                           ggggATACCACTCTCCGCTGGGGGCGATCTCTCCACCGCTGGCGCAGTCCACCGCCAAAGCCCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACttt                              <  2:424685/251‑1 (MQ=255)                                                                                                                                                                                                            gcgATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACtttt                             <  1:624491/249‑1 (MQ=255)                                                                                                                                                                                                             ggATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAACAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATCGCAAAGTGTGACGCCGTGCAAATAATAAATGTGGATTTTTc                            >  2:593484/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggctaTCTCTTCACCGCTGGCGAAGTCCTCCGCCAACGCTCGCACAGAATCACTGCCACAATCGAGGCCAATTGCAATCGCCATCGTATCACTCCATCCAAAAAAACGAGAATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGAAGGATCCGCTAAACTTATGGGTAAAAATGATCTGGCAGAGCAAGGTGTGAAGCCGTGCAAATAATCAATGTGGACTTTTATGACGAGATTAGAGACACtt        >  2:163747/5‑251 (MQ=255)                                                                                                                                                                                                                                         ctTCACCGCTGGCGCATCCCACCGCCAAATCTCGCACAGAATCACTGCCAAAAACGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGTTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGTTATGATCCGCTAATCTTATTGATAAAAATGCTATGGCATAGCAAAGTGTTACGCCGTGCAAATAATCAATGTTGTCTTTTCTGCCGTGATTATAGACTCTTTTGTTa  <  1:231178/249‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTA  >  NC_000913/69741‑70221

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