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

TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCC  >  NC_000913/1652086‑1652545                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                        tATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCag                                                                                                                                                                                                                     <  1:214545/250‑1 (MQ=255)     gtgtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGTCATTCATCTTATTAAGCCAGCACTCG‑ACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGt                                                                                                                                                                                                                >  1:576566/1‑251 (MQ=255)               gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                                                                                                                                                         <  1:66843/248‑1 (MQ=255)                          aGTTCGCCCTGGGGAGCTGTGTGCACTTGGATGGGAGGACATAGAGCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTGTCACGGTACCAAAAACACAGGCGGGCACCAAGCGGGTCATTCATCTTATTAAGCCAACAATCG‑ACGCTCTCCGGAGTCAGATGTCATTAAAGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaacccaat                                                                                                                                                                                           <  2:628777/251‑7 (MQ=255)                               gCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGc                                                                                                                                                                                     >  1:84194/1‑251 (MQ=255)                                    ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCAcc                                                                                                                                                                                  <  2:632883/249‑1 (MQ=255)                                        aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCCGCAATCG‑ACTCTCTCCGGAGGCAAAGGTAATTAACCGGACTGAGCAAAGAGACA‑ATCTTTTATGTTAACCTCAGAGAGTATGGGAGAAAAGAAAAA‑CAAAAATGAACCTTTGtt                                                                                                                                                                            >  1:232315/1‑251 (MQ=255)                                          cTGTGTGCACTTGGATGGGAGGACATATATCTGAAAAATGGAAACATGATGATCAGGAGAAATTTCGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCACTCAGGGCATTAATCTTATTTAACCAGCAATCGGACG‑TCTCCGGAGCCAGATGACATTAACGAGGATGAGCAAAGAG‑CATATCATTGATGGTCACCTCCGAGAGGCTGGCCGAACAGAAAAAACaaaaagccactttttttt                                                                                                                                                                           >  2:205718/1‑238 (MQ=255)                                                  aCTTGTATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAATAACACAGGCGGGAGCCAAGCTGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑GATATCAGTGATGTTGACCTCAGAGAGTATGTCAAAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                                                                   <  2:412240/250‑1 (MQ=255)                                                            gaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                                                        <  2:543880/251‑1 (MQ=255)                                                             aggaCATAGAGCTGAAAAATGGAACAATGATGAGCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAAGCG‑ACGCTCTCCGGAGTCAGATGACAGTAACGATACTGAGCCAAGAG‑CATATCATTGATGTTCATCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCATACTTAAGTGTCAGCg                                                                                                                                                       <  2:441878/251‑1 (MQ=255)                                                                    aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCATCGAGAGTa                                                                                                                                                 <  1:147366/250‑1 (MQ=255)                                                                          gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACATAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                               <  1:466436/246‑1 (MQ=255)                                                                          gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                                                                                                                                              <  1:376094/247‑1 (MQ=255)                                                                            aaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTAt                                                                                                                                        <  1:163569/251‑1 (MQ=255)                                                                            aaaaTGAAACAAAGATGAGCATGAGAAAAGTAGCAAATGAACGTGACAAGGTACCAAAAACACAGGCGGGAACCAATCGCGTCATTCATCTTATTAAGACAGCAAAAG‑ACGCTCACCGGAGTAAGAAGACATTATAGAGACTGATCAAAAAG‑AATATCATTGATGATCAACTCAGAGAGTATGGCAGAAAAGAAAAA‑CAAAAAAGCACCATAGTTTTTAAACATGAAGCGACAGCAAGAGTAAACAATTAt                                                                                                                                        <  2:454167/251‑1 (MQ=255)                                                                                        tgatgaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTTAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtttt                                                                                                                            >  1:461079/1‑251 (MQ=255)                                                                                          atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCCTCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAA‑AAAAAATGCAc                                                                                                                                                                                   >  1:182485/1‑194 (MQ=255)                                                                                          atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCAc                                                                                                                                                                                   <  2:182485/194‑1 (MQ=255)                                                                                          atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAAACG‑ACGCTCTCCGGAGTCCGATGACAATAAAGAGACTGAGCAAAAAA‑CCTATCATTGATGTTCACCTCACAGAGGATGGCAGAAAAGAAAAA‑CAACAATGAAACTTTGTTTTTCAACCTGAAATGGTAGGGCGAGTAAAAAATTATGGGGAACAATTTAc                                                                                                                          >  2:138962/1‑251 (MQ=255)                                                                                            gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATATTATTAAGCCAACAATCG‑CCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAA‑CATATCATTGATGTTCACCCCAGAGAGAATGGCAGAACAAAAAAA‑CAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGAGAGTAAAAAATTAATGTGACCATTTTAccc                                                                                                                        >  2:187784/1‑250 (MQ=255)                                                                                               cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGttt                                                                                                                     >  1:131827/1‑250 (MQ=255)                                                                                               cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTg                                                                                                                     <  1:485464/251‑1 (MQ=255)                                                                                                  gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                               <  1:278126/222‑1 (MQ=255)                                                                                                  gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTAAGGCCGAACAGAAAAA‑AAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                               >  2:278126/1‑222 (MQ=255)                                                                                                   agaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCAATCATCTTATTAAGCCAGCAATCG‑ACGCCCTCCGGAGTCAGATTAGATTAACGAGAATGAGCCAAAAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAA‑CAAAAATGATCCTTTGTTTTTCAACCTGAAATGTCAACGAGAGAAAAAAATTGAGTTGACCACTTTACCCTTTACTc                                                                                                                 >  1:601825/1‑251 (MQ=255)                                                                                                                aaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTTATCTTATTAAGCCAGCAATCG‑ACGCTCTCCCGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAAAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAATGTCAGCGAGAGTAAAAAATTATGGTGACCCTTTTACCGTTTACTCAATa                                                                                                             >  2:11835/1‑242 (MQ=255)                                                                                                                aaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATa                                                                                                             <  1:11835/242‑1 (MQ=255)                                                                                                                         tCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGGTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGc                                                                                            <  2:576566/250‑1 (MQ=255)                                                                                                                           aCGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGAAATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTTACCTCCGAGAGTATGGCAGAAAAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGGGAGGAAAAAAATATGGGGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGc                                                                                         >  1:275598/1‑251 (MQ=255)                                                                                                                                aCCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaa                                                                                    <  2:601825/251‑1 (MQ=255)                                                                                                                                      aaCACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCCGCGGTAAAACGTg                                                                              >  1:276501/1‑251 (MQ=255)                                                                                                                                                      aaTCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCg                                                              >  2:344019/1‑251 (MQ=255)                                                                                                                                                               attcatCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                    >  1:60222/1‑172 (MQ=255)                                                                                                                                                               attcatCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                    <  2:60222/172‑1 (MQ=255)                                                                                                                                                               attcatCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCata                                                     <  2:93708/251‑1 (MQ=255)                                                                                                                                                                  cattttattaAGCCAGCAATCG‑ACGCTCTCCTGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACATAAAAA‑CAAAAATGCGCCTTTGTTTTTCAACCTGAAGTGTCAGCGTGAGTAAAAAAGTATGGTGACCATTTTTCCGTTGACTCAATAAGGCAGATTGGGGACGCAGCGATAAAACGTGCCTGACTCCTCCAGCTAAAATCATATc                                                   <  2:538843/246‑1 (MQ=255)                                                                                                                                                                     cTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTc                                               >  1:386762/1‑251 (MQ=255)                                                                                                                                                                        attaAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCga                                             <  2:84194/250‑1 (MQ=255)                                                                                                                                                                               cAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACt                                     <  1:507522/251‑1 (MQ=255)                                                                                                                                                                                         gcGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                            <  2:71365/250‑1 (MQ=255)                                                                                                                                                                                         aCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                            >  2:569878/1‑251 (MQ=255)                                                                                                                                                                                                ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGc                     <  1:344019/251‑1 (MQ=255)                                                                                                                                                                                                     gTCAGATGACCTTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGAc                 >  1:583431/1‑250 (MQ=255)                                                                                                                                                                                                             aCATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGc        >  1:231433/1‑251 (MQ=255)                                                                                                                                                                                                              cATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCt       <  1:493080/251‑1 (MQ=255)                                                                                                                                                                                                               aTTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCt       >  2:183878/1‑250 (MQ=255)                                                                                                                                                                                                                   aCGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAccc  >  1:473967/1‑251 (MQ=255)                                                                                                                                                                                                                    cGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAccc  <  1:178899/250‑1 (MQ=255)                                                                                                                                                                                                                    cGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAccc  >  2:178899/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                        TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCC  >  NC_000913/1652086‑1652545

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: ‑