Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	364,170	2 bp→AA	coding (2135‑2136/3075 nt)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,170	0	G	A	100.0%	123.4 / NA	39	A712V (GCC→GTT)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)
*	NC_000913	364,171	0	G	A	100.0%	125.6 / NA	39	A712V (GCC→GTT)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)

GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGAT  >  NC_000913/363927‑364410                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                 gtctgtATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGTCGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTg                                                                                                                                                                                                                                            <  1:66450/246‑1 (MQ=255)   cGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑Ac                                                                                                                                                                                                                                         <  1:16674/251‑1 (MQ=255)                    cGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCgg                                                                                                                                                                                                                        <  1:5465/251‑1 (MQ=255)                                                  agTTGTTTTTTATCGCCAAGCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATGGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCGGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTGCCGCCAGACGTCGCTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACGACGCGTACTGTGAGCCAGAGTTg                                                                                                                                                                                          <  2:335419/251‑1 (MQ=255)                                                      gttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCg                                                                                                                                                                                      <  1:222934/251‑1 (MQ=255)                                                                         caccTGTGAAAGAAAGCATGACTGGCGGTTAAATTGCTATCGCTTAGTGCCCATCTCGAGGCAATAATCCATGGCGCTGGTGGTCAGATGAGGGATGGCGCTGGACTCGGCGGGGATCGTCACACTGAGGGTTGCCGGCAGTCGCCACTGCTGCCAGGCGCGGAGGTGCCCAACTGCTG‑TCCATGCGGTCGCGTTCGGTTGCTCTACGCGAACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCGGCGGtc                                                                                                                                                                   <  2:60136/248‑2 (MQ=255)                                                                         caTCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTa                                                                                                                                                                   <  1:138441/251‑1 (MQ=255)                                                                                 aaagaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGGCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGc                                                                                                                                                           <  1:219755/251‑1 (MQ=255)                                                                                                       aaaTTGCCAACGCTTATTACCCAGCTCGTTGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGGACGTCACACTGAGGGTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCCTGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGGAGGTCAGGCAGTTCACTCAACTGTTTACCtt                                                                                                                                     >  1:73180/1‑251 (MQ=255)                                                                                                            gCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGa                                                                                                                                <  2:70180/251‑1 (MQ=255)                                                                                                                    ttattatacAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGGTCGTCACCATGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTCCTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGGTTACCTTGCGGCGCGACATc                                                                                                                        >  1:93983/9‑251 (MQ=255)                                                                                                                          cccAGCTCGATGCAAAATTCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCTGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGTCGCTCTCCGGCTGCGGTAGTTCAGGCATGTCAATCTACTGTTTACCTTGTGGAGCGACATCCAGAgg                                                                                                                  <  1:154621/251‑1 (MQ=255)                                                                                                                                  gATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGGGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCCGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCAc                                                                                                          >  1:287657/1‑251 (MQ=255)                                                                                                                                   aTGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGCTGGCGTGGGACGCGGCGGGGGGCGTCCAACTGAGGTTTTCCGCCCGACGCCACTGCTGCCAGGCGCTTATGTGCCCAACTTCTGTACC‑TGCGGTCGCGTTCGGTTGCCCTACGCGGCCTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCCGGCAGTTCAATCAACCTTTTACCTTGTGGAGCGGCATCCAGAGGCCCTTcccc                                                                                                         >  2:21594/1‑248 (MQ=255)                                                                                                                                   aTGCAAAAATCCATTTCGCTGGTGGGCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGATGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCGCCGGCTGCGGTAGTTCAGGCAGTTCAAGCAACTGTGTACCTTGTGGAGCTACATCCAGAGGCACTTCAcc                                                                                                         <  2:160877/251‑1 (MQ=255)                                                                                                                                      cAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGACGTGGGACGCGACGGGCAGCGTCACACTGAGGTTTTCCGCCAGAAGACCCTGCTGACAGGCGCTGACGTGACCAACTTCCG‑AACATGCCGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGGTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCCATCAACTGTTTACCTTGTGGCGCGACATCCAGAGGCACTTCACCGCt                                                                                                      >  1:357653/1‑251 (MQ=255)                                                                                                                                               aTTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGAGAGCGTCACACTGAGGGTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTTCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCCGAGGCACTTCACCGCTTGCCAGCgg                                                                                             >  1:251026/1‑251 (MQ=255)                                                                                                                                                       ggtggtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCa                                                                                      <  1:121871/250‑1 (MQ=255)                                                                                                                                                         gggtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGGCACACCGGGGTTTTCCGCCAGACGCCACTGCTGACAGGCGCTGGTATGCCCAACTTCTG‑ACCCTGCGGGCGCGTTCGGGTGCAATACGCGTACGGAGAGCCCGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGGTCCATCCACTGGGTACCTTGTGGCGCGACATCCAGAGGCACTTCACCGCtatcacgcggcttaacaacc                                                                                   >  2:286398/2‑232 (MQ=255)                                                                                                                                                          ggtCAGAAGCGGGATGGCGTGGGACGCGGCGAGCAGCGTAACACTAAGGGTTTCCGCCAGACCCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCCCTACCCGTACCGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGTGGTAGTTCAGGCAGTTCAATCAACCGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCCGCGGCTTACCATcc                                                                                   >  1:301630/1‑250 (MQ=255)                                                                                                                                                           gtCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAg                                                                                 <  1:284975/251‑1 (MQ=255)                                                                                                                                                                aTGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCAAACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCCCTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGGTCAATCAACTGTTTa                                                                                                                                         >  1:240331/1‑190 (MQ=255)                                                                                                                                                                 gtcgtgATGGCGTGGGACGCGGCGGGGAGCGTCACTCTGAAGTTTTCCGCCAGACGCCACTGCTGCCAGGCGTTGATGTGCCCAACTTCTG‑CCCATGCGGTCGCGTTCGGTTGCACTTCGCGTACTGTGAGCCAGAGTTGCCCTGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCACCTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCGCCGCTTGCCAGCGGCTTACCATCCAGCGcca                                                                            <  2:195032/245‑1 (MQ=255)                                                                                                                                                                     ggATGGCGTGGGACGCGGCGGGGAGCGTCACACTGCGGGTTTCCGCCAGACGCCACTGCTGCCAGGAGCTGATGTGCCAAAATTCTG‑CCCATGCCGTCGCGGTCGGCTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGGAGTTCAGGCAGGTCAATCAACTGTTTACCTTGTGGAGGGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATc                                                                       >  1:233145/1‑251 (MQ=255)                                                                                                                                                                     ggATGGCGTGGGACGCGGCGGGGAGCGTACCAACGAGGTTTTCCGCCCGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTCCGCCTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTTCCTTGTGGGGCGACATCCAGGGGCACTTCACCGCTTGCCAGCGGCTTACCATCCACCGCCACCATc                                                                       >  2:20810/1‑251 (MQ=255)                                                                                                                                                                       aTGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCa                                                                     <  1:104259/251‑1 (MQ=255)                                                                                                                                                                       aTGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCa                                                                     <  1:110626/251‑1 (MQ=255)                                                                                                                                                                              gggACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTGCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAgg                                                              <  2:216004/251‑1 (MQ=255)                                                                                                                                                                                gACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAg                                                            <  1:248808/251‑1 (MQ=255)                                                                                                                                                                                        ggggAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGCGCTCGTTAt                                                    >  2:302698/1‑251 (MQ=255)                                                                                                                                                                                                tCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTTAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTAt                                              <  1:261305/249‑1 (MQ=255)                                                                                                                                                                                                  acacTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCa                                                                                                                                                          <  1:347904/139‑1 (MQ=255)                                                                                                                                                                                                  acacTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCa                                                                                                                                                          >  2:347904/1‑139 (MQ=255)                                                                                                                                                                                                                        gCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACt                    <  2:76091/251‑1 (MQ=255)                                                                                                                                                                                                                          cACTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTc                  >  1:172322/1‑251 (MQ=255)                                                                                                                                                                                                                           aCTGCTGCCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCg                 >  1:158754/1‑251 (MQ=255)                                                                                                                                                                                                                              gctgcCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGCACAGGTATTCGCTGGTCCCTTCgctg              >  2:348266/1‑248 (MQ=255)                                                                                                                                                                                                                               ctgcCAGGCGCTGATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATgg             >  1:314191/1‑251 (MQ=255)                                                                                                                                                                                                                                           gATGTGCCCAACTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTCGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGAt  >  1:238384/1‑250 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                 GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTG‑ACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGAT  >  NC_000913/363927‑364410

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