Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,802	A→C	5.9%	K208Q (AAG→CAG) ‡	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,802	0	A	C	5.9%	214.2 / 4.5	68	K208Q (AAG→CAG) ‡	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base C (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101657‑101950                                                                                                                                                   |                                                                                                                                                     gagCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   <  5:127793/149‑1 (MQ=255)    cAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGcc                                                                                                                                                >  5:291980/1‑149 (MQ=255)      aGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTg                                                                                                                                              <  6:70284/149‑1 (MQ=255)       gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   >  6:276727/1‑143 (MQ=255)       gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAgg                                                                                                                                                   <  5:276727/143‑1 (MQ=255)       gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTg                                                                                                                                              >  6:239954/1‑148 (MQ=255)       gCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTg                                                                                                                                              >  5:362825/1‑148 (MQ=255)         tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           >  7:233151/1‑149 (MQ=255)         tCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           >  7:134841/1‑149 (MQ=255)          cAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                               >  5:125569/1‑144 (MQ=255)          cAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCt                                                                                                                                               <  6:125569/144‑1 (MQ=255)             cGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCt                                                                                                                                       <  6:291980/149‑1 (MQ=255)                    ggCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCATGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGa                                                                                                                                 <  5:239954/148‑1 (MQ=255)                       aGCCGTTTCTCCAATCAGCTGCTAGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAc                                                                                                                             >  7:51818/1‑149 (MQ=255)                               ctcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGaa                                                                                                                      <  3:9014/148‑1 (MQ=255)                               ctcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGaa                                                                                                                      <  8:134841/148‑1 (MQ=255)                                tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTAGCCGACGAAGATCCGCTGGCCGGCCTGACCGATACGGCCAAGCCGCAGATGGCCGATGCAGCCAAGGCCAAGGGCCTGGACGGATGGTTGATCACCCTTGaa                                                                                                                      >  1:252896/1‑148 (MQ=255)                                tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAg                                                                                                                                                    >  8:264317/1‑117 (MQ=255)                                tcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAg                                                                                                                                                    <  7:264317/117‑1 (MQ=255)                                  cattctgctgctGGACGCCACCCTGGCCTGGACCATGCATGTTA‑CCGAGTACGGTGCGGTGGGCGGCCTGACCTTTTCGGCCATGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTc                                                                                                                   <  8:255107/142‑1 (MQ=255)                                  cAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTcc                                                                                                                  <  8:233151/149‑1 (MQ=255)                                   aaTCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTgggc                                                                                                                                           <  2:298152/123‑3 (MQ=255)                                   aaTCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTgggc                                                                                                                                           >  1:298152/1‑121 (MQ=255)                                   aaTCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTcc                                                                                                                  >  7:260491/1‑148 (MQ=255)                                     tccgctgctGGACGCCACCCAGGCCTGCACAAACCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGTTTCGGCCATGGCGCATATGGCCGCTTCAGCCTTGTCTATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                               <  8:51818/146‑1 (MQ=255)                                        ggtgctgGAGGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCTGGCCAAGGGCCTGGACGGCTGGTTg                                                                                                                                  <  2:124247/125‑1 (MQ=255)                                               gtcGCCACCCAGGCCGGGACCAAGCATGTTA‑CCGACGAAGCTGCGGTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgc                                                                                                     <  6:362825/147‑1 (MQ=255)                                                 cGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCgg                                                                                                   <  7:10878/149‑1 (MQ=255)                                                  gCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                  <  7:276145/149‑1 (MQ=255)                                                    cACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGa                                                                                                                       <  5:239392/126‑1 (MQ=255)                                                    cACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGa                                                                                                                       >  6:239392/1‑126 (MQ=255)                                                       ccAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           >  3:79349/1‑103 (MQ=255)                                                       ccAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                           <  4:79349/103‑1 (MQ=255)                                                       ccAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCt                                                                                                                          >  3:285052/1‑120 (MQ=255)                                                       ccAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCt                                                                                                                          <  4:285052/120‑1 (MQ=255)                                                           gCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 <  5:169042/125‑1 (MQ=255)                                                           gCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTccc                                                                                                                 >  6:169042/1‑125 (MQ=255)                                                           gCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGAGGACCTa                                                                                         >  1:293360/1‑149 (MQ=255)                                                                 aCCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                            <  2:203498/92‑2 (MQ=255)                                                                    aaGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGa                                                                                >  3:282033/1‑149 (MQ=255)                                                                    aaGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGa                                                                                >  3:172487/1‑149 (MQ=255)                                                                    aaGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGa                                                                                >  1:269683/1‑149 (MQ=255)                                                                        aTGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACcgc                                                                            <  2:293360/149‑1 (MQ=255)                                                                           ttA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                    >  8:188766/1‑122 (MQ=255)                                                                                cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgc                                                                     >  2:214016/1‑149 (MQ=255)                                                                                cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgc                                                                     <  6:267558/149‑1 (MQ=255)                                                                                cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgc                                                                     >  1:327413/1‑149 (MQ=255)                                                                                 gacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcg                                                                    <  4:172487/149‑1 (MQ=255)                                                                                 gacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcg                                                                    >  4:269449/1‑149 (MQ=255)                                                                                      aGCTGCGCTGTCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCTCGCCCTGCGCGAag                                                               <  4:282033/149‑1 (MQ=255)                                                                                           cgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                          <  2:252896/149‑1 (MQ=255)                                                                                            gcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctact                                                                                                         <  6:364309/101‑1 (MQ=255)                                                                                            gcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctact                                                                                                         >  5:364309/1‑101 (MQ=255)                                                                                                       tGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTAc                                              >  5:37108/1‑149 (MQ=255)                                                                                                         aCCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGTCTGGTTGATCAGCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGGGCCGTGCGCGAAGAGCTATACGCCGCCTACt                                             <  8:260491/148‑1 (MQ=255)                                                                                                          ccGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTg                                            >  8:240396/1‑148 (MQ=255)                                                                                                            gATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAc                                         <  2:327413/149‑1 (MQ=255)                                                                                                            gATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAc                                         >  4:111145/1‑149 (MQ=255)                                                                                                                             gATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAgg                        <  2:269683/149‑1 (MQ=255)                                                                                                                                   cGCTGCCCCCCCGCCCATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                          <  2:166017/109‑1 (MQ=255)                                                                                                                                   cGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                          >  1:166017/1‑109 (MQ=255)                                                                                                                                        ccgtccACGCCCATGGCCTGGACCGCTGGTTGATCACCCTTGAATTAACCAGCTACTACCCGGTGATGACCTACCCCAGCGACCGCGCCCTGCGCGAACAGCAATACCCCGCCTACTGCACCCGCGCCTCCGACCACCGCCCGAATGcc             >  2:146608/5‑149 (MQ=255)                                                                                                                                           cccAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc          <  3:111145/149‑1 (MQ=255)                                                                                                                                           cccAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc          <  3:269449/149‑1 (MQ=255)                                                                                                                                           cccAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc          <  7:240396/149‑1 (MQ=255)                                                                                                                                             cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                    <  3:343004/57‑1 (MQ=255)                                                                                                                                             cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTAcgcg                                                                                                    >  4:343004/1‑57 (MQ=255)                                                                                                                                                  cAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc    >  3:179353/1‑148 (MQ=255)                                                                                                                                                   aaGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGa  >  4:59660/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GAGCAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTA‑CCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGA  >  NC_002947/101657‑101950

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