Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	436,899	C→A	8.0%	W130L (TGG→TTG)	PP_0358 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	436,899	0	C	A	8.0%	185.2 / 4.4	59	W130L (TGG→TTG)	PP_0358	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base A (3/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAT  >  NC_002947/436751‑437041                                                                                                                                                     |                                                                                                                                                caGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGG‑TCAGCCCATAGCAGCCGAAGGTCAACGCCAGCACCAGTGATACCCCCGGCAGGCTGCCCCCTTGCCCCACCTGCTGCGCCACCCCCACGGCAGCCATGCCCACCgccagcca                                                                                                                                                >  7:114595/1‑149 (MQ=255)  aGGCCGGGCAGCGCAGCCACCGTTGCCTTATTACGGATCTGCCCTTTGATGCCGATGCTCAACGCCAGCACCAGTGATACGCACGTCAGGCTGCCCACTTGCCACACCTGCTGCTCCACACCCTCGGCAGCCATGCCCACCgcctgcca                                                                                                                                                <  2:320414/149‑1 (MQ=255)  aGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCgccagccg                                                                                                                                                <  4:231323/149‑2 (MQ=255)     ccGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCGACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTg                                                                                                                                             >  6:175407/1‑149 (MQ=255)       gggCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCa                                                                                                                                           <  5:351590/149‑1 (MQ=255)         gcagcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTgg                                                                                                                                         <  3:99713/149‑1 (MQ=255)                agcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCa                                                                                                                                  <  5:223681/149‑1 (MQ=255)                     ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCCACCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGACGgcgc                                                                                                                             >  4:67071/1‑149 (MQ=255)                     ccGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGgcgc                                                                                                                             <  5:15848/149‑1 (MQ=255)                      cGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCt                                                                                                                            <  1:193075/149‑1 (MQ=255)                            cTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCaa                                                                                                                      <  8:210787/149‑1 (MQ=255)                               cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGaa                                                                                                                   >  7:262681/1‑149 (MQ=255)                                   cGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                               <  6:337811/149‑1 (MQ=255)                                      aTCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCAt                                                                                                            <  1:303214/149‑1 (MQ=255)                                      aTCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCAt                                                                                                            <  3:187535/149‑1 (MQ=255)                                             cATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCGCACCAAGAATCAGCATGC‑Ccag                                                                                                     >  1:208985/1‑148 (MQ=255)                                                       gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGa                                                                                          >  6:84271/1‑149 (MQ=255)                                                       gAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGa                                                                                          >  1:187231/1‑149 (MQ=255)                                                         aGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAgg                                                                                                                                         <  4:254133/101‑1 (MQ=255)                                                         aGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAgg                                                                                                                                         >  3:254133/1‑101 (MQ=255)                                                            tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGc                                                                                     >  6:150392/1‑149 (MQ=255)                                                               aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCggg                                                                                  <  5:175407/149‑1 (MQ=255)                                                                    aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                             >  8:277708/1‑149 (MQ=255)                                                                    aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                             <  8:262681/149‑1 (MQ=255)                                                                    aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                             <  7:277708/149‑1 (MQ=255)                                                                     gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACTTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATg                                                                            <  1:301094/149‑1 (MQ=255)                                                                     gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATg                                                                            <  2:220372/149‑1 (MQ=255)                                                                      cACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCATACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGc                                                                                                              <  2:123993/115‑1 (MQ=255)                                                                       aCCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCAACTGCAGCCATGCCCACAGCCAGCAATTGCAGGCGGCGCAGGCGCTCACGAAGAAACAGCATGC‑CCAGCAGTAAATTGATCAGCGGGTTGATGt                                                                           >  6:54021/1‑148 (MQ=255)                                                                          aGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATa                                                                       <  5:150392/149‑1 (MQ=255)                                                                              aTACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGcc                                                                    >  5:217806/1‑148 (MQ=255)                                                                                  ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGc                                                               >  1:359390/1‑149 (MQ=255)                                                                                  ccACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCACGCCCACCGCCTGCAATGGCTggcggtgg                                                                                                                                   >  2:105661/1‑79 (MQ=37)                                                                                   cACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                >  7:317609/1‑147 (MQ=255)                                                                                   cACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                <  8:317609/147‑1 (MQ=255)                                                                                         aGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTc                                                        <  8:278110/149‑1 (MQ=255)                                                                                             tGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAg                                                     >  4:4065/1‑148 (MQ=255)                                                                                               cccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAAGTAATAGCCCAGGCTGGCCTCGAGc                                                    <  2:208985/147‑1 (MQ=255)                                                                                                 cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                             >  2:356973/1‑120 (MQ=255)                                                                                                 cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGAt                                                                             <  1:356973/120‑1 (MQ=255)                                                                                                    ttGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGATGGCCTCGAGCATACGGc                                             >  6:218252/1‑149 (MQ=255)                                                                                                         acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                        <  5:84271/149‑1 (MQ=255)                                                                                                            actgctgCGCCACACCCCCGGCATCCATGCCCACCGCCAGACATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATTC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAg                                                                      <  6:178660/115‑1 (MQ=255)                                                                                                                ctgcGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAcc                                  <  2:187231/148‑1 (MQ=255)                                                                                                                   cgcCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGa                                <  8:114595/147‑1 (MQ=255)                                                                                                                     ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCTGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGCGCATACGGCCGTTGTTCACCGACcac                            <  1:8680/149‑1 (MQ=255)                                                                                                                     ccACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACca                             >  1:238092/1‑148 (MQ=255)                                                                                                                      cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                           >  7:262841/1‑149 (MQ=255)                                                                                                                      cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcac                            >  6:210786/1‑148 (MQ=255)                                                                                                                      cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcac                            <  5:210786/148‑1 (MQ=255)                                                                                                                         acCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGt                        <  2:359390/149‑1 (MQ=255)                                                                                                                           ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                      <  7:88151/149‑1 (MQ=255)                                                                                                                            cACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGGTGTTCACCGCCCACACGTAgg                     >  4:85741/1‑149 (MQ=255)                                                                                                                               ggCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCa                  <  6:217806/149‑1 (MQ=255)                                                                                                                                 cAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGcc                                            >  4:88344/1‑121 (MQ=255)                                                                                                                                  aGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                >  6:177221/1‑148 (MQ=255)                                                                                                                                    ccATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCGCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCa              >  1:279316/1‑149 (MQ=255)                                                                                                                                        gCCCACCGCCAGCCAATGCAGGCGGCGCAGGCGCTCAGCAAGAATGTGCATGC‑CCAGCAGGACAGTGATCAGCGGGTTGAAGTAATAGCCCAGGCTCGCCGCGAGCATACGGCCGTTGGTCACCGCCCACACGTAGGGCAGCCGGTTGc         >  4:258280/1‑149 (MQ=255)                                                                                                                                         tgcgccGCCTGCAATGGCTGGCGGTGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGcc        >  4:234620/5‑149 (MQ=255)                                                                                                                                         cccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGcc        >  5:310851/1‑149 (MQ=255)                                                                                                                                           cACCGCCAGACATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCATCGGGTTGATGTAATAGCCCAGGATGGCCTCGAGCATACGGACGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg      >  8:108656/1‑149 (MQ=255)                                                                                                                                            gccGCCTGCAATGGCTGGCGGTGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACTCGTAGGTCAGCCAGTTGCCGGc     <  7:227252/148‑1 (MQ=255)                                                                                                                                               gccagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt  <  5:177221/149‑1 (MQ=255)                                                                                                                                                 cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt  >  1:227413/1‑147 (MQ=255)                                                                                                                                                   gccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                <  7:150264/83‑1 (MQ=255)                                                                                                                                                   gccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAgg                                                                >  8:150264/1‑83 (MQ=255)                                                                                                                                                    ccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                           >  3:340856/1‑119 (MQ=255)                                                                                                                                                    ccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACcaca                           <  4:340856/119‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CAGGCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGC‑CCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAT  >  NC_002947/436751‑437041

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