Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,258,559	Δ1 bp	8.3%	coding (180/1005 nt)	asd →	aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,258,559	0	C	.	8.3%	186.1 / 11.6	60	coding (180/1005 nt)	asd	aspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/22);  new base . (3/2);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAG  >  NC_002947/2258414‑2258707                                                                                                                                                  |                                                                                                                                                       ccACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCt                                                                                                                                                    >  1:1363556/1‑150 (MQ=255) ccACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCt                                                                                                                                                    >  2:298545/1‑150 (MQ=255)  cACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCt‑t                                                                                                                                                  <  1:245753/150‑1 (MQ=255)   aCCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑Tc                                                                                                                                                 <  1:1884596/150‑1 (MQ=255)       gCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAATTCGACAGCT‑TC‑GACt                                                                                                                                            >  1:856456/1‑150 (MQ=255)           cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGA‑AGCTGTC‑GACTTCGc                                                                                                                                        <  2:29932/150‑1 (MQ=255)               ggTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAg                                                                                                                                    <  2:1640408/150‑1 (MQ=255)                      ccTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGc                                                                                                                             >  2:232470/1‑150 (MQ=255)                           gTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGcc                                                                                                                        <  1:1712436/150‑1 (MQ=255)                                gATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCcttctt                                                                                                                   <  2:856456/150‑1 (MQ=255)                                 aTCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCttcttc                                                                                                                  >  2:172100/1‑150 (MQ=255)                                      cgaggcaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTACGCGAAGTCGACAGCT‑TC‑GCCTTCGCCCAGGTCCAGCTGGCCTTCTTCGCCAc                                                                                                             >  2:1765673/7‑150 (MQ=255)                                       gaagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTg                                                                                                                           >  1:1392042/1‑135 (MQ=255)                                       gaagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTg                                                                                                                           <  2:1392042/135‑1 (MQ=255)                                        aagaaCTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCg                                                                                                           >  1:1863188/1‑150 (MQ=255)                                             cTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGcgcc                                                                                                      >  2:523974/1‑150 (MQ=255)                                              tGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAc                                                                                                                                                        <  2:1257156/101‑1 (MQ=255)                                              tGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAc                                                                                                                                                        >  1:1257156/1‑101 (MQ=255)                                                  aTTGCCGGTTGCCACTTTGCACCTGCTGGCCAGCATGGAATGGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAGACTGAGGGGGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCGTCTTCGCCACCGGCGCCGCCGt                                                                                                 <  2:1745972/150‑1 (MQ=255)                                                   ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GAc                                                                                                                                             >  2:399396/1‑105 (MQ=255)                                                   ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GAc                                                                                                                                             <  1:399396/105‑1 (MQ=255)                                                      ccGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGc                                                                                             <  2:1015701/150‑1 (MQ=255)                                                              cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAc                                                                                                                                                        <  2:226331/85‑1 (MQ=255)                                                              cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGAc                                                                                                                                                        >  1:226331/1‑85 (MQ=255)                                                              cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGAACTGTTAGGACGGATAAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCtt                                                                                     >  2:1911045/1‑150 (MQ=255)                                                               aCTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                    <  1:298545/150‑1 (MQ=255)                                                                   tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                <  1:1294044/150‑1 (MQ=255)                                                                        ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                <  1:2591/145‑1 (MQ=255)                                                                        ctgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                                >  2:2591/1‑145 (MQ=255)                                                                             ggCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGGCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCt                                                                      >  1:385809/1‑150 (MQ=255)                                                                               ccAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGc                                                                    >  1:858914/1‑150 (MQ=255)                                                                                  gCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGTCAAGGCGCTGCAgg                                                                 >  1:1626398/1‑150 (MQ=255)                                                                                    ctggatccggcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GCCTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCg                                                               <  1:1461657/143‑1 (MQ=255)                                                                                     tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                              <  1:1126557/150‑1 (MQ=255)                                                                                     tGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCACCGCCGTCAGACGCTGCTGTGTCGGCGAGGCGCTGCAGGCgg                                                              >  2:575742/1‑150 (MQ=255)                                                                                          tCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCGAGGGCAAGCTTGCCGTCTTCGCCACCGGCGCCTCCGGCATCCGCAGCTTTGCCGTCAATGCGCTGCAGGCGGGCTgg                                                         >  2:1841630/1‑149 (MQ=255)                                                                                             gcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACg                                                      <  2:1363556/150‑1 (MQ=255)                                                                                               gggCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGt                                                    <  1:523974/150‑1 (MQ=255)                                                                                                   agcagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATc                                                <  2:1863188/150‑1 (MQ=255)                                                                                                      agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                             >  2:1130037/1‑150 (MQ=255)                                                                                                       gcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAACTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAcc                                            >  1:1724947/1‑150 (MQ=255)                                                                                                           gATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTc                                        >  2:755655/1‑150 (MQ=255)                                                                                                               ttCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCAACGTCACCTCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCcggc                                    >  1:1554002/1‑150 (MQ=255)                                                                                                                    tGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCt                               <  2:858914/150‑1 (MQ=255)                                                                                                                    tGGCAAGAAACAGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGATCATCGACCTGTCCGGCGGccc                               >  2:715876/1‑149 (MQ=255)                                                                                                                      gCAAGAAACTGAAGGTGCGCGAAGTCGACAGAT‑TCAGGATTATGCTAGG‑CAACATAAGTTCCTTCTCCACCGACACCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTgg                             >  1:514203/1‑150 (MQ=255)                                                                                                                        aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgac                           >  1:1645874/1‑150 (MQ=255)                                                                                                                        aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgac                           >  1:1520562/1‑150 (MQ=255)                                                                                                                          gAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGgacga                         >  2:359512/1‑150 (MQ=255)                                                                                                                              cTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgcc                     >  2:185586/1‑150 (MQ=255)                                                                                                                               tGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccc                    >  2:1804638/1‑150 (MQ=255)                                                                                                                                     tgggcgAAGTCGA‑AGCTGTC‑GACTTCGCGCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGGCATCGACCGTTCCGGCGGCCTGGACGACgccctggcac              >  2:1260300/4‑148 (MQ=255)                                                                                                                                     tGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccc              >  2:806264/1‑150 (MQ=255)                                                                                                                                       cgcgAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg            <  1:172100/150‑1 (MQ=255)                                                                                                                                        gcgAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGccctggccctgg           >  1:159997/1‑150 (MQ=255)                                                                                                                                          gAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTg         >  2:614233/1‑150 (MQ=255)                                                                                                                                          gAAGTCGA‑AGCTGTC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                   <  2:953364/140‑1 (MQ=255)                                                                                                                                          gAAGTCGA‑AGCTGTC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccct                   >  1:953364/1‑140 (MQ=255)                                                                                                                                          gAAGTCGA‑AGCTGTC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTg         >  1:1759830/1‑150 (MQ=255)                                                                                                                                             gTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGATGTACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCg      <  1:1573590/150‑1 (MQ=255)                                                                                                                                                gACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa   <  1:614233/150‑1 (MQ=255)                                                                                                                                                gACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa   <  1:1911045/150‑1 (MQ=255)                                                                                                                                                gACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa   >  1:1687627/1‑150 (MQ=255)                                                                                                                                                 aCAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAg  >  1:5742/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       CCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TC‑GACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAAG  >  NC_002947/2258414‑2258707

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