Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,673,048	G→T	M334I (ATG→ATT)	hcaD →	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,048	0	G	T	100.0%	223.7 / NA	63	M334I (ATG→ATT)	hcaD	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (38/25);  total (38/25)

CGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATC  >  NC_000913/2672913‑2673183                                                                                                                                        |                                                                                                                                        cGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTg                                                                                                                                     <  1:1178381/140‑1 (MQ=255)                 gCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTg                                                                                                                     <  2:1388330/140‑1 (MQ=255)                          cAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCaa                                                                                                            <  1:1197492/140‑1 (MQ=255)                                               gcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAt                                                                                       <  2:1176142/140‑1 (MQ=255)                                               gcAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGAt                                                                                       <  2:293353/140‑1 (MQ=255)                                                    gTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                  <  2:1422888/140‑1 (MQ=255)                                                       aGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                  <  2:1587594/137‑1 (MQ=255)                                                       aGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGt                                                                                  >  1:1587594/1‑137 (MQ=255)                                                         ggCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAAt                                                                             >  2:272997/1‑140 (MQ=255)                                                             accgctaccgctACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGc                                                                         >  2:762948/1‑140 (MQ=255)                                                             accgctaccgctACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGGTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGc                                                                         >  2:1804164/1‑140 (MQ=255)                                                              ccgctaccgctaCTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGTTATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCa                                                                        >  2:699585/1‑140 (MQ=255)                                                               cgctaccgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCaa                                                                       >  2:767733/1‑140 (MQ=255)                                                                   accgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                                                                   <  1:1334331/140‑1 (MQ=255)                                                                   accgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                                                                   >  1:1624767/1‑140 (MQ=255)                                                                   accgctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACg                                                                   <  2:1624767/140‑1 (MQ=255)                                                                       ctacTGCCGCCGCCGTGGTTCTGGAGCGATCAGGACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTAGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGt                                                               >  2:1179110/1‑140 (MQ=255)                                                                           tGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCtt                                                           >  2:1470756/1‑140 (MQ=255)                                                                            gccgccgccGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTa                                                          >  2:1482312/1‑140 (MQ=255)                                                                             ccgccgccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTAt                                                         <  1:1482312/140‑1 (MQ=255)                                                                                 cgccgTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGt                                                     <  1:1119714/140‑1 (MQ=255)                                                                                     gTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCgg                                                 <  1:922889/140‑1 (MQ=255)                                                                                      tGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGt                                                <  1:767733/140‑1 (MQ=255)                                                                                            tGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCt                                          <  1:1552813/140‑1 (MQ=255)                                                                                              gAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAACCTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGa                                        <  1:1804164/140‑1 (MQ=255)                                                                                                   aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGTCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                    >  1:485396/1‑139 (MQ=255)                                                                                                   aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGTCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                    <  2:485396/139‑1 (MQ=255)                                                                                                   aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAg                                   <  1:699585/140‑1 (MQ=255)                                                                                                   aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAg                                   >  1:590656/1‑140 (MQ=255)                                                                                                   aTCAGTACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAg                                   >  1:1577698/1‑140 (MQ=255)                                                                                                        tACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                    <  2:1666691/134‑1 (MQ=255)                                                                                                        tACAGTGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCa                                    >  1:1666691/1‑134 (MQ=255)                                                                                                             tGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGTAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGt                             >  1:1041201/1‑136 (MQ=255)                                                                                                             tGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGTAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGt                             <  2:1041201/136‑1 (MQ=255)                                                                                                             tGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTgaga                         >  2:1249208/1‑140 (MQ=255)                                                                                                             tGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTgaga                         <  1:272997/140‑1 (MQ=255)                                                                                                             tGATAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTgaga                         >  2:1602679/1‑140 (MQ=255)                                                                                                               aTAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGAtt                       >  2:1380879/1‑140 (MQ=255)                                                                                                               aTAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGAtt                       >  2:562150/1‑140 (MQ=255)                                                                                                               aTAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGAtt                       >  1:147622/1‑140 (MQ=255)                                                                                                               aTAACTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGACGTGAGAtt                       >  2:1636582/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                    >  1:1016598/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                    >  2:728986/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                    >  2:250317/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                    >  1:778215/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGc                    >  2:24668/1‑140 (MQ=255)                                                                                                                  aCTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTAAATCAGGGGCGTGAGATTCGc                    >  2:933757/1‑140 (MQ=255)                                                                                                                   cTTACAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCTGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGcc                   <  1:1636582/140‑1 (MQ=255)                                                                                                                       cAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAAt               >  2:91816/1‑140 (MQ=255)                                                                                                                       cAGTTTATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAAt               >  1:1283412/1‑140 (MQ=255)                                                                                                                          tttATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCg            >  1:899253/1‑140 (MQ=255)                                                                                                                          tttATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCg            <  2:590656/140‑1 (MQ=255)                                                                                                                            tATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCa          >  1:724771/1‑140 (MQ=255)                                                                                                                            tATTGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCa          <  1:1380879/140‑1 (MQ=255)                                                                                                                              ttGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCaaa        <  1:91816/140‑1 (MQ=255)                                                                                                                               tGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAAtt              <  1:752711/133‑1 (MQ=255)                                                                                                                               tGGCGATATTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAAtt              >  2:752711/1‑133 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  2:1211524/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  2:1253014/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  1:1056930/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  1:95965/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  2:1391582/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  1:1633761/1‑140 (MQ=255)                                                                                                                                    atatTCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATc  >  2:1507385/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        CGAAAGCTGGGAAAACGCCAATAACCAGGCGCAAATTGCCGCTGCCGCAATGTTAGGGCTACCGCTACCGCTACTGCCGCCGCCGTGGTTCTGGAGCGATCAGTACAGTGATAACTTACAGTTTATTGGCGATATGCGTGGCGATGACTGGCTTTGTCGTGGCAACCCGGAAACTCAGAAGGCGATTTGGTTTAATCTGCAAAACGGCGTGCTTATCGGTGCGGTCACGCTGAATCAGGGGCGTGAGATTCGCCCAATTCGCAAATGGATC  >  NC_000913/2672913‑2673183

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