Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,600,215	+C	intergenic (‑26/+23)	yjjA ← / ← dnaC	DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,600,214	1	.	C	100.0%	221.8 / NA	65	intergenic (‑25/+24)	yjjA/dnaC	DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (29/36);  total (29/36)

GCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑CAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTG  >  NC_000913/4599966‑4600456                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   ggggcgtgcgtgtggtcaccggtatcagtgtcggttttgtcttctacttactggaccagatcttcggcccgctgacgttgttttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcgtcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaaggaacgaataaagggagcgaaagcgctctttattattgttagcaaagtgtgcgtcgttcattcctgaTGTGTCg‑cc                                                                                                                                                                                                                                                  <  1:239314‑M2/9‑2 (MQ=255)  gggcgtgcgtgtggtcaccggtaccagtttcggttttgtcttctacgtcctggaccagatcttcggcccgctgacgttggtttagggcaccccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctggtaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctcccttgattattgttagcaaagtgtgcttcgttcattcctgaTGTGTcg‑cca                                                                                                                                                                                                                                                 <  2:243669‑M2/10‑3 (MQ=255)      gtgcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGtat                                                                                                                                                                                                                                              <  1:102673‑M2/14‑1 (MQ=255)            gtggtcaccggtatcagtttcggtttgttcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGcc                                                                                                                                                                                                                                        <  1:424567‑M2/20‑1 (MQ=255)              ggtcaccgttaccagtttcggttttgccttctacgtacgggaccagatcgtcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgtgcattcctgaTGTGTCGCCAGTATACCGCCa                                                                                                                                                                                                                                       <  2:87139‑M2/21‑1 (MQ=255)              ggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgacagcttcttcttaatcagcctgtggctgttcatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCcat                                                                                                                                                                                                                                      >  1:546244‑M2/230‑250 (MQ=255)                 caccggtatcagtttcggttttgtcttctacgtactggaacagatcttcggcccgctggcgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgcgaaaatcgcaaccgaatgaacgaataacgggagcgaaagctccctttattattgttagcaaagtgagcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGaa                                                                                                                                                                                                                                    >  2:518521‑M2/227‑250 (MQ=255)                   ccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACg                                                                                                                                                                                                                                 >  1:454199‑M2/225‑251 (MQ=255)                               tcggttttgtcttctccgtactggaccagatctgcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcatgcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATAc                                                                                                                                                                                                                      <  1:60686‑M2/38‑1 (MQ=255)                                cggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatctgcgcactgttgccaagcgccagcttcttcttactcagcctgtgactgttaatgagaaaatcgtaaccgaatgaacgactaaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATActc                                                                                                                                                                                                                    >  2:396153‑M2/212‑251 (MQ=255)                                         cttctacgtactggaccagatcttcggcccgctgacgttggtttatggcctcccgccgatcatcggcgcactgttgccaagcgcccgcttcttcttaatcagcctgtggctgttaatgagaaaatagtaaccgaatgaacgactaaagggcgcgaaagctacctttattctagttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGCTTAATACTCTTTACCTGt                                                                                                                                                                                                           >  1:570257‑M2/203‑251 (MQ=255)                                           tctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                            <  1:593309‑M2/48‑1 (MQ=255)                                           tctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGtt                                                                                                                                                                                                          <  1:93356‑M2/50‑1 (MQ=255)                                               cgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCg                                                                                                                                                                                                     >  1:139920‑M2/197‑251 (MQ=255)                                                gtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCgg                                                                                                                                                                                                    >  2:443441‑M2/196‑251 (MQ=255)                                                                 cgccccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCa                                                                                                                                                                                   <  1:54752‑M2/73‑1 (MQ=255)                                                                  ggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                                                                    >  1:475964‑M2/178‑249 (MQ=255)                                                                    cccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGCTTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGtt                                                                                                                                                                                >  1:114443‑M2/176‑251 (MQ=255)                                                                      cgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgcTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGATAGCTATCCCAGTTg                                                                                                                                                                               >  2:445883‑M2/174‑250 (MQ=255)                                                                          gacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAt                                                                                                                                                                          <  1:315569‑M2/82‑1 (MQ=255)                                                                                 gtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcac                                                                                                                                                                   <  1:417641‑M2/89‑1 (MQ=255)                                                                                          atcccgccgctcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGtt                                                                                                                                                           >  1:491813‑M2/154‑250 (MQ=255)                                                                                                         ggcgcactgttgccaatcgccagcttcttcttactcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACgg                                                                                                                                           >  2:176402‑M2/139‑251 (MQ=255)                                                                                                            gcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTcc                                                                                                                                        <  2:427477‑M2/116‑1 (MQ=255)                                                                                                                      caagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGtt                                                                                                                              <  1:425485‑M2/126‑1 (MQ=255)                                                                                                                      caagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAGCTCCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGtt                                                                                                                              <  1:556958‑M2/126‑1 (MQ=255)                                                                                                                            ccagcttcttcttcatcagcctgtggctgttaatgggaacctcgtaaccgaatgaacgaataaagggagcgcaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcgc                                                                                                                                 >  2:604761‑M2/120‑242 (MQ=255)                                                                                                                            ccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcgc                                                                                                                                 <  1:604761‑M2/123‑1 (MQ=255)                                                                                                                            ccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCa                                                                                                                        <  1:32172‑M2/132‑1 (MQ=255)                                                                                                                              agcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGc                                                                                                                      <  1:340390‑M2/134‑1 (MQ=255)                                                                                                                                 ttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAAc                                                                                                                   >  1:277838‑M2/115‑251 (MQ=255)                                                                                                                                        taatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTGCCTGTTACCCGGCTACGGTAGCTATCCGAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCGGCAACTTGGTc                                                                                                             <  1:98942‑M2/143‑1 (MQ=255)                                                                                                                                        taatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCTCGCGTTCGCCCAGCAACTTGGt                                                                                                              <  1:176402‑M2/142‑1 (MQ=255)                                                                                                                                         aatcagactgtggctgttaatgagaaaatcgtagccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCCACTTGGTCAt                                                                                                           >  2:513858‑M2/107‑251 (MQ=255)                                                                                                                                            cagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtct                                                                                                         <  1:468730‑M2/147‑1 (MQ=255)                                                                                                                                                             tgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTgg                                                                                       >  1:436675‑M2/87‑251 (MQ=255)                                                                                                                                                              gagaaaagcgtaaccgaatgaacgaataaagggagcgaaatctccctttattattggtagcaaagtgtgcttcgttcattcctgaTGTGTGGCCAGTATACCGCCAGAAACGGGGATTAATTCTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTTAAGATCACCCGCAAACTGTTTCCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAGCTCTTCCATATTACTGTTGGt                                                                                      <  2:182392‑M2/166‑1 (MQ=255)                                                                                                                                                               agaaaatcgtaaccgaatgaccgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcatgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTc                                                                                     >  2:276362‑M2/85‑251 (MQ=255)                                                                                                                                                                  aaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACTGTCCATCACGCGTTCGCCCAGCAACTTGTTCATCTCTTCCATATTACTGTTGGTCAGc                                                                                  <  2:575443‑M2/170‑1 (MQ=255)                                                                                                                                                                          accgaatgaacgaataatgggagcgaaagctccctttattattgttagccaagtgtgcttcgttcattcttgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCGAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGt                                                                          >  2:98047‑M2/74‑251 (MQ=255)                                                                                                                                                                               atgaaagaataaagggagcgaaagctccctttcttattgttagccaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAACCGGGGATTCATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGCCGCATACGGTCCCTCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGc                                                                     >  2:48612‑M2/69‑251 (MQ=255)                                                                                                                                                                                 gaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGt                                                                   <  2:454199‑M2/185‑1 (MQ=255)                                                                                                                                                                                         aaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                                           <  2:277838‑M2/193‑1 (MQ=255)                                                                                                                                                                                         aaagggagcgaaaaatccctttattattgttagccaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCCGACACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGga                                                              >  2:238998‑M2/59‑248 (MQ=255)                                                                                                                                                                                             ggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAagag                                                       <  1:26721‑M2/197‑1 (MQ=255)                                                                                                                                                                                               agcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcg                                                     <  1:27120‑M2/199‑1 (MQ=255)                                                                                                                                                                                                  gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                                  >  2:492127‑M2/50‑251 (MQ=255)                                                                                                                                                                                                     agctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATc                                               >  2:318308‑M2/47‑251 (MQ=255)                                                                                                                                                                                                           ctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGAt                                         >  2:365201‑M2/41‑251 (MQ=255)                                                                                                                                                                                                             ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACGGTTACCCAGGCGCATACGGTCCAGCACGCGTGCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCAGCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGAt                                         <  1:596611‑M2/211‑1 (MQ=255)                                                                                                                                                                                                                  attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGATCGGCGATCGACGATCTGGTtg                                  <  1:505387‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                  attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGtt                                   <  1:491158‑M2/217‑1 (MQ=255)                                                                                                                                                                                                                    tgttagccaaatgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAACTGTTACCCGGCTCCGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgat                                >  2:534051‑M2/32‑251 (MQ=255)                                                                                                                                                                                                                       tagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTgatgat                             >  1:140224‑M2/29‑251 (MQ=255)                                                                                                                                                                                                                           aaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACt                         <  2:361505‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                     tcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAt                <  2:475964‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                     tcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAt                <  1:383736‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                                      cgttccttcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTACTACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAttt              >  2:545480‑M2/14‑251 (MQ=255)                                                                                                                                                                                                                                      cgttcagtgctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTGTGGAAGAAGATCGGCGATCGACGATCTGGTTTATGATCACTTTTTCGTAttt              <  2:588646‑M2/238‑1 (MQ=255)                                                                                                                                                                                                                                       gttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGGTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTc             <  2:229222‑M2/239‑1 (MQ=255)                                                                                                                                                                                                                                           attccggaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTGCCCGGCTACGGTAGCTATCCCAGTTGAAGATCTCCCGCAAACTGTTTCCCAGGCGCATGCGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGGAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGAtt         <  1:303632‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                             tcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTc        <  1:396153‑M2/244‑1 (MQ=255)                                                                                                                                                                                                                                                tgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAttt              <  2:462391‑M2/238‑1 (MQ=255)                                                                                                                                                                                                                                                tgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAttt              >  1:462391‑M2/4‑241 (MQ=255)                                                                                                                                                                                                                                                tgaTGTGTCGCCAGTATACAGCAAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTc    >  1:172397‑M2/4‑251 (MQ=255)                                                                                                                                                                                                                                                 gaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCt   <  1:342594‑M2/249‑1 (MQ=255)                                                                                                                                                                                                                                                  aTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGACCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTg  >  1:589487‑M2/2‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   GCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑CAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTG  >  NC_000913/4599966‑4600456

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