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%	151.0 / NA	41	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/20);  total (21/20)

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATT  >  NC_000913/1652087‑1652551                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGCCTGAGCCAAAAG‑CATATCATTGATGTTCACCTCag                                                                                                                                                                                                                            >  2:13263/1‑249 (MQ=255)   ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCagaga                                                                                                                                                                                                                         <  1:67815/250‑1 (MQ=255)                    gctcttcgatctccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                                                                                                                                                       <  2:673457/239‑1 (MQ=255)                                ccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTTTTTAGCCAGCAATCGACGCTCTCCGGCGGCAGATTACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCCCCTCAGAGAGGATGGCCGAACAGaaaa                                                                                                                                                                                                       >  1:673457/1‑239 (MQ=255)                                   ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                                                                                                                                                        <  1:91423/251‑1 (MQ=255)                                                     ggATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTTAACC‑TGAAgg                                                                                                                                                                     >  1:620745/1‑250 (MQ=255)                                                               aCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAAACGGGTCATTCCTATTATTAAGCCAGCAAACGACGCTCTCCGGAGTCAGATGACATTAACGGGACTGAGCAAAAAA‑CAAATCATTGATGTTTAACTCCGCGGGGATGGCAGAAAAGAAAAAAAAAAAAGAACCTTTGTTTTTTAACC‑TGAAGTGTCAGCgggg                                                                                                                                                           >  1:271170/1‑248 (MQ=255)                                                                         gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTaaaaaa                                                                                                                                                   <  1:726227/249‑1 (MQ=255)                                                                             aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGAACAAACGGGTCATTCATATTATTAAACCAGCAACCGACGCCCTCCGGAGTCAGGTGACATTAACGAGAATGAGCAAAGAG‑CATATAATTTATTTTCAACTCAGAGAGGATGGCAGAACAGAAAAAAAAAAATGCCCCTTTGTTTTTCTACC‑TGAAATGTCCGCGAGAGTaaaaaaaaat                                                                                                                                               >  2:513691/1‑245 (MQ=255)                                                                             aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGGGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑Tgg                                                                                                                                            >  1:50401/1‑251 (MQ=255)                                                                              aTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGt                                                                                                                                           <  1:13263/251‑1 (MQ=255)                                                                                  aaCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAGCCAGCAATCGCCGCTCTCCGCAGGCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCATTGTTTTTTAAAC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGAcc                                                                                                                                       >  1:485027/1‑251 (MQ=255)                                                                                     aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCAtt                                                                                                                                    <  1:458742/251‑1 (MQ=255)                                                                                            atCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                                                                                                                                                     <  1:324134/197‑1 (MQ=255)                                                                                                     aaTTTAGCAAAAGACCGTTTCACGGTACCACAAACACAGGCGGGAACCAAACGGGGCATTCATCTTATTAAACCAACAAAAGACGCTCACCCGAGGACGATGACATTAACGAGACTGAACAAA‑AGACCTATCATTTATGTTCACCTCACAGAGTATGGCAGAAACGAAAAACAAAAAAGCACCTTTGTTATTCAACC‑TGAAATGTCAGCGAGAGAAAAAAATTAATGG‑GACCATTTTATCGTTggatcaat                                                                                                                    >  2:232275/1‑244 (MQ=255)                                                                                                              aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCACACGGGGCCTTCCTCTTTTTAAGCCAGCAATCCACGGCCCCCGGAGTCAGATGACATTAAAAAGACTGTGAAAAAAG‑CATATAAATGATGTTCACCTCAGAGAGTATGGGAGAAAAAAAAAAAAAAAATGCACCTTTGTTTTTCAACCCTGAA‑TATCAGCGAGAGTAAAAAAATA‑AGGGGAACATTTTACCCGTGACTCaaaaaagaaaa                                                                                                            >  2:761157/1‑241 (MQ=255)                                                                                                              aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGAc                                                                                                                                        <  1:332654/222‑1 (MQ=255)                                                                                                                aGACCGTTTCACGGTACTCAAAACGTATGCGGGAACCAATCTGGTCATTCATCTTATTAAGCCAGTCATCGACGCTATCCGTAGTCAGATTACATTAACGTGGCTGAGCAAAGAG‑CTTAGTATTGATGTTCGTCTCAGAGAGTATGGCAGAACAGAAAATGAAAAATGCATCTTTGTTTTTAAACC‑TGAATTATCAGCAAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATTAGGCACAtgt                                                                                                         <  2:512101/251‑1 (MQ=255)                                                                                                                  cccGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTgg                                                                                                       <  2:529077/250‑1 (MQ=255)                                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTAACGTTGACTCAATAAGGCAGATGTGGGa                                                                                                     <  2:46551/251‑1 (MQ=255)                                                                                                                                      acacAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGc                                                                                   <  1:91543/251‑1 (MQ=255)                                                                                                                                        ccaggctgtaaCCAATCGTGTCATTCATCTTNTTAAGCCAGCAAGCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATAGCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCg                                                                                 <  2:355181/242‑1 (MQ=255)                                                                                                                                            gCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACt                                                                             >  1:691991/1‑251 (MQ=255)                                                                                                                                                gAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGc                                                                         >  2:95128/1‑251 (MQ=255)                                                                                                                                                  aCCAAACGGATCATTCATATTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCACCC‑TGAAGTGTCAGCGTGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTTACTCAATAAGGCAGATGTGGNACGCAGCGTTAAAACGTGCCGGACCCCgacc                                                                       >  2:339382/1‑248 (MQ=255)                                                                                                                                                         gagTCATTTATCTTAGTAAGCCAGCAATCGACGCCATCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAGAGAAAAACAAAAATGCACATTTGTTTTTCAACC‑TGAAGAGTAAGCGAGAGTAAAAAATTC‑TGGTTACCATTTTACCGTAGACTCAATAAGGAAGATGTGGGACGCAGCGATAAAACGTGCGGGCCTCCGCCATCGaac                                                                 >  2:531341/3‑249 (MQ=255)                                                                                                                                                             cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCat                                                            >  1:72663/1‑251 (MQ=255)                                                                                                                                                                      tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                    <  2:50401/250‑1 (MQ=255)                                                                                                                                                                         taAGCCAGCAATCGAAGCTCTCCGGAGTCANATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTCATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACGCCGCCATCGAAAATCATATCAGTCgaga                                                 <  2:373099/250‑1 (MQ=255)                                                                                                                                                                                  aaTCGACTCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAATTAtc                                       >  1:21856/1‑250 (MQ=255)                                                                                                                                                                                   aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGc                                      <  1:730310/251‑1 (MQ=255)                                                                                                                                                                                         gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAACGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                >  2:596813/1‑251 (MQ=255)                                                                                                                                                                                               cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCt                          <  1:95128/251‑1 (MQ=255)                                                                                                                                                                                                ggAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCACCTGGTCGCTg                         >  2:468216/1‑251 (MQ=255)                                                                                                                                                                                                     cAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGGCTCCGCCATCGAAAATCATATCAGTCGNGACATACTTATGCCTGCTGGTCGCTGACAGc                    >  2:337644/1‑251 (MQ=255)                                                                                                                                                                                                     cAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGc                    <  1:411564/251‑1 (MQ=255)                                                                                                                                                                                                         tGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGt                >  1:117596/1‑251 (MQ=255)                                                                                                                                                                                                                   gagaCTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCgg      <  1:513691/251‑1 (MQ=255)                                                                                                                                                                                                                       cTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt  >  1:140225/1‑251 (MQ=255)                                                                                                                                                                                                                            cAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCATCC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCa    <  2:403676/244‑1 (MQ=255)                                                                                                                                                                                                                            cAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCa    >  1:403676/1‑244 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACC‑TGAAGTGTCAGCGAGAGTAAAAAATTA‑TGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATT  >  NC_000913/1652087‑1652551

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