Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	674,518	G→A	A89V (GCG→GTG)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	674,518	0	G	A	95.7%	247.9 / NA	71	A89V (GCG→GTG)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/30);  minor base C (1/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATC  >  NC_000913/674282‑674757                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                tcgtcagtccggtgtggcggtgggggCACCCTTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACACCTCGGTGAATAATTTCTTTTCCCAACGGTATTTTTCCGGCGTACAGGTTGCCGGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATCCGCGATGTTGTCGTCCGTCCATGGTGCCG‑GGGCGGTGTTGTTTTTCTCCGCCACGCCTTCCgcagg                                                                                                                                                                                                                                    <  2:379292/227‑1 (MQ=255)   gCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATCCCAGTCCTTTTTTATACTGCCCGGTGAAGATTTTCTGTCCCCAGCGGAAGTATTCCGGCGTACGGGTTGCCAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCTTTTTGAGCTGGTTTTTCATAGACGCGATGTTGCCGTGCGTCCACGGGGCGG‑GGGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggcag                                                                                                                                                                                                                                 <  2:152470/251‑1 (MQ=255)   gCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGGGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCgcaggc                                                                                                                                                                                                                                   >  2:434056/1‑249 (MQ=255)             tCTGGACGTTCGGGCCCAAGTTGCCCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTCCCCACCGGTAGTATTCCGGCTTACAGGTGGCCAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGTTTTTTCATATACGAGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAAcgc                                                                                                                                                                                                                        <  1:544882/250‑1 (MQ=255)                       cGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCaatacc                                                                                                                                                                                                             >  1:329687/1‑246 (MQ=255)                       cGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCGCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCg                                                                                                                                                                                                              <  1:131428/250‑1 (MQ=255)                                    cccGCAGAAGTCTTCTTATATCCCAGGCCTTTTTTATACAGCTCGGTAAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGAc                                                                                                                                                                                                <  1:280968/250‑1 (MQ=255)                                      cGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGt                                                                                                                                                                                              >  1:515808/1‑251 (MQ=255)                                          gaattcttctTATATCCCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGCCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGttttt                                                                                                                                                                                          <  1:129373/247‑1 (MQ=255)                                                 tcttATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAg                                                                                                                                                                                    <  1:323079/250‑1 (MQ=255)                                                     atataCCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATAc                                                                                                                                                                               <  1:401803/251‑1 (MQ=255)                                                       ataCCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGc                                                                                                                                                                             >  1:5291/1‑251 (MQ=255)                                                         aCCAGGCTTTTTTTATACAGCGCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCTTTTTGAGCTGGTTTTTCATATACGCGATTTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCt                                                                                                                                                                            <  1:508097/250‑1 (MQ=255)                                                           cAGGCCTTTTTTATACAGCTCGGTAAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCAAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGTGGTTGTCGTACGTCCACGGTGCCG‑GGGCGGTGTTGTTTTCCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGt                                                                                                                                                                         <  1:75868/251‑1 (MQ=255)                                                             ggCCATTTTTCTACAGCTAGGTGAAGAATTTCTGCTCACAACGGTAGTATTCCGGCGTACAGGGTGACAGCTCGCGGATCCAGGCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGGCGTACGGCCACGGTGCCG‑GAACGGAGTTGTTTTTCACCTCCATGCATTACGCAGGTACACCAACCGTGTACCATCCGATAGGCTGCAGGCCGTATTTGCCAAGAATAcacaggt                                                                                                                                                                         >  2:106683/1‑243 (MQ=255)                                                                cTTTTTTATACAGCTAGGAGAAGAATTTCTGCTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATCACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCCATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGAGTTGTTTTTACCCGCCCCGCCTTCCGCAGGCAGACCCGACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCAAGCATGCGCCGGTAgagg                                                                                                                                                                    >  2:454584/1‑248 (MQ=255)                                                                    tttATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCg                                                                                                                                                                 >  1:216557/1‑250 (MQ=255)                                                                        aaCAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCATTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCAc                                                                                                                                                            <  1:444016/250‑1 (MQ=255)                                                                         aCAGCTCGGTGAAGATTTTCTGTTCCCACCGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATCCGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCAc                                                                                                                                                            <  1:434056/250‑1 (MQ=255)                                                                          cagttCGGTGAAGAATTTCTCTTCCCACCGTTGTTCTTCCGCCTTCCCGTTTCCCACCTCGCGGTTCCAGTCATAGCCAACCCCCAGCCTTTTGACCTGGTTTTTCTTTGCCCCGATGTTGTCGTTCATCACCGTTGCCG‑GCGCGTTTTTGTTTTTCCCCGCCCCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGt                                                                                                                                                          <  2:390097/247‑1 (MQ=255)                                                                            gCTCGGTGAAGAATTTATGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGACAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGGTTTTCATATACGCGATGGTGTCGTACGTCCACGGTGCCG‑GAGCGGATTTGTTTTTCACCGCCACACCTTCCGCAGGCAGACCAAACGCGTCCCCGCCGATCGACTGCAGGACGTTTTTGACCAGCCTACGCTGGCAGCGGGCGATCACGTc                                                                                                                                                         >  2:68586/1‑250 (MQ=255)                                                                             cTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCAc                                                                                                                                                       >  2:113741/1‑251 (MQ=255)                                                                              ccGGTGAATAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCGCGCGGCTCCAGTCATAACCAAAGCCCAGCTTTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCAcc                                                                                                                                                      <  1:7055/250‑1 (MQ=255)                                                                               cGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCg                                                                                                                                                     >  1:225883/1‑251 (MQ=255)                                                                                     agaaTTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGGTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATgg                                                                                                                                                 >  1:305691/1‑249 (MQ=255)                                                                                              gTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTta                                                                                                                                         <  1:268504/248‑1 (MQ=255)                                                                                                       ggTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGccc                                                                                                                             >  1:313010/1‑251 (MQ=255)                                                                                                              ttCCGGCGTCCAGGTTGCCAGCTCGCGGCTCCAGTCATACCCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATgtgt                                                                                                                       <  1:61067/250‑1 (MQ=255)                                                                                                                ccGGCGTACAGGTGGCCAGCTCGGGGCTCAAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGTTGTTGTCGTACGTCCACGGGGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGt                                                                                                                    <  1:166480/251‑1 (MQ=255)                                                                                                                   gCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCg                                                                                                                  >  1:119708/1‑250 (MQ=255)                                                                                                                     gTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTCGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGCGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGCTGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGAcc                                                                                                               >  1:134406/1‑251 (MQ=255)                                                                                                                       aCAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTACAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGAcc                                                                                                               >  1:331674/1‑249 (MQ=255)                                                                                                                             tccctgctcgggtccccAGACATAACCAAACCCCTGTATTTTGGTCGGGTTTTTCATATGCGCGATGTTGGCGGACGTCCAGGGGGCCG‑GAGGGGGGTTGTTTTTACCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGa                                                                                                        <  2:25036/235‑1 (MQ=255)                                                                                                                               ccAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTT‑GAGCTGTTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTCCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAgg                                                                                                    <  1:107755/251‑1 (MQ=255)                                                                                                                                  gcgcgcgGCTCCAGTAAACACAAAAGCCCAGATTTTTGAGGTGGTTTTTCATATACGCGATGTGGTCGTACGTCCAGGGGGCCG‑GAGCGGTGTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGa                                                                                                  <  2:7185/248‑1 (MQ=255)                                                                                                                                       cgGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGGTGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCAt                                                                                             >  1:497231/1‑251 (MQ=255)                                                                                                                                              agtgaTAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATCCGCTATGTTTGCGTACGTCC‑CGGGGCTGTGAGGGGTGTTGTTTTTCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAgg                                                                                      <  1:415945/247‑1 (MQ=255)                                                                                                                                               gTCATACCCAAAGCCCAGCATTTGGAGCTGGTTTTTCATATACGCGTTGTGGTCGTACGTCCACGGTGCGG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGc                                                                                     <  1:513247/251‑1 (MQ=255)                                                                                                                                                        aaaGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGAAGTTACGTACGTGGCCCATGAGTAGTCGACCCGAAGGCTCGGGAAGAATAGACAGGCCGTAAtcat                                                                            >  2:468752/1‑248 (MQ=255)                                                                                                                                                        aaaGCCCAGCATTTTGAGCTGGTTTTTAATATACGCGATGTTGTCGTACGTCCACGGAGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCGGCAGGACGTTTTTGCCCAGCATACGTTGGGAGCGGGCGATCACGTCACCGCTGGTGTAGTTACGTACGTAGCCCATGTGCAGTCGACCACAAGGATAGGGGAGCATAGACAGGCAGTaacgct                                                                            >  2:248961/1‑247 (MQ=255)                                                                                                                                                             ccAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTctct                                                                       >  2:340945/1‑251 (MQ=255)                                                                                                                                                              cAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GGGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCtt                                                                      <  1:241090/251‑1 (MQ=255)                                                                                                                                                                gCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCCCGGTGCCG‑GGGCGGTTTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTg                                                                    <  1:255754/251‑1 (MQ=255)                                                                                                                                                                       gAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGt                                                             >  2:88285/1‑251 (MQ=255)                                                                                                                                                                           ggttttttCATATGCGCGATGTGGTCTTACGCCCACGGTGCCG‑GAGGGGTGTTGTTTTTCACCGCCTCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTc                                                         <  1:483552/248‑1 (MQ=255)                                                                                                                                                                            ggTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTCTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGGAATACTTATCTTTGCTCTCGTCTTCg                                                        >  2:264867/1‑251 (MQ=255)                                                                                                                                                                               ttttCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGtt                                                     >  2:538080/1‑251 (MQ=255)                                                                                                                                                                                tttCATATACGCGATGTTGACGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATCCTTCTCTTTGCTCTCGTCTTCgggt                                                     >  2:95111/1‑248 (MQ=255)                                                                                                                                                                                 ttCATATGCGCGATGTTGTCGTACGTCCACGGTGCGG‑GGGCGGGGTGGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTa                                                    <  1:69026/250‑1 (MQ=255)                                                                                                                                                                                 ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTAc                                                   >  1:456095/1‑251 (MQ=255)                                                                                                                                                                                 ttCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTAc                                                   >  1:139902/1‑251 (MQ=255)                                                                                                                                                                                  ccatatcggcgaggttttcgtccatCCACGTTGCGG‑GAGGGGTGTTTTTTTCCACCGCCTCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAATTATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGGCTTCGGTTAc                                                   <  1:76435/226‑1 (MQ=255)                                                                                                                                                                                    atatACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATAATTATCTTTGCTCTCGTCTTCGGTTACttt                                                >  2:6232/1‑250 (MQ=255)                                                                                                                                                                                    atatACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTAACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGAACGTGGCCCATGTGTAGTCGAACAGAAGGATAGGGAAGAATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTc                                                >  2:402376/1‑251 (MQ=255)                                                                                                                                                                                       taCGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGCCAGGCAGTAATACTTCTCTTTGCTCTCGGCTTCGGTTACTTCaa                                              >  2:35478/1‑250 (MQ=255)                                                                                                                                                                                       taCGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCaa                                              >  1:208284/1‑250 (MQ=255)                                                                                                                                                                                        ccgcgATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAt                                            <  1:27952/250‑1 (MQ=255)                                                                                                                                                                                        aCGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTTTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCaaa                                             <  1:454337/250‑1 (MQ=255)                                                                                                                                                                                         cgcgATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGCGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAAtg                                           >  1:73730/1‑251 (MQ=255)                                                                                                                                                                                            gATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGCCAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTgc                                        >  2:414992/1‑251 (MQ=255)                                                                                                                                                                                             aTGTTGTCGTACGTCCACGGTGCCG‑GTGCGGTGTTGTTTTTCCCCGCCACGCCTTACGCAGGCAGACCAAAAGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTAAGTGGCCCATGTGTAGTCGACCAGAAGGGTAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGcgcg                                       >  2:404661/1‑251 (MQ=255)                                                                                                                                                                                                 tgtCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTtc                                   >  1:43791/1‑251 (MQ=255)                                                                                                                                                                                                        cGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATccc                            >  2:435966/1‑251 (MQ=255)                                                                                                                                                                                                          ccctcGGTGCCG‑GAGCGGTTTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCaa                          <  1:234680/247‑1 (MQ=255)                                                                                                                                                                                                           ccgcGGTGCCG‑GGGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAAt                         <  1:206149/248‑1 (MQ=255)                                                                                                                                                                                                           ccACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTTGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAAt                         >  1:194601/1‑251 (MQ=255)                                                                                                                                                                                                               ggTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGa                       >  2:73949/1‑249 (MQ=255)                                                                                                                                                                                                               ggTGCCG‑GAGCGGTGTTGTTTTTCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGa                       >  2:4333/1‑249 (MQ=255)                                                                                                                                                                                                                                  tttttCACCGCCCCGCCCTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCAGAGACAGGGAGTAAGACTTCTCTTTGCTCTCGTCTGCGGTTAGTTCAAATGTGCGCTTCTGATCCCAATGAAGCTGTTCTTTGGATTCTa    <  2:142762/250‑1 (MQ=255)                                                                                                                                                                                                                                   ttttCCCCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTAt   <  1:113741/250‑1 (MQ=255)                                                                                                                                                                                                                                    tttCACCGCCACGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTCCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTAtc  <  1:340945/250‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                TCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCG‑GAGCGGTGTTGTTTTTCACCGCCGCGCCTTCCGCAGGCAGACCAAACGCGTCCCAGCCGATCGGCTGCAGGACGTTTTTGCCCAGCATACGCTGGTAGCGGGCGATCACGTCACCGATGGTGTAGTTACGTACGTGGCCCATGTGTAGTCGACCAGAAGGATAGGGAAGCATAGACAGGCAGTAATACTTCTCTTTGCTCTCGTCTTCGGTTACTTCAAATGTGCGCTTCTCATCCCAATGAAGCTGTACTTTGGATTCTATC  >  NC_000913/674282‑674757

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: ‑