Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,594	G→T	6.1%	E237D (GAG→GAT)	phhA →	phenylalanine‑4‑hydroxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,102,594	0	G	T	6.1%	198.8 / 3.9	66	E237D (GAG→GAT)	phhA	phenylalanine‑4‑hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

AGACCGTCTACAGTCTGTCTGACGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAA  >  NC_002947/5102449‑5102741                                                                                                                                                   |                                                                                                                                                    agaCCGTCTACAGTCTGTCTGACGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGAc                                                                                                                                                  <  6:89326/149‑1 (MQ=255)              tctgtctgACGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCATGATGTCCTCGTGGGccagg                                                                                                                                     <  3:54415/149‑5 (MQ=255)              tctgtctgACGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATgg                                                                                                                                     >  8:121057/1‑149 (MQ=255)                      aCGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGAc                                                                                                                                                  <  1:243781/128‑1 (MQ=255)                      aCGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGAc                                                                                                                                                  >  2:243781/1‑128 (MQ=255)                           ccTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCAt                                                                                                                        >  5:90994/1‑149 (MQ=255)                             tGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGTAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATg                                                                                                                       >  6:50088/1‑148 (MQ=255)                             tGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCa                                                                                                                                 <  6:322881/138‑1 (MQ=255)                             tGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCa                                                                                                                                 >  5:322881/1‑138 (MQ=255)                              gAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCATACGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATg                                                                                                                       >  1:82945/1‑148 (MQ=255)                                   ccAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGTAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTgg                                                                                                                <  5:50088/149‑1 (MQ=255)                                   ccAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTgg                                                                                                                >  1:168620/1‑149 (MQ=255)                                     aGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGCCATTCTGCAACCGGTGTATTACGTACTGCCGAACATGAAGCGCCTGTTCGACATGGCACACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGt                                                                                                              >  2:101288/1‑149 (MQ=255)                                      ggCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTc                                                                                                             >  6:88716/1‑149 (MQ=255)                                      ggCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTc                                                                                                             >  6:114963/1‑149 (MQ=255)                                       gCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCt                                                                                                            <  2:233295/149‑1 (MQ=255)                                          ttCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCtgca                                                                                                         >  6:109505/1‑149 (MQ=255)                                            cGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCtgcat                                                                                                        <  5:326841/148‑1 (MQ=255)                                            cGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCtgca                                                                                                         <  4:14252/147‑1 (MQ=255)                                               cccGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCAc                                                                                                    >  4:243824/1‑149 (MQ=255)                                                 cGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATg                                                                                                                       <  6:157734/128‑1 (MQ=255)                                                 cGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATg                                                                                                                       >  5:157734/1‑128 (MQ=255)                                                     gaaGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGt                                                                                              <  1:240295/147‑1 (MQ=255)                                                        ggCCATGCGTACACCCTCCCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTc                                                                                           <  6:90994/149‑1 (MQ=255)                                                            aTGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCAcc                                                                                       >  7:210235/1‑149 (MQ=255)                                                               cGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCAccc                                                                                      >  2:110181/1‑147 (MQ=255)                                                                     cccTACCGCAT‑CGACCTTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGc                                                                              <  6:27561/149‑1 (MQ=255)                                                                      ccTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGACGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCt                                                                             >  4:193527/1‑149 (MQ=255)                                                                           cGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCtcaggc                                                                        <  2:82945/149‑6 (MQ=255)                                                                            gCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCt                                                                         >  2:28946/1‑147 (MQ=255)                                                                            gCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGa                                                                       <  1:28946/149‑1 (MQ=255)                                                                              aT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATaaa                                                                                                                             <  1:114308/93‑1 (MQ=255)                                                                              aT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATaaa                                                                                                                             >  2:114308/1‑93 (MQ=255)                                                                                  gACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGAc                                                                  <  3:243824/149‑1 (MQ=255)                                                                                      ttCTGCAACCGGTGTATTTCTTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGTCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCtt                                                               <  5:20051/148‑1 (MQ=255)                                                                                      ttCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCACACGAGGACATCATGGGCATGGTCCATAAACCCATGCAGCTGGGTCAGCATGCACCCAAGTTTCCACCAAAGCTCGCTGCCTGAGCGACCtt                                                               >  6:20051/1‑148 (MQ=255)                                                                                        cTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGc                                                             <  5:114963/148‑1 (MQ=255)                                                                                           cAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTCAGGCAGCGACCTTgggtg                                                          <  1:110181/148‑5 (MQ=255)                                                                                              ccGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaac                                                      <  8:210235/149‑1 (MQ=255)                                                                                              ccGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaa                                                       >  5:213078/1‑148 (MQ=255)                                                                                              ccGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaa                                                       <  6:213078/148‑1 (MQ=255)                                                                                                 gtgtATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCAGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaacaac                                                   >  2:212009/1‑149 (MQ=255)                                                                                                    tATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCg                                                <  5:109505/149‑1 (MQ=255)                                                                                                      tttCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGaa                                              >  2:249642/1‑149 (MQ=255)                                                                                                          gTACTGCCGAACATGAAGCGCCTGTTCGACCTGTCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTcc                                                                                          <  6:161492/101‑1 (MQ=255)                                                                                                          gTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTcc                                                                                          >  5:161492/1‑101 (MQ=255)                                                                                                               gCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGaaaa                                      >  1:223450/1‑148 (MQ=255)                                                                                                               gCCGAACATGAAGCGCCTATTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGaaaa                                      >  5:251402/1‑148 (MQ=255)                                                                                                                ccGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGc                                                                     >  2:307417/1‑116 (MQ=255)                                                                                                                ccGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGc                                                                     <  1:307417/116‑1 (MQ=255)                                                                                                                ccGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAaca                                    >  6:124049/1‑149 (MQ=255)                                                                                                                 cGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAacac                                   <  1:249642/149‑1 (MQ=255)                                                                                                                  gAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGc                                                             >  3:150840/1‑122 (MQ=255)                                                                                                                  gAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGc                                                             <  4:150840/122‑1 (MQ=255)                                                                                                                  gAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACAcc                                  <  3:193527/149‑1 (MQ=255)                                                                                                                      aTGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCAt                              <  2:168620/149‑1 (MQ=255)                                                                                                                              ccTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCt                      <  6:236493/149‑1 (MQ=255)                                                                                                                               cTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCtt                     >  2:264380/1‑149 (MQ=255)                                                                                                                                  ttCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGaa                  <  5:88716/149‑1 (MQ=255)                                                                                                                                     gACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGCACCGGAAAACCCCTCATGAATGCCTTGAAcc                >  4:43861/1‑148 (MQ=255)                                                                                                                                       ccTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAg             <  6:251402/149‑1 (MQ=255)                                                                                                                                        catgCCCATGATGTCCTCGTGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGc            >  2:95331/4‑149 (MQ=255)                                                                                                                                          tgCCCATGATGTCCTCGTGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaa                                                       <  4:186531/103‑1 (MQ=255)                                                                                                                                          tgCCCATGATGTCCTCGTGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATaa                                                       >  3:186531/2‑104 (MQ=255)                                                                                                                                             ccACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCAtt       >  7:250640/1‑149 (MQ=255)                                                                                                                                                  aGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGaa  <  8:250640/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGACCGTCTACAGTCTGTCTGACGAGCCTGAGCACCAGGCCTTCGACCCGATCGAGGCCATGCGTACACCCTACCGCAT‑CGACATTCTGCAACCGGTGTATTTCGTACTGCCGAACATGAAGCGCCTGTTCGACCTGGCCCACGAGGACATCATGGGCATGGTCCATAAAGCCATGCAGCTGGGTCTGCATGCACCGAAGTTTCCACCCAAGGTCGCTGCCTGAGCGACCTTGCCGATAACAACTCGAACCGGAAAACACCTCATGAATGCCTTGAACCAAGCCCATTGCGAA  >  NC_002947/5102449‑5102741

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