Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	20,298	T→C	6.0%	L34P (CTG→CCG)	PP_0016 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	20,298	0	T	C	6.0%	196.7 / 2.1	68	L34P (CTG→CCG)	PP_0016	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CTACTCCCCTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACTCATCACG  >  NC_002947/20153‑20441                                                                                                                                                  |                                                                                                                                                ctACTCCCCTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCt                                                                                                                                              >  1:284424/1‑149 (MQ=255)    cTCCCCTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACCGCTCAc                                                                                                                                           >  2:180584/1‑149 (MQ=255)    cTCCCCTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAAACCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACCGCTCAc                                                                                                                                           >  6:197185/1‑149 (MQ=255)        ccTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGcc                                                                                                                                       >  2:272792/1‑149 (MQ=255)        ccTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACAACACCAAGACGAAAGTCTTTTGGGGTAACTGCAGCGCATGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGcc                                                                                                                                       >  5:65003/1‑149 (MQ=255)            aTTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCttt                                                                                                                                    <  1:166777/148‑1 (MQ=255)                ttATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCAc                                                                                                                                <  6:247057/148‑1 (MQ=255)                  aTGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGa                                                                                                                              <  4:31560/148‑1 (MQ=255)                  aTGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAg                                                                                                                             >  4:193793/1‑149 (MQ=255)                     tCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACATGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGcc                                                                                                                           >  7:107192/1‑148 (MQ=255)                       ttCGTTTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACGCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCcgac                                                                                                                        <  2:17367/149‑1 (MQ=255)                            tttCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGAc                                                                                                                   <  3:14285/149‑1 (MQ=255)                                        gtAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGcc                                                                                                                           >  8:217491/1‑129 (MQ=255)                                        gtAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGcc                                                                                                                           <  7:217491/129‑1 (MQ=255)                                         tAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGAc                                                                                                      <  6:115004/149‑1 (MQ=255)                                         tAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGAc                                                                                                      <  1:322502/149‑1 (MQ=255)                                         tAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGAc                                                                                                      >  4:57774/1‑149 (MQ=255)                                         tAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGAc                                                                                                      >  5:117588/1‑149 (MQ=255)                                           ggAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGAc                                                                                                      >  2:322502/1‑147 (MQ=255)                                            gAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTgg                                                                                                   >  5:127713/1‑149 (MQ=255)                                               tGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTg                                                                                                >  3:92399/1‑149 (MQ=255)                                                    ggTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAga                                                                                            <  3:57774/148‑1 (MQ=255)                                                         tttAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGatgatg                                                                                       >  7:290764/1‑148 (MQ=255)                                                           tAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGa                                                                                     <  3:193793/148‑1 (MQ=255)                                                           tAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGa                                                                                     <  2:326614/148‑1 (MQ=255)                                                            aaTCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGaa                                                                                    <  6:316041/148‑1 (MQ=255)                                                               cAACCCACGCCCCCACCAAGCCGCAAGTCTTTTGGGGTACCTGCCGCGCCTGGCGACAGCCAACGGATTACCCAGAACGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGCCCGATGTTGCGGACTGGCTGCAGATGATGCAAGGAc                                                                                <  4:92399/149‑1 (MQ=255)                                                                  cccACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCAc                                                                             >  5:197841/1‑149 (MQ=255)                                                                           acacCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTTAGCTGGc                                                                    <  8:290764/149‑1 (MQ=255)                                                                           acacCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGc                                                                    <  6:197841/149‑1 (MQ=255)                                                                               cAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGc                                                                >  3:204932/1‑149 (MQ=255)                                                                               cAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGc                                                                >  3:319183/1‑149 (MQ=255)                                                                                agggcGAAAGTATTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTTCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCa                                                                              <  8:26196/130‑1 (MQ=255)                                                                                accaCGCAACTCATTTGCGGAACCTCCAGCGCCTGGCGAACGCAAACGGCTTCCCCCCAAAGGAACAGCTCCCCGCCTTTCCACGACCCGACGAGACCGATGTTGCCGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCg                                                               <  8:107192/146‑1 (MQ=255)                                                                                aaGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCg                                                               <  6:127713/149‑1 (MQ=255)                                                                                 aGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACCGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCgg                                                              <  1:180584/149‑1 (MQ=255)                                                                                 aGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACCGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCgg                                                              <  5:197185/149‑1 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATg                                                                                                              >  2:120884/1‑97 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATg                                                                                                              <  1:120884/97‑1 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGTGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATc                                                          >  1:127983/1‑149 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAga                                                                                            <  3:315450/115‑1 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAga                                                                                            >  4:315450/1‑115 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATc                                                          >  1:293745/1‑149 (MQ=255)                                                                                     gAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATc                                                          >  6:193370/1‑149 (MQ=255)                                                                                      aaaGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAg                                                                                   <  2:37258/123‑1 (MQ=255)                                                                                      aaaGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAg                                                                                   >  1:37258/1‑123 (MQ=255)                                                                                        aGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGTCTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATctgc                                                       <  1:244144/149‑1 (MQ=255)                                                                                         gTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATctgct                                                      <  6:117588/149‑1 (MQ=255)                                                                                              ttGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGa                                                    <  1:272792/146‑1 (MQ=255)                                                                                              ttGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGa                                                    <  3:3120/146‑1 (MQ=255)                                                                                              ttGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGtt                                                 >  4:303633/1‑149 (MQ=255)                                                                                                 gggTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGTGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCg                                              <  2:127983/149‑1 (MQ=255)                                                                                                 gggTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTc                                               <  4:319183/148‑1 (MQ=255)                                                                                                  ggTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGt                                             <  2:293745/149‑1 (MQ=255)                                                                                                  ggTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGt                                             <  6:65003/149‑1 (MQ=255)                                                                                                    tACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGa                                                    >  5:252709/1‑140 (MQ=255)                                                                                                    tACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGa                                                    <  6:252709/140‑1 (MQ=255)                                                                                                         gCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATcc                                        >  4:3361/1‑147 (MQ=255)                                                                                                                            aaCGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGtt                   <  3:303633/149‑1 (MQ=255)                                                                                                                                   tACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGt                                             >  6:285276/1‑116 (MQ=255)                                                                                                                                   tACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGt                                             <  5:285276/116‑1 (MQ=255)                                                                                                                                     cccAGAAAGGAACTGCTCACCGCCTTTGTACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCg                                              >  6:39702/1‑113 (MQ=255)                                                                                                                                      ccAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACt         >  5:307708/1‑149 (MQ=255)                                                                                                                                      ccAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACt         >  5:114965/1‑149 (MQ=255)                                                                                                                                      ccAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATATGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACt         >  6:211558/1‑149 (MQ=255)                                                                                                                                       cAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACtc        >  6:236835/1‑149 (MQ=255)                                                                                                                                        aGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACtca       >  3:139264/1‑149 (MQ=255)                                                                                                                                             ggAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACTCATCACg  >  7:73371/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                CTACTCCCCTGATTTTTATGTCTTCGATTTCATGGATGTGTAGGAATGTTTGGTCTTTTAATCAACCCACGACCACACCAAGACGAAAGTCTTTTGGGGTACCTGCAGCGCCTGGCGAAAGCCAACGGATTACCCAGAAAGGAACTGCTCACCGCCTTTGCACGAGCCGACGAGACCGATGTTGCGGACTGGCTGCAGATGATGGAAGGACCACTGAGCTGGCCTGCGGTATCTGCTGAGTTTCGTGATCCAAGGCCCAAGGGGGTCAAGTTATGGACGACTCATCACG  >  NC_002947/20153‑20441

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