Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,469	G→A	L141L (CTG→CTA)	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,672,469	0	G	A	100.0%	231.4 / NA	65	L141L (CTG→CTA)	hcaD	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/34);  total (31/34)

AAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCAC  >  NC_000913/2672336‑2672602                                                                                                                                      |                                                                                                                                      aaaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGccc                                                                                                                                 >  2:728361/1‑140 (MQ=255)   aGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGa                                                                                                                               <  1:579433/140‑1 (MQ=255)        gCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGt                                                                                                                          <  2:1131845/140‑1 (MQ=255)        gCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGt                                                                                                                          <  1:424351/140‑1 (MQ=255)           cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAg                                                                                                                       <  1:1438717/140‑1 (MQ=255)                 tCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCAGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGa                                                                                                                 <  2:917459/140‑1 (MQ=255)                  cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGt                                                                                                                          <  1:510684/130‑1 (MQ=255)                  cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGt                                                                                                                          >  2:510684/1‑130 (MQ=255)                  cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGAt                                                                                                                >  2:1603489/1‑140 (MQ=255)                      tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATc                                                                                                            <  2:1295514/140‑1 (MQ=255)                      tttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATc                                                                                                            <  2:652126/140‑1 (MQ=255)                              gCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCgg                                                                                                    <  1:807840/140‑1 (MQ=255)                               cAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                                                                   >  2:397742/1‑140 (MQ=255)                               cAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                                                                   >  1:269645/1‑140 (MQ=255)                                aaCCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGaa                                                                                                  <  1:1353893/140‑1 (MQ=255)                                  ccGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACt                                                                                                >  1:771498/1‑140 (MQ=255)                                    ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                                                              <  2:1530815/140‑1 (MQ=255)                                    ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                                                              <  2:164396/140‑1 (MQ=255)                                    ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                                                              >  1:180735/1‑140 (MQ=255)                                      cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTg                                                                                            >  1:1252393/1‑140 (MQ=255)                                                cGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                                                              >  1:712967/1‑128 (MQ=255)                                                cGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                                                              <  2:712967/128‑1 (MQ=255)                                                      cTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                                                                            <  1:1343879/140‑1 (MQ=255)                                                      cTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                                                                            <  2:269645/140‑1 (MQ=255)                                                      cTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                                                                            <  1:1603489/140‑1 (MQ=255)                                                              gCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCa                                                                    <  1:874568/140‑1 (MQ=255)                                                                   aTGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                                                               >  2:1193920/1‑140 (MQ=255)                                                                     gCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                                                             <  1:380798/140‑1 (MQ=255)                                                                     gCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                                                             <  1:941755/140‑1 (MQ=255)                                                                     gCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                                                             <  1:1548742/140‑1 (MQ=255)                                                                      cACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgc                                                            <  1:481737/140‑1 (MQ=255)                                                                         tGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa                                                         <  1:507796/140‑1 (MQ=255)                                                                         tGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa                                                         >  2:946428/1‑140 (MQ=255)                                                                           ggAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg                                                       <  2:180735/140‑1 (MQ=255)                                                                                 cGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt                                                 >  1:1498848/1‑140 (MQ=255)                                                                                     gctTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt                                                 >  2:974864/1‑136 (MQ=255)                                                                                     gctTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGt                                                 <  1:974864/136‑1 (MQ=255)                                                                                     gctTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACa                                             >  2:1329619/1‑140 (MQ=255)                                                                                     gctTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACa                                             <  1:415133/140‑1 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGAc                                              <  2:618712/128‑1 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGAc                                              >  1:618712/1‑128 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGaa                                    <  1:1052138/138‑1 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGaa                                    >  2:1052138/1‑138 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACt                                  >  2:1158888/1‑140 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACt                                  >  2:1431286/1‑140 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACt                                  >  2:1490970/1‑140 (MQ=255)                                                                                                cGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACt                                  >  2:28504/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  2:1116955/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  1:1597169/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  1:1597948/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  2:270946/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  1:637361/1‑140 (MQ=255)                                                                                                     tGCCGGTGATGCCACCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcgg                             >  1:555130/1‑140 (MQ=255)                                                                                                      gCCGGTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTggcggc                            >  2:1091266/1‑140 (MQ=255)                                                                                                          gTGATGCCGCCAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGTCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAAcc                        >  1:1030784/1‑140 (MQ=255)                                                                                                                  gccAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATggg                <  1:1329619/140‑1 (MQ=255)                                                                                                                  accAGACTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATggg                <  2:555130/139‑1 (MQ=255)                                                                                                                       aCTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTa           >  2:620508/1‑140 (MQ=255)                                                                                                                       aCTGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTa           >  1:1160157/1‑140 (MQ=255)                                                                                                                         tGCGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAAt         <  2:1597169/140‑1 (MQ=255)                                                                                                                           cGAGAAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGc       <  1:1116955/140‑1 (MQ=255)                                                                                                                              gaAGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGcacc    <  1:1193920/140‑1 (MQ=255)                                                                                                                                aGTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGcaccac  <  1:519720/140‑1 (MQ=255)                                                                                                                                 gTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATg                  >  1:1609301/1‑123 (MQ=255)                                                                                                                                 gTTCTACAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATg                  <  2:1609301/123‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      AAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCAC  >  NC_000913/2672336‑2672602

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