Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,517,171	G→A	L46F (CTT→TTT)	damX ←	cell division protein DamX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,517,171	0	G	A	97.8%	135.9 / NA	45	L46F (CTT→TTT)	damX	cell division protein DamX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/24);  minor base T (0/1);  total (20/25)
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

GCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCT  >  NC_000913/3516939‑3517404                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          gcAGGACCGATACCGATGATCAACAGCAGTAGATCCAGGATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGTCGCTCTACGCGCTCTTCATCCACGGTTTCATTTTCGGATGCGTCAATTGCTTCTTCGAGGTCCTGTTTCTCATTGCGCTCTTTTTGCTTACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAtt                                                                                                                                                                                                                         <  1:196648/251‑1 (MQ=255)    gAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGAt                                                                                                                                                                                                                      >  1:316217/1‑251 (MQ=255)      accgataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCt                                                                                                                                                                                                                    <  1:164423/251‑1 (MQ=255)      accgataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCAATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAGATTCAATATCATCAAAATTGATCt                                                                                                                                                                                                                    >  1:294830/1‑251 (MQ=255)                atgatCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCCTATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGcc                                                                                                                                                                                                          <  2:295750/251‑1 (MQ=255)                     aaaCAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTTTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGGCAGTGTCATCAAATTCAATATCATCAAAATTGATCTTCGGTTCGCCACGtt                                                                                                                                                                                                     <  1:95591/250‑1 (MQ=255)                      aaCAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGGTTCATCTTCGGATTCGTCAATTTCTTATTAGATTTCCGGTTTCTCATTGCGCTCTTTTTGCTCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGGTCGCCACATTc                                                                                                                                                                                                    >  2:207569/1‑251 (MQ=255)                               aGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCTGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCt                                                                                                                                                                                            <  1:23719/250‑1 (MQ=255)                                   ccAGAATGCCTACGCCCAGCATCATATTCTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGaa                                                                                                                                                                                       <  1:151921/251‑1 (MQ=255)                                       aaTGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTAATGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGGTCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGtt                                                                                                                                                                                   >  1:110409/1‑251 (MQ=255)                                           ccGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCaga                                                                                                                                                                               <  1:58707/251‑1 (MQ=255)                                                                 gACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGTCGCTCTACGCGCTCTTCATCCACGGTTTCATCTGCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTa                                                                                                                                                         <  2:309058/251‑1 (MQ=255)                                                                    gagaAGCGGGTTTACTGGCTGCTTTTTTGCGCATACGCGGACGACGCTCTACGCGATCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATGGTCAGTGTCATCAAATTCAATATCATCAAAATTGCTCTGCGGTTCTCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAAACCGACCacaatg                                                                                                                                                       >  1:98303/1‑245 (MQ=255)                                                                     agaAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGTCGAGAACGACCAGTacgac                                                                                                                                                     >  1:192797/1‑251 (MQ=255)                                                                      gaAGCGGGTTTTCTGGCTGCTTTTTTGCCCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTacgacg                                                                                                                                                    <  1:349564/251‑1 (MQ=255)                                                                         gcgggtgtatctgcttctTTTTTGCGCTTACGCCGAACACGCTCTTCGCGCTCTTCATCCTCGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTGGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATAGCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATc                                                                                                                                                 <  2:348820/235‑1 (MQ=255)                                                                                   tGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCg                                                                                                                                       >  1:361149/1‑251 (MQ=255)                                                                                        gctTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGttt                                                                                                                                  >  1:245247/1‑251 (MQ=255)                                                                                        gctTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATGCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGttt                                                                                                                                  <  1:336860/251‑1 (MQ=255)                                                                                                 cgcTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGtc                                                                                                                         >  1:174337/1‑251 (MQ=255)                                                                                                      ccgcgGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTg                                                                                                                    <  1:227644/250‑1 (MQ=255)                                                                                                             cgacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGa                                                                                                              >  1:29853/1‑250 (MQ=255)                                                                                                              gatccgtcaCGCGCTCTGCACCCACGGTTTCTTCTGCGGATTCGTCAATTTCGTCTCCGATTTCCGGGTACTCATTGCGCTCTTTTTGCGCACGAGGCGGACGTCGATCGTCAGTGGCATCAAATTCAATACCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCATTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAAt                                                                                                            <  2:169379/243‑1 (MQ=255)                                                                                                              gacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGAt                                                                                                                                                  <  2:142923/213‑1 (MQ=255)                                                                                                              gacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGAt                                                                                                                                                  >  1:142923/1‑213 (MQ=255)                                                                                                                 gCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGAAAGAAGATCGCTGGGAGCGGGTTTCAGCCCGACTTCTGGTTTGAAttac                                                                                                         >  2:123077/1‑249 (MQ=255)                                                                                                                             ctTCATCCACGGTTTCATCTTCTGATTCGGCAATTTCTTCGTCGATTTCCGGGTCCTCATTTCGCTCTGGTTGCGCACGAGTCGGACGGCGATCGTCAGGGTCATCAATGTCAATATCATCAAAATTGATCTGCGGGTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGATGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAac                                                                                               <  2:361646/249‑1 (MQ=255)                                                                                                                               tCATCCACGGTTTAATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCAACAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAacac                                                                                             <  1:371529/249‑1 (MQ=255)                                                                                                                                         tttCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGTCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAgg                                                                                 >  1:314612/1‑251 (MQ=255)                                                                                                                                                 tCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAAt                                                                            <  1:347184/248‑1 (MQ=255)                                                                                                                                                     aTTCGTCAATTTCGTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTTAATGCATCCATTGAGCTCCCCACGAAAAGGTTAATGCTTAcc                                                                     <  1:105747/251‑1 (MQ=255)                                                                                                                                                         gTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCCAATTCAATATCATCAACATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTACAAGGTTAATGCTTACCACGt                                                                 >  1:353369/1‑251 (MQ=255)                                                                                                                                                                   tcttcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGACTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAAc                                                       >  1:300169/1‑251 (MQ=255)                                                                                                                                                                     ttcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCt                                                     >  1:295789/1‑251 (MQ=255)                                                                                                                                                                      tcGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGGACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGAGCGCTGGGATCGGGTTTCAGCTCGTCTTCGGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCt                                                     <  1:371940/250‑1 (MQ=255)                                                                                                                                                                             ccGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTa                                              <  1:104318/250‑1 (MQ=255)                                                                                                                                                                                  tcatcactgacatctTTTTGCGCACGCGTCGGCCGGCGATCGTCAGTGTCATCAAATTCACTATCCTCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGCTCGCTGGGATCGGGTTTGCGCTCGTGTGCTGGTTTGAATTCATCCATTTAAACCCCCACTAACAGGTTAATGCTTACCACGTTGCATTTACCCTGCGCCTAATAGAc                                        >  2:170656/13‑251 (MQ=255)                                                                                                                                                                                  tCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGAc                                        >  1:25219/1‑251 (MQ=255)                                                                                                                                                                                          gcgcTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGAt                                <  1:119222/251‑1 (MQ=255)                                                                                                                                                                                          gcgcTATTTTTGCGCACGCGTCTGACCGCGATCGTCAGTTTCATCAAATTCAATATCAGCAAAATGTAGCTGCGGTTCGCCACGTTCAGTACGCTCATAACTTTCAGAAGATTGACGAGAACGACCAGTACGATGAGCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAAGGCTTACCACGTTGCCATTAACCTGAAGCTAATAGACCGCTTGa                                 <  2:224617/250‑1 (MQ=255)                                                                                                                                                                                           cgcTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGGTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTGCTGGTTTGAATTCATCCATTTAACACCCCAc                                                                                        <  1:227249/194‑1 (MQ=255)                                                                                                                                                                                           cgcTCTTTTTGCGCACGAGTCGGACGGAGATCGTCAGTGTCATCCAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGGTCAGTACGCTCAGAACGATCAGAAGATTGACGAGAACGACCAGTACGTCGATCGCTGGGATCGAGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCAc                                                                                        >  2:227249/1‑194 (MQ=255)                                                                                                                                                                                                               cGGACTGCGATCGTAAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGCTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGATTTAAGCTCGCATTCTGGTTTGAATTAATCCATTTAACACACCACTAAAAGGTTAATGCTTACCACGTTAAAATTAACCTGAAGATAATAAACCGCTTGATaaacgacgtgaccataattgt           >  2:284328/1‑232 (MQ=255)                                                                                                                                                                                                                   cGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCGCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTgttgtt       <  1:128511/251‑1 (MQ=255)                                                                                                                                                                                                                        aTCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTTCGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAATCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCt  <  1:305774/251‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          GCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCT  >  NC_000913/3516939‑3517404

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