Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	657,688	A→G	7.0%	H10R (CAT→CGT) ‡	PP_0565 →	Nudix hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	657,688	0	A	G	7.0%	205.4 / 5.0	71	H10R (CAT→CGT) ‡	PP_0565	Nudix hydrolase
Reads supporting (aligned to +/- strand):  ref base A (29/37);  new base G (2/3);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGCTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATTGTGT  >  NC_002947/657542‑657836                                                                                                                                                   |                                                                                                                                                     tGCTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATc                                                                                                                                                    <  4:185247/149‑1 (MQ=255) tGCTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATc                                                                                                                                                    >  3:60131/1‑149 (MQ=255)  gCTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGTTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCg                                                                                                                                                   >  3:182113/1‑149 (MQ=255)   cTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCgg                                                                                                                                                  <  3:177124/149‑1 (MQ=255)    tCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCggg                                                                                                                                                 >  6:186858/1‑149 (MQ=255)           cAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTc                                                                                                                                          <  1:114833/149‑1 (MQ=255)            aGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTcc                                                                                                                                         <  7:33517/149‑1 (MQ=255)            aGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTcc                                                                                                                                         <  7:58874/149‑1 (MQ=255)                gCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGGTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGAc                                                                                                                                      <  4:182113/148‑1 (MQ=255)                 cTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACcg                                                                                                                                    >  6:724/1‑149 (MQ=255)                 cTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAAGCCGCCCGATGGGCGGCCTcggttgg                                                                                                                                    <  1:4381/149‑7 (MQ=255)                  ttCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTcc                                                                                                                                         >  1:220073/1‑143 (MQ=255)                  ttCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTcc                                                                                                                                         <  2:220073/143‑1 (MQ=255)                         tgtgGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTgg                                                                                                                            <  5:330833/149‑1 (MQ=255)                             ggCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCa                                                                                                                         >  6:282146/1‑148 (MQ=255)                             ggCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCAc                                                                                                                        >  3:49466/1‑149 (MQ=255)                               ccAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGATGCCGCCCATCGGGCGGCTTCCGACCCCGAGCTGGTCACCt                                                                                                                      <  4:128749/149‑1 (MQ=255)                               ccAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCt                                                                                                                      >  3:79651/1‑149 (MQ=255)                                 aaCCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTgg                                                                                                                    <  8:166914/149‑1 (MQ=255)                                 aaCCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTgg                                                                                                                    >  3:331409/1‑149 (MQ=255)                                  aCCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTggg                                                                                                                   <  4:331409/149‑1 (MQ=255)                                   ccACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGTGGCAGCTGCTGTCGCAGGGTTTATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGt                                                                                                                  <  6:298204/149‑1 (MQ=255)                                   ccACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGTCTTCCGACCGCGAGCTGGTCACCTGGGt                                                                                                                  <  5:724/149‑1 (MQ=255)                                   caaCTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGTGGCTTCCGACCGCGAGCTGGTCACCTGGGt                                                                                                                  >  1:361120/3‑149 (MQ=255)                                    cACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCAcc                                                                                                                       <  7:48705/143‑1 (MQ=255)                                    cACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCAcc                                                                                                                       >  8:48705/1‑143 (MQ=255)                                         gTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGAc                                                                                                              <  1:265848/147‑1 (MQ=255)                                          ttGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGTGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGaa                                                                                                           <  2:361120/149‑1 (MQ=255)                                                aTGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGAc                                                                                                     >  8:18679/1‑149 (MQ=255)                                                 tGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGAcc                                                                                                    <  4:380599/149‑1 (MQ=255)                                                   gCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAg                                                                                                  <  1:62144/149‑1 (MQ=255)                                                   gCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAg                                                                                                  <  8:195064/149‑1 (MQ=255)                                                    ccAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAgg                                                                                                 <  4:60131/149‑1 (MQ=255)                                                      aaCTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTa                                                                                               >  7:17134/1‑149 (MQ=255)                                                        cTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACt                                                                                             <  8:272494/149‑1 (MQ=255)                                                         tACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTg                                                                                            <  5:60442/149‑1 (MQ=255)                                                           cATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTggg                                                                                          <  8:226980/149‑1 (MQ=255)                                                            aTTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGt                                                                                                <  6:347965/142‑1 (MQ=255)                                                            aTTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGt                                                                                                >  5:347965/1‑142 (MQ=255)                                                            aTTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGgc                                                                                         <  5:186858/149‑1 (MQ=255)                                                            aTTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGGGGACGAAACCGACCAGGTACTGGgc                                                                                         >  4:101389/1‑149 (MQ=255)                                                             tttCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAAGCCGCCCGATGGGCGGCCTc                                                                                                                                          <  1:145005/99‑1 (MQ=255)                                                             tttCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAAGCCGCCCGATGGGCGGCCTc                                                                                                                                          >  2:145005/1‑99 (MQ=255)                                                               tCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGgcgc                                                                                       <  3:101389/148‑1 (MQ=255)                                                                               cGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAg                                                                       <  4:49466/148‑1 (MQ=255)                                                                                  ttATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTgc                                                                   <  8:17134/149‑1 (MQ=255)                                                                                     tCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTgcgc                                                                 >  2:205949/1‑148 (MQ=255)                                                                                         ccTGCGACAGCAGCTGACGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCg                                                            >  5:260702/2‑149 (MQ=255)                                                                                             cGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTg                                                          <  7:18679/147‑1 (MQ=255)                                                                                                 aGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGa                                                    >  5:294476/1‑149 (MQ=255)                                                                                                    agctgACGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCg                                                 <  6:260702/148‑1 (MQ=255)                                                                                                      ctgTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCACATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACGAGGTACTGGTCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGt                                               <  6:375826/149‑1 (MQ=255)                                                                                                         tCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGt                                            >  1:147752/1‑149 (MQ=255)                                                                                                         tCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGt                                            >  3:215957/1‑149 (MQ=255)                                                                                                                 ttGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGtt                                    <  5:282146/149‑1 (MQ=255)                                                                                                                 ttGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCCGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGtt                                    >  1:154027/1‑149 (MQ=255)                                                                                                                    atgatgGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGtgtg                                                          >  6:215910/1‑124 (MQ=255)                                                                                                                    atgatgGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGtgtg                                                          <  5:215910/124‑1 (MQ=255)                                                                                                                    atgatgGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCa                                  <  1:203872/148‑1 (MQ=255)                                                                                                                             aTCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGtt                        >  6:394628/1‑149 (MQ=255)                                                                                                                                aGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCaa                     <  1:205949/149‑1 (MQ=255)                                                                                                                                 gCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCaa                     >  5:14253/1‑148 (MQ=255)                                                                                                                                 gCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCaa                     >  1:20716/1‑148 (MQ=255)                                                                                                                                 gCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCaa                     >  8:226076/1‑148 (MQ=255)                                                                                                                                     aaCCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTg                >  1:348077/1‑149 (MQ=255)                                                                                                                                          aagccgccCGATGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCgg                                                <  8:42926/110‑1 (MQ=255)                                                                                                                                          aagccgccCGATGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCgg                                                >  7:42926/3‑112 (MQ=255)                                                                                                                                          aGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGt           >  5:202285/1‑149 (MQ=255)                                                                                                                                          aGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGt           <  7:16171/149‑1 (MQ=255)                                                                                                                                                 ccATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATtgt    <  2:348077/149‑1 (MQ=255)                                                                                                                                                  cATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATtgtg   >  7:141112/1‑149 (MQ=255)                                                                                                                                                   aTCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATtgtgt  <  7:226076/149‑1 (MQ=255)                                                                                                                                                   aTCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATtgtgt  <  5:394628/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGCTCAACGGCAGTGGCTTCGGCATGTGGGCCAACCACTCGTTGGTCATGGCCAACTACATTTCCCGGGCCAACGCTGCGGTTATCGATCTGCGTCAGCTGCTGTCGCAGGGTTGATGATGGCTATCAGCCCAACCGAGGCCGCCCATCGGGCGGCTTCCGACCGCGAGCTGGTCACCTGGGTGGACGAAACCGACCAGGTACTGGGCGCCTTGCCCCGTGCCGAGCTGCGCGAGCGTGGCCTGATCGGTCGTTGCACGTTCATTCTGCTGTTCAACAGTGCCGGTGAATTGTGT  >  NC_002947/657542‑657836

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: ‑