Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	156.4 / NA	42	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/27);  total (15/27)

GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATG  >  NC_000913/1652092‑1652574                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      gtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTAt                                                                                                                                                                                                                                            <  1:349663/251‑1 (MQ=255)                                      gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGATAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGttttt                                                                                                                                                                                                      <  1:517768/250‑1 (MQ=255)                                        gtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCaa                                                                                                                                                                                                   <  2:331488/251‑1 (MQ=255)                                                          aCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                                 <  1:87431/251‑1 (MQ=255)                                                            aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGt                                                                                                                                                                               <  2:41116/251‑1 (MQ=255)                                                            aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTTAAGTGTCAGCGAGAGt                                                                                                                                                                               <  2:297075/251‑1 (MQ=255)                                                              aGATCTGAAAAATGGAACAATGATGATCAGGAGAAAATTAGCAAAAGACCGTTTCACGGTACCAAAAACACATGCGGGAACCAATCGGGTCATTCATCTTATTAAGCAAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATTGATGTTCGCCTCATAGAGTATGG‑CAGAACAGAAAAGCAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                                                                                                                                                                              <  2:388498/250‑1 (MQ=255)                                                                    gAAAAATGGAACAATGATGATCAGGAGAAATNTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑GGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAtt                                                                                                                                                                       <  2:266807/251‑1 (MQ=255)                                                                      aaaaTGGAACAATGATGATCAGGAGAAATTTCGCAAAAGACCGTTTAACGGAACCAAAAACACAGGCGGAAACCAATAGAAGCATTCCTCTTATTCAGCCAGCAATCACCGCACTCAGGGGTCAGGTGACATTAACGAGGCAGAACAAAAAGACTATAATTAATGGTAACCTCACAAAGTGTGGGCGGAA‑AGAAAAACAAAAA‑TAGACAGTTGTTTTTCAACCTGAAGTTACATCGAAAGTaaaaaacaat                                                                                                                                                                     >  2:97143/1‑247 (MQ=255)                                                                          tGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGACCCAATCGGGGCATTCATCTTATTAAGCCAGCAAACGACGCTCTCCGGAGTTAGATGACATTAAAGAGCCTGAGCAAAAAGCATATCATTTATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAAACAAAA‑TGCACCTTTGTTTTTCAAACTGAAGTGGAAGCGAGAGTAAAAAAtaaagggg                                                                                                                                                                 >  1:77477/1‑244 (MQ=255)                                                                                    atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCt                                                                                                                                                                                                <  1:22998/210‑1 (MQ=255)                                                                                      gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACTCAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCGGCAATCGACGCTCTCTGGAGTCATATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCTCCTCAGAGAGTATGG‑CAGGACAGAAGGACAGAAA‑TGCACCTTTGTTGTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc                                                                                                                                                      <  1:495449/250‑1 (MQ=255)                                                                                             agaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTc                                                                                                                                              <  2:116492/251‑1 (MQ=255)                                                                                                   ttAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAACCAGCAATCGACGCTCTCCCGAGTCAGAAGACATTAAAGAGACTGAGCAAAAAGCATATCATTGATGATCACCTCAGAGAGTATGG‑CAGAAAAGAAAAACAAAAAATG‑ACCTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGGAAAAAATTATTGTGACCATTTTACCGTTGAATCAATaa                                                                                                                                         >  1:384632/1‑250 (MQ=255)                                                                                                     aGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGc                                                                                                                                      <  1:361458/251‑1 (MQ=255)                                                                                                     aGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGGGTCAAATGACATTAACGAGACTGAGCAAAAAGGATATCATTGATGTTCAACTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCCTTGTTTTTCAAACTGAAGTGTCAGCGAGAGGAAAAAAATATGGTGACCATTTTACCGGTGACTCAATAAGGc                                                                                                                                      >  1:63515/1‑251 (MQ=255)                                                                                                       cAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGTTGTTTACCTCAGAGAGGATGG‑CAGAACAGAAAAAAAAAAA‑TGCACCTTTGTTTTTTAACATGAAGTGTCAGCGAGAGTAAAAAATTATTGTGACCATTTTTACGTTGGCTCAAAAAGGc                                                                                                                                      >  1:504727/1‑249 (MQ=255)                                                                                                              ccGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCCATCATCCTATTTAGCCAGCAATCGCCGCTCTCCGGAGTCAGAAGACATTAACGAGACTGAGCAAAAAGCATATCATTGAAGTTCACCTCCAAGAGGATAG‑CAGAACCGAAAAACAAAAA‑TGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGCGAGGAAAAAATTATGGGGACCATTTTACCGATGACTAAAATAGGCAGATGgggg                                                                                                                             >  1:120743/1‑251 (MQ=255)                                                                                                                    cACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGTCTGAGCAAAGAGCATATCATTGATGTTCGCCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAg                                                                                                                       <  1:465184/251‑1 (MQ=255)                                                                                                                    cACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAg                                                                                                                       <  2:305723/251‑1 (MQ=255)                                                                                                                           ccAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaa                                                                                                                <  1:56929/251‑1 (MQ=255)                                                                                                                                       gCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAcg                                                                                                    <  2:354410/251‑2 (MQ=255)                                                                                                                                             aCCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGcc                                                                                               >  2:532998/1‑250 (MQ=255)                                                                                                                                                      gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCGTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                                     <  1:491468/251‑1 (MQ=255)                                                                                                                                                           caagcgtattcaaccAGCAATCGAAGCTCTACGCAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAGTGATGTTCACCTCAGATGTTATTG‑CAGACCAGAAAAACAAAAA‑TGCACCGTTTTTTTTCAACCGTAAGTGTCAGCGAGAGGAAAAAATTATGGTGACCATTTTACCGTTGCCTCAATAACGCCGATGTGGGAGGCAGCGATAAAACGTACCGGACTCCTCTATCGAAATTCatag                                                                                 <  2:439762/237‑2 (MQ=255)                                                                                                                                                             atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                                >  1:465132/1‑249 (MQ=255)                                                                                                                                                                    taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGAc                                                                       >  1:482287/1‑251 (MQ=255)                                                                                                                                                                           gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTa                                                                <  1:214946/251‑1 (MQ=255)                                                                                                                                                                                 gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                          >  1:94372/1‑251 (MQ=255)                                                                                                                                                                                       ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTgg                                                      <  1:333371/249‑1 (MQ=255)                                                                                                                                                                                       atcCGGAGTCAGAGGACAATAACGAGACTGAGCCAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAACATTATGGTGACCATTTTACCGTTGACTCCATAAGGCAGAAGTGGGACGCAGCGATAAAACGTGCAGGACTCAGCAATCGAAAATCATATCAGTCGAGAAATACTTATGCCTGCTGGTc                                                    >  2:448191/2‑251 (MQ=255)                                                                                                                                                                                         ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGc                                                  <  2:332096/251‑1 (MQ=255)                                                                                                                                                                                         ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAATATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGATTAAAAAATTATGGTGACCATTTTACCGTGGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                                    <  2:19922/249‑1 (MQ=255)                                                                                                                                                                                            gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTg                                                >  1:412721/1‑250 (MQ=255)                                                                                                                                                                                                 aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATAGCAGTCGAGATATACTTATGCCTGCTGGTCGCTGACAGCt                                          <  2:214479/251‑1 (MQ=255)                                                                                                                                                                                                           aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                 <  2:56445/250‑1 (MQ=255)                                                                                                                                                                                                                    gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAttt                        <  1:249775/250‑1 (MQ=255)                                                                                                                                                                                                                               atatCATTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc              >  1:8390/1‑249 (MQ=255)                                                                                                                                                                                                                                 atCATTGATGTTCACCTCAGAGTGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc              >  1:157627/1‑247 (MQ=255)                                                                                                                                                                                                                                    aTTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                         <  2:107552/233‑1 (MQ=255)                                                                                                                                                                                                                                    aTTGATGTTCACCTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                         >  1:107552/1‑233 (MQ=255)                                                                                                                                                                                                                                          gTTCACCTCAGAGAGTATGG‑CAGAACAGAANAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCGGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATg  <  2:265329/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGG‑CAGAACAGAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATG  >  NC_000913/1652092‑1652574

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