Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,140,124	C→A	D700Y (GAC→TAC)	ptsA ←	putative PTS multiphosphoryl transfer protein PtsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,140,124	0	C	A	100.0%	244.0 / NA	69	D700Y (GAC→TAC)	ptsA	putative PTS multiphosphoryl transfer protein PtsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/42);  total (27/42)

TAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCAC  >  NC_000913/4139989‑4140258                                                                                                                                        |                                                                                                                                       tAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCtt                                                                                                                                    <  2:150504/140‑1 (MQ=255)  aaaGCCCACAAAGGTAGAGAAACCGACCTCGCGCCCCCCCAAGTCGGCTTCCATTTTACGGGGATAGCGGCACCGCCCCGCCAGCAGCAGGTTGTCGGTCATGCCTTTGAGCTCTTCTTCTTTGCTGCGCCAGTAGCttt                                                                                                                                   <  2:1579953/140‑1 (MQ=255)    aGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTcc                                                                                                                                 >  2:1478789/1‑140 (MQ=255)    aGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTcc                                                                                                                                 >  1:1112736/1‑140 (MQ=255)           aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGa                                                                                                                          >  1:943385/1‑140 (MQ=255)           aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGa                                                                                                                          >  1:1579953/1‑140 (MQ=255)           aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGa                                                                                                                          >  2:1185403/1‑140 (MQ=255)           aCCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGa                                                                                                                          >  2:972356/1‑140 (MQ=255)               gTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCa                                                                                                                      >  1:1332157/1‑140 (MQ=255)                          ggCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGa                                                                                                           >  1:161656/1‑140 (MQ=255)                          ggCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCACCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGa                                                                                                           >  2:718364/1‑140 (MQ=255)                            ccTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTcggcg                                                                                                               <  1:1548004/134‑1 (MQ=255)                            ccTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTcggcg                                                                                                               >  2:1548004/1‑134 (MQ=255)                             cTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCa                                                                                                        <  2:1112736/140‑1 (MQ=255)                               cgcgcgCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGc                                                                                                      <  1:782754/140‑1 (MQ=255)                               cgcgcgCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGc                                                                                                      >  1:1267303/1‑140 (MQ=255)                                gcgcgcCGACAAGTCGGTTTCAAGTTTATGCGGATGGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCg                                                                                                     <  2:144360/140‑1 (MQ=255)                                  gcgcCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGt                                                                                                   <  1:608255/140‑1 (MQ=255)                                  gcgcCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGt                                                                                                   <  1:1033017/140‑1 (MQ=255)                                  gcgcCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGt                                                                                                   <  1:483988/140‑1 (MQ=255)                                     cccACAAGTCGGTTTCCAGTTTACGCGGATAGCGGCAGGGCCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCATCGGTGCg                                                                                                <  2:400250/140‑1 (MQ=255)                                     cccACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCGCTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCg                                                                                                <  2:381262/140‑1 (MQ=255)                                     cccACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCg                                                                                                <  2:1332157/140‑1 (MQ=255)                                     cccACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGAGACCAGCGGTGCg                                                                                                <  2:131052/140‑1 (MQ=255)                                      ccACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGt                                                                                               <  2:1512603/140‑1 (MQ=255)                                      ccACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGt                                                                                               <  1:872025/140‑1 (MQ=255)                                         caAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCtt                                                                                            >  2:677372/1‑140 (MQ=255)                                            gTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCttgtt                                                                                         <  1:1048861/140‑1 (MQ=255)                                              cGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGt                                                                                               <  2:445144/132‑1 (MQ=255)                                              cGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGt                                                                                               >  1:445144/1‑132 (MQ=255)                                               ggCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTggg                                                                                      <  2:943385/140‑1 (MQ=255)                                               ggCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTggg                                                                                      <  2:585680/140‑1 (MQ=255)                                                    ccAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATg                                                                                 >  2:234208/1‑140 (MQ=255)                                                             gcgGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTgc                                                                        <  1:972356/140‑1 (MQ=255)                                                             gcgGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTgc                                                                        <  1:1185403/140‑1 (MQ=255)                                                                 aTAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCca                                                                    >  1:384382/1‑140 (MQ=255)                                                                    ggggctgcgggCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGTTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagca                                                                 <  1:803354/138‑1 (MQ=255)                                                                     cggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagcag                                                                <  1:677372/140‑1 (MQ=255)                                                                     cggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagcag                                                                <  2:39249/140‑1 (MQ=255)                                                                     cggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagcag                                                                <  2:1267303/140‑1 (MQ=255)                                                                     cggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagcag                                                                <  1:1614162/140‑1 (MQ=255)                                                                     cggcagcggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCcagcag                                                                <  1:177356/140‑1 (MQ=255)                                                                           cggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTc                                                          >  1:1147548/1‑140 (MQ=255)                                                                            ggcCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCt                                                         <  1:458858/140‑1 (MQ=255)                                                                               ccGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACt                                                      >  1:408222/1‑140 (MQ=255)                                                                                 gCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTc                                                    <  2:384382/140‑1 (MQ=255)                                                                                      cagcagGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCg                                               <  2:408222/140‑1 (MQ=255)                                                                                               aTCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGc                                      <  1:818675/140‑1 (MQ=255)                                                                                               aTCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGc                                      <  1:618151/140‑1 (MQ=255)                                                                                                   gTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGc                                  <  2:161656/140‑1 (MQ=255)                                                                                                      aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCa                                >  1:287941/1‑139 (MQ=255)                                                                                                      aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCa                                <  2:287941/139‑1 (MQ=255)                                                                                                      aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCAt                               >  2:1165091/1‑140 (MQ=255)                                                                                                      aTGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCAt                               >  2:579506/1‑140 (MQ=255)                                                                                                       tGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCAtt                              >  2:1541771/1‑140 (MQ=255)                                                                                                        gCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTg                             <  2:1517732/140‑1 (MQ=255)                                                                                                          cTTTGAGCACTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGcc                           <  2:374218/140‑1 (MQ=255)                                                                                                                  aCTTCTTCTTTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAg                   <  2:1147548/140‑1 (MQ=255)                                                                                                                    ttcttcttTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCa                 <  1:1478789/140‑1 (MQ=255)                                                                                                                         cttTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAg                   >  1:1574361/1‑133 (MQ=255)                                                                                                                         cttTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAg                   <  2:1585331/133‑1 (MQ=255)                                                                                                                         cttTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAg                   >  1:1585331/1‑133 (MQ=255)                                                                                                                         cttTGCTGCGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAg                   <  2:1574361/133‑1 (MQ=255)                                                                                                                              ctgcGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCa       >  2:119419/1‑140 (MQ=255)                                                                                                                              ctgcGCCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAAATTGCGGCa       >  2:1481252/1‑140 (MQ=255)                                                                                                                                 cgcCAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTc    >  2:1348177/1‑140 (MQ=255)                                                                                                                                   ccAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCAc  <  1:1329190/140‑1 (MQ=255)                                                                                                                                   ccAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCAc  <  2:114278/140‑1 (MQ=255)                                                                                                                                   ccAGTAGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTACCTGGTTGAGCAACTTGCGGCACTCAc  >  2:150482/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       TAAAGCCCAGACCGGTAGAGAAAACGGCCTCGCGCGCCCACAAGTCGGCTTCCAGTTTACGCGGATAGCGGCAGCGGCCCGCCAGCAGCAGGTTATCGGTCATGCCTTTGAGCACTTCTTCTTTGCTGCGCCAGTCGCTTTCCAGTGTGATGCACTCGGCGGTGACCAGCGGTGCGTCTTGTTGGGTCATGCGGAATTGCGCCAGCAGGTGTTCTACTTCCAGCGAAGTACGGCAGGCCATTGCCTGGTTGAGCAACTTGCGGCACTCAC  >  NC_000913/4139989‑4140258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑