Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,183	A→C	intergenic (+163/‑51)	thrC → / → yaaX	threonine synthase/DUF2502 domain‑containing protein YaaX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,183	0	A	C	100.0%	173.6 / NA	55	intergenic (+163/‑51)	thrC/yaaX	threonine synthase/DUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (31/24);  total (31/24)

GTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGC  >  NC_000913/4941‑5428                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTcc                                                                                                                                                                                                                                               >  2:525338/1‑251 (MQ=255)       aTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGatat                                                                                                                                                                                                                                         >  2:473191/1‑251 (MQ=255)           aCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAc                                                                                                                                                                                                                                     <  1:409543/251‑1 (MQ=255)             ccTTGTTTCCACAAAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCAGCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGTAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAccc                                                                                                                                                                                                                                   <  2:623051/251‑1 (MQ=255)                tgtttttaaaaagtgtggccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAAGCATCAGTAAAATCTATTCATTATCACAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAGCGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGGGAGTTGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCAt                                                                                                                                                                                                                                 <  2:236404/232‑1 (MQ=255)                          aaaTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATGAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATGACGGGCAATg                                                                                                                                                                                                                      <  2:595073/250‑1 (MQ=255)                            aTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                                                                                                                                                    >  1:537329/1‑251 (MQ=255)                                gCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                                                                                                                                                <  1:558458/251‑1 (MQ=255)                                  ccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAATCATTATCTCAATCAGGCCGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                                                                                                                                              <  1:229716/251‑1 (MQ=255)                                     ccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAACCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                                                                                                                                           <  1:473191/251‑1 (MQ=255)                                       gATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                                                                                                                                          >  1:581382/1‑250 (MQ=255)                                          tttGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑t                                                                                                                                                                                                     >  1:487551/1‑251 (MQ=255)                                            tgctgcGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑tat                                                                                                                                                                                                   >  1:329579/1‑249 (MQ=255)                                                     gcgtAAATTGAGGAGGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAATAACGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTGCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGAt                                                                                                                                                                                          <  1:36974/251‑1 (MQ=255)                                                       gtAAATTGATGATGAATAATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGc                                                                                                                                                                                        >  2:87016/1‑251 (MQ=255)                                                         aaaTTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCa                                                                                                                                                                                       >  1:346485/1‑250 (MQ=255)                                                              gatgatgaATCATCAGTAAAAACTATTCATTATCTCAATCAGGCCGGGTTTGCTTGTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCAACGCGTAATCTCCGGATAACGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATcga                                                                                                                                                                                  <  2:525317/250‑3 (MQ=255)                                                                 gatgaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCg                                                                                                                                                                               <  1:400548/250‑1 (MQ=255)                                                                    gaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTAc                                                                                                                                                                            <  1:164638/250‑1 (MQ=255)                                                                        catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCa                                                                                                                                                                       <  2:166611/251‑1 (MQ=255)                                                                           caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACtt                                                                                                                                                                    <  2:2479/251‑1 (MQ=255)                                                                                 actCTAGTCATTATCTCACTCAGGCCGGGTTTACTTATATGCAGCC‑GGCTTTTTTTTCAAGAAATTATGGAGAAAAACGACAGGGAAATAGGTGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGAGAATAACAACCGCCGTACTCATCGCGTAATCTCCGGATATAAACCCATAACGGGCAATGCTAACAGGA‑TCACCATATGAAAAAGCTGAAATCTATCGTACTCGAACTTTCCCTg                                                                                                                                                              >  2:109328/3‑251 (MQ=255)                                                                                  aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                                                                                                                                             >  2:178695/1‑251 (MQ=255)                                                                                        tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                                                                                                                                                       >  1:411235/1‑251 (MQ=255)                                                                                          aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTGTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                                                                                                                                                     <  2:212253/251‑1 (MQ=255)                                                                                                     cAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                                                                                                                                          >  2:358045/1‑251 (MQ=255)                                                                                                         ccGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGATAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAATATGCAATCTATCGTACTCGCACTTTCCCTCGTTCTGGCCGCTCCCATGgcagca                                                                                                                                      <  2:207439/251‑1 (MQ=255)                                                                                                                   ttttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg                                                                                                                            >  1:448274/1‑251 (MQ=255)                                                                                                                    tttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg                                                                                                                            >  1:329905/1‑250 (MQ=255)                                                                                                                                 gCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGt                                                                                                               <  2:305462/251‑1 (MQ=255)                                                                                                                                 gCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGACAGCACAGGCTGCGGAAATTACGTTAGt                                                                                                               >  1:329642/1‑251 (MQ=255)                                                                                                                                         aTGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTc                                                                                                         >  1:143343/1‑249 (MQ=255)                                                                                                                                                aaTTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAAtt                                                                                                <  1:178695/251‑1 (MQ=255)                                                                                                                                                     tGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTAc                                                                                              >  1:550629/1‑248 (MQ=255)                                                                                                                                                          aaaaaCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGc                                                                                      <  2:512233/251‑1 (MQ=255)                                                                                                                                                          aaaaaCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGc                                                                                      <  1:59013/251‑1 (MQ=255)                                                                                                                                                              acgaccgggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGCTATCGACCCATAACGGGCCATGATAAAAGGAGTAACC‑TATGCAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCAATGGCAGCACGGGCTGCGGAAATTAAGTTAGTCCCGTCAGTAAAATTAAAGATAGGCGATc                                                                                  >  1:312192/7‑251 (MQ=255)                                                                                                                                                                 aCAGGGAAAAAGGAGACATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTg                                                                               >  1:130876/1‑251 (MQ=255)                                                                                                                                                                       aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATc                                                                         >  2:179172/1‑251 (MQ=255)                                                                                                                                                                         aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGAGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATACCGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGt                                                                       >  1:89584/1‑251 (MQ=255)                                                                                                                                                                          aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTg                                                                      <  1:132884/251‑1 (MQ=255)                                                                                                                                                                                 aaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTAtt                                                                >  2:181712/1‑250 (MQ=255)                                                                                                                                                                                  aTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTa                                                               >  1:219553/1‑250 (MQ=255)                                                                                                                                                                                       cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TAGGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTg                                                            <  1:363592/248‑1 (MQ=255)                                                                                                                                                                                              tGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGt                                                  >  1:583388/1‑251 (MQ=255)                                                                                                                                                                                                      cTTAGAGCTTAGGACTGCGGAGAATAACAACAGTCGTTCTCATAGCGTAATCTCCGGATATCGACCAATAACGGGCACTGATAAAAAGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCAATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGACAAATTACAGCACGGCGAACGTGATAATCGTGGCTCTTACTGGGATGGAGGCCACTGgcg                                          >  2:275291/1‑251 (MQ=255)                                                                                                                                                                                                             aTTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATACTCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCAc                                   >  2:499542/1‑251 (MQ=255)                                                                                                                                                                                                                 ggATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCTGATATCGACCCATAACGGGCAATGATAAAATGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCTCACTTTCCCTTGTTCTTGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCTTCAGTAAAATTACAGATAGGCGATCGTg                                                                               <  2:60704/203‑1 (MQ=255)                                                                                                                                                                                                                 ggATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGAACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTg                                                                               >  1:60704/1‑203 (MQ=255)                                                                                                                                                                                                                      gCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACGGGGATGGAGGTCACTGGCGCGACCACGGCtggtg                           <  1:165200/250‑1 (MQ=255)                                                                                                                                                                                                                       cGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACACATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGACTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGAGCGTGATAATCGTGGCTATTATTGGGTTGGAGGTCACTGGCGCGACCACGACtggtgg                          >  1:478228/1‑250 (MQ=255)                                                                                                                                                                                                                        ggAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGaa                        <  1:401410/251‑1 (MQ=255)                                                                                                                                                                                                                                       tcgTTCTCATCGCGAAATCTCCGGATATAGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTAGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATAGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATg         >  2:286175/2‑251 (MQ=255)                                                                                                                                                                                                                                         gTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGc       >  1:227347/1‑251 (MQ=255)                                                                                                                                                                                                                                               aTCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATACAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGc  >  2:12607/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGC  >  NC_000913/4941‑5428

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