Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,935	C→T	Q38Q (CAG→CAA)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,935	0	C	T	98.3%	187.8 / NA	60	Q38Q (CAG→CAA)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/35);  minor base G (0/1);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATC  >  NC_000913/69688‑70181                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGACAACGCCAATCCAGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCCTTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGc                                                                                                                                                                                                                                                     >  2:599957/1‑251 (MQ=255)     ccctcatacgagcgtggtgTTCGGGTTCTCGGCATGCGCCGGCCGCAGCACCAGCCCGTTCCCGTCGGCATCAACCGGTGCGGGCGTCGAGCCGGTTCTGTCAACGCCAATCCCGACCCCAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTTACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATGCGGGGCAGCGCAAAATTGCCCTTGTTGCCAAc                                                                                                                                                                                                                                                 <  2:486132/233‑1 (MQ=255)      cacaATACGAACATCGCGTTCGGGTTTTCGGCAAACCCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACg                                                                                                                                                                                                                                                >  1:589552/1‑251 (MQ=255)         aaTACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACgggg                                                                                                                                                                                                                                             <  1:158367/251‑1 (MQ=255)             cGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGGACTGGGTATTCGGGGAATCACAGAATTGCCCTTTTTGCCAACGGGGATAc                                                                                                                                                                                                                                         >  2:489844/1‑251 (MQ=255)                                          cctgggCGCGGCGCCGGCACGGTACCGTCGGCGTCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCGCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTGCGCTGGTGGCGATCTCTTCAcc                                                                                                                                                                                                            <  2:223932/248‑1 (MQ=255)                                           ccGGGCGCAGCGCCAGCACGGTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCAcc                                                                                                                                                                                                            <  2:386380/250‑1 (MQ=255)                                                cgcAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTAGCAAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCgatggc                                                                                                                                                                                                      >  1:516151/1‑246 (MQ=255)                                                   aGCGCCAGCACGTTACCGCCGGCATCAATACGTGCGGACATCGCGCCGGTCCTGTCAACCCCCATCCCACCCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGACGCTTCAATTGACTAAATGTAGTCACGCGGATGATGACGGACCTGGTTATTCGGGGCATCACAAAATTGCCCTCTTTGCAAACGGGGATACAACTCTACTCTGGTGGCGATCTCTTCAccactggcgcc                                                                                                                                                                                                   >  2:468087/1‑242 (MQ=255)                                                      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)                                                                                                                                   gcgcTGTTCGACGCTAAGCTCTGCAAGCACGGTTTGCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGcc                                                                                                                                                                                                                                                    <  1:539474/122‑1 (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)                                                                                                                                                                                   cATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTc                                                                                                                                                                                                               <  1:330628/111‑1 (MQ=255)                                                                                                                                                                                   cATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTc                                                                                                                                                                                                               >  2:330628/1‑111 (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)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATC  >  NC_000913/69688‑70181

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