Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,415,898	T→C	*542Q (TAA→CAA)	aidB →	putative acyl‑CoA dehydrogenase AidB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,415,898	0	T	C	100.0%	213.3 / NA	58	*542Q (TAA→CAA)	aidB	putative acyl‑CoA dehydrogenase AidB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/25);  total (33/25)

CGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCG  >  NC_000913/4415761‑4416033                                                                                                                                           |                                                                                                                                        cgcAAATGTTAAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCaa                                                                                                                                       <  1:29292/140‑1 (MQ=255)      aTGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCaa                                                                                                                                       >  2:369734/1‑135 (MQ=255)         ttGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGta                                                                                                                                  >  2:279912/1‑137 (MQ=255)         ttGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGta                                                                                                                                  <  1:279912/137‑1 (MQ=255)                  gCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGc                                                                                                                      >  1:944859/1‑140 (MQ=255)                   cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCg                                                                                                                     >  2:952027/1‑140 (MQ=255)                   cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCg                                                                                                                     >  1:72719/1‑140 (MQ=255)                     tccccGCCAATGGCGCAGGCGTGGTTGCAGGGTATGGTAGATTCGCGCGGGGGGGGTCGGTTTGTAGAGCAGATCCAGGATGGTTTAATTCCGGGGGGGGAC‑GGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGAc                                                                                                                   <  2:661196/137‑1 (MQ=255)                      ctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACg                                                                                                                  <  1:291806/140‑1 (MQ=255)                      ctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACg                                                                                                                  <  2:944859/140‑1 (MQ=255)                           ccAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGt                                                                                                             >  2:301093/1‑140 (MQ=255)                                 gcgcAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGAt                                                                                                       <  1:73229/140‑1 (MQ=255)                                   gcAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTa                                                                                                     <  2:897825/140‑1 (MQ=255)                                          tggtgTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAAt                                                                                              >  2:565337/1‑140 (MQ=255)                                              gtCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGaaa                                                                                          <  2:72719/140‑1 (MQ=255)                                                      aTGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaa                                                                                  >  1:356839/1‑140 (MQ=255)                                                      aTGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaa                                                                                  >  2:70507/1‑140 (MQ=255)                                                       tGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                 >  1:811164/1‑140 (MQ=255)                                                        gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGACTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaaa                                                                                >  1:916300/1‑140 (MQ=255)                                                                 gcgcgGTGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGcc                                                                       >  2:964814/1‑140 (MQ=255)                                                                 gcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGcc                                                                       >  1:671051/1‑140 (MQ=255)                                                                 gcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGcc                                                                       >  1:209287/1‑140 (MQ=255)                                                                   gcggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCgg                                                                     >  2:401707/1‑140 (MQ=255)                                                                    cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGACTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGa                                                                    <  2:916300/140‑1 (MQ=255)                                                                                 ttGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTcgcg                                                       >  1:631070/1‑140 (MQ=255)                                                                                   gTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgcc                                                     >  2:840039/1‑140 (MQ=255)                                                                                   gTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgcc                                                     <  2:811164/140‑1 (MQ=255)                                                                                        agCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATc                                                >  1:901182/1‑140 (MQ=255)                                                                                         gCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATcc                                               >  2:865447/1‑140 (MQ=255)                                                                                           aGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCgg                                             >  2:258157/1‑140 (MQ=255)                                                                                           aGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCgg                                             <  2:631070/140‑1 (MQ=255)                                                                                           aGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCgg                                             <  1:301093/140‑1 (MQ=255)                                                                                             aTCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATc                                                <  2:76087/135‑1 (MQ=255)                                                                                             aTCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATc                                                >  1:76087/1‑135 (MQ=255)                                                                                             aTCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCa                                           >  1:320098/1‑140 (MQ=255)                                                                                                  gAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGtt                                      >  1:371195/1‑140 (MQ=255)                                                                                                   aaTGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTc                                     >  2:28298/1‑140 (MQ=255)                                                                                                    aTGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgcc                                                     >  2:272274/1‑123 (MQ=255)                                                                                                    aTGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgcc                                                     <  1:272274/123‑1 (MQ=255)                                                                                                      gATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAtt                                  <  2:320098/140‑1 (MQ=255)                                                                                                       aTTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCa                                                  >  2:147107/1‑123 (MQ=255)                                                                                                       aTTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCa                                                  <  1:147107/123‑1 (MQ=255)                                                                                                       aTTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAttt                                 >  2:499442/1‑140 (MQ=255)                                                                                                       aTTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAttt                                 >  2:290289/1‑140 (MQ=255)                                                                                                                ctgcGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGa                         >  1:311123/1‑140 (MQ=255)                                                                                                                  gcGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATg                       <  2:356839/140‑1 (MQ=255)                                                                                                                  gcGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATg                       >  2:618106/1‑140 (MQ=255)                                                                                                                   cGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGa                      >  2:751981/1‑140 (MQ=255)                                                                                                                   cGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGa                      >  2:506868/1‑140 (MQ=255)                                                                                                                   cGGGCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGa                      <  1:401707/140‑1 (MQ=255)                                                                                                                   cGGGCGACGGCGGGAGTGTGAGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGa                      >  2:661692/1‑140 (MQ=255)                                                                                                                     ggCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                    <  2:209287/140‑1 (MQ=255)                                                                                                                     ggCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                    <  1:865447/140‑1 (MQ=255)                                                                                                                     ggCGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGaca                    >  1:635032/1‑140 (MQ=255)                                                                                                                       cGACGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACt                  <  1:565337/140‑1 (MQ=255)                                                                                                                          cGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCAt               <  1:751981/140‑1 (MQ=255)                                                                                                                          cGGGGGGAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCAt               <  1:506868/140‑1 (MQ=255)                                                                                                                              gggAGTTTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCAt           <  2:173996/140‑1 (MQ=255)                                                                                                                              gggAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCAt           <  1:952027/140‑1 (MQ=255)                                                                                                                              gggAGTGTGTGTGCAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCAt           <  1:70507/140‑1 (MQ=255)                                                                                                                                   tgtgtgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCa            <  2:32052/134‑1 (MQ=255)                                                                                                                                   tgtgtgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCa            >  1:32052/1‑134 (MQ=255)                                                                                                                                    gtgtgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTa     >  1:56936/1‑140 (MQ=255)                                                                                                                                     tgtgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTAt    >  2:487684/1‑140 (MQ=255)                                                                                                                                      gtgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATc   >  1:183126/1‑140 (MQ=255)                                                                                                                                       tgtgcAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCg  >  2:489206/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                        CGCAAATGTTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTT‑ATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCG  >  NC_000913/4415761‑4416033

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