Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,448,947	A→C	100%	T477P (ACG→CCG)	feaB →	phenylacetaldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,448,947	0	A	C	100.0%	283.6 / NA	86	T477P (ACG→CCG)	feaB	phenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (41/45);  total (41/45)

ATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAC  >  NC_000913/1448853‑1449043                                                                                               |                                                                                                 atatAGCGATCGCTTACAGGCAGGGACGGTGTTGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGc                                                                                             <  2:1797225/100‑1 (MQ=255)   ataGCGATCGATTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCg                                                                                           <  2:935956/100‑1 (MQ=255)    taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGt                                                                                          >  2:307727/1‑100 (MQ=255)     aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTg                                                                                         <  1:598445/100‑1 (MQ=255)     aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTg                                                                                         <  2:105427/100‑1 (MQ=255)        gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAAGACCGGGCCGTGAtt                                                                                      >  1:2943584/1‑100 (MQ=255)           cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                   <  2:2748470/100‑1 (MQ=255)           cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                   >  2:2713974/1‑100 (MQ=255)           cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                   <  2:706752/100‑1 (MQ=255)            gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTgg                                                                                  >  2:2304780/1‑100 (MQ=255)             cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTgg                                                                                  <  2:2943552/99‑1 (MQ=255)             cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGc                                                                                 >  1:1239246/1‑100 (MQ=255)              ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTg                                                                                   <  1:1308005/97‑1 (MQ=255)              ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcc                                                                                <  1:1912806/100‑1 (MQ=255)              ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcc                                                                                <  2:1636361/100‑1 (MQ=255)               tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGccc                                                                               <  2:1131348/100‑1 (MQ=255)                 cttgcagaGACGGTGTGGTTAAACAGCAATACCTTAATTGACGCTAACTTAACGGTTAGGGCGATGAAGCAGTCAGGACCCGGCCGTGAattt                                                                                    >  2:1777597/4‑93 (MQ=255)                     caggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAAATTACAGTTTGGTGGAATGAAGCAGTCAGGACAGGGCAGTGATGTTGGc                                                                                 >  2:2252829/1‑92 (MQ=255)                      aggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActg                                                                         <  2:2081584/99‑1 (MQ=255)                        ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                      >  1:1107243/1‑100 (MQ=255)                         gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggctg                                                                     <  1:491370/100‑1 (MQ=255)                            ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                  >  1:1771851/1‑100 (MQ=255)                             gtgtGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                 >  2:1128962/1‑100 (MQ=255)                              tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg                                                                >  1:1178420/1‑100 (MQ=255)                              tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCAGTGATTTTGGACCCGACTGGCTGGACgg                                                                >  1:664983/1‑100 (MQ=255)                               gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                               <  2:1107234/100‑1 (MQ=255)                               gtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                               >  2:2582627/1‑100 (MQ=255)                               gtgGGTAAACAGCCATACCTTAATTGACGCTAACGTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                               <  2:2103236/100‑1 (MQ=255)                                tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt                                                              <  2:1052553/100‑1 (MQ=255)                                tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt                                                              <  1:2201608/100‑1 (MQ=255)                                 gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                             >  1:2403203/1‑100 (MQ=255)                                 gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                             >  2:1523966/1‑100 (MQ=255)                                 gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                             >  2:1204310/1‑100 (MQ=255)                                 gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                             >  2:467642/1‑100 (MQ=255)                                  ggTAAACAGCCATACCTTAATTGACGTTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                            <  1:2808440/100‑1 (MQ=255)                                  ggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                            <  1:3000904/100‑1 (MQ=255)                                  ggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                            <  2:2827472/100‑1 (MQ=255)                                  ggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                            <  1:2703211/100‑1 (MQ=255)                                    tAAACAGCCATACCTTAATTGACGCTAACTTCCCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                          >  1:1807990/1‑100 (MQ=255)                                    tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                          >  1:900849/1‑100 (MQ=255)                                    tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGAACCGGGCCATGATTTTGGCTCCGACTGGCTGGACGGTTGgtg                                                          >  1:219057/1‑100 (MQ=255)                                    tAAACAGCCATACCTTAATTAACGATAACATACCGATTGGTGAGATTAATCAGACAGGACATTGCCGTGATTATGGCACAGACTCGCTGCACGGt                                                               >  2:927104/1‑95 (MQ=255)                                     aaaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                         <  2:1178409/100‑1 (MQ=255)                                      aaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGtgtg                                                        >  1:2775201/1‑100 (MQ=255)                                        cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                      <  2:2940727/100‑1 (MQ=255)                                        cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGGGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                      <  2:2435215/100‑1 (MQ=255)                                          gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTAGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                    >  1:250275/1‑100 (MQ=255)                                            cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                  >  1:2749477/1‑100 (MQ=255)                                                ccTTAATTGACGCTAACTTACCGTTTGGTGGGATTAAGCAGTCAGGACCGGGCCATGAATTTGGTCCCGACTGGCTGGACGGTTGGTGTAAAACTAAGTc                                              >  1:1718486/1‑100 (MQ=255)                                                ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                              <  2:1930130/100‑1 (MQ=255)                                                ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                              <  1:307730/100‑1 (MQ=255)                                                ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGACCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                              >  2:1693773/1‑100 (MQ=255)                                                  ttAATTGACGCTAGCTTACCGTTTGGTGGTATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCgg                                            <  1:1761887/100‑1 (MQ=255)                                                  ttAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCgg                                            <  1:2521830/100‑1 (MQ=255)                                                    aaTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                          <  2:2200220/100‑1 (MQ=255)                                                        gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCTTGATTTTGGCCCCGACTGGCTGGACGAATGATATAAAACTAAGTCGgtgtg                                        >  1:332655/1‑98 (MQ=255)                                                          cGCTAACTTACCGTTTGGTAGGATGAAGCAGTCAGGACCGGGACGTGATTTTGGCCCCGACGGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                    >  2:522783/1‑100 (MQ=255)                                                            cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACg                                  >  1:2423414/1‑100 (MQ=255)                                                             tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGTACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                      <  2:2803571/95‑1 (MQ=255)                                                             tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACgg                                 <  2:2468285/100‑1 (MQ=255)                                                               aCTTACCGATTGGTGGGATGAAGCAGTCAGGACGGGGCAGCGATTTTGGCCACGACTGACTGGACGGTTGGTGTGACAGTAAGTCGGTGTGTGTCCg                                  >  2:2960368/1‑97 (MQ=255)                                                                cTTACCGTTTGGTGGGATGAAGCAGTCAGGACTGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                              <  1:213450/100‑1 (MQ=255)                                                                cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGATACTAAGTCGGTGTGTGTACGGTa                               <  2:728384/99‑1 (MQ=255)                                                                cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                              >  1:1600616/1‑100 (MQ=255)                                                                cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                              >  1:1338122/1‑100 (MQ=255)                                                                cTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTAGACGGTTGATGTGAAACTAAGTCGGTGTGTGTACGGTAt                              >  1:245704/1‑100 (MQ=255)                                                                  tACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                            <  2:543663/100‑1 (MQ=255)                                                                    ccGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAAt                          <  1:2155687/100‑1 (MQ=255)                                                                    ccGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGACTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAAt                          >  1:768355/1‑100 (MQ=255)                                                                     cGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATc                         <  2:1440406/100‑1 (MQ=255)                                                                      gTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCt                        <  1:2999865/100‑1 (MQ=255)                                                                       tttGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTg                       <  1:1012723/100‑1 (MQ=255)                                                                         tggtggGATGAAGCTGTCAGGACCGGGCCGTGATTTGGGTCCCGACTGGCTGTACGGTTGGCGTGAAACTAAGTCGGTGTTGGTATGGTATTAATCTGGt                     <  2:1018855/100‑1 (MQ=255)                                                                          ggtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                    <  2:2591262/100‑1 (MQ=255)                                                                                aTGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa              <  1:2213391/100‑1 (MQ=255)                                                                                  gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa            <  2:1454148/100‑1 (MQ=255)                                                                                  gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa            >  2:101708/1‑100 (MQ=255)                                                                                   aaGCAGTCAGGACCGGGCCGTGATTTTTGCTCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa           <  1:1452299/100‑1 (MQ=255)                                                                                   aaGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATaa           >  2:1290172/1‑100 (MQ=255)                                                                                    aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAg          >  1:1674118/1‑100 (MQ=255)                                                                                    aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAg          >  2:331013/1‑100 (MQ=255)                                                                                    aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGATCATAAg          <  2:249374/100‑1 (MQ=255)                                                                                     gCAGCCAGGACCAGGCCGTGACTTTGTCCCCGACTGGCTGGACGGTTGATGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCACTCATAAg          >  2:1089114/1‑99 (MQ=255)                                                                                       aGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGATGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaa       <  1:417/100‑1 (MQ=255)                                                                                        ntCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa      >  2:2337493/2‑100 (MQ=255)                                                                                            ggACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc  >  2:2838756/1‑100 (MQ=255)                                                                                               |                                                                                                 ATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAC  >  NC_000913/1448853‑1449043

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