Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,998,848	T→C	8.5%	Q47R (CAG→CGG)	PP_2622 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,998,848	0	T	C	8.5%	197.2 / 6.5	70	Q47R (CAG→CGG)	PP_2622	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (3/3);  total (36/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 = 9.58e-01

TGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCG  >  NC_002947/2998720‑2998996                                                                                                                                 |                                                                                                                                                        tgcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCa                                                                                                                                   >  3:170907/1‑149 (MQ=255) tgcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCa                                                                                                                                   <  2:162562/149‑1 (MQ=255) tgcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGAATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCa                                                                                                                                   <  1:40400/149‑1 (MQ=255)    gtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTc                                                                                                                                <  4:229432/149‑1 (MQ=255)    gtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCCGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                                 >  5:294942/1‑148 (MQ=255)    gtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCCGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                                 >  5:299872/1‑148 (MQ=255)       ccTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCAAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                              >  8:230388/1‑148 (MQ=255)        cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGTTATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCt                                                                                                                                             <  1:264948/134‑1 (MQ=255)        cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCt                                                                                                                                             >  2:264948/1‑134 (MQ=255)          gTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCagcag‑c                                                                                                                         >  6:179791/1‑149 (MQ=255)            aTGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGt                                                                                                                       >  7:294567/1‑149 (MQ=255)               cgcgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTc                                                                                                                    <  4:181245/149‑1 (MQ=255)                gcgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTc                                                                                                                    >  8:9478/1‑148 (MQ=255)                       gAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTgg                                                                                                             <  7:9478/148‑1 (MQ=255)                        aaGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGgcg                                                                                                           >  2:209774/1‑149 (MQ=255)                                gAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGAt                                                                                                     <  4:291131/147‑1 (MQ=255)                                     tGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCCGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGgcggcg                                                                                              <  6:294942/149‑1 (MQ=255)                                     tGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCCGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGgcggcg                                                                                              <  6:299872/149‑1 (MQ=255)                                      gcggcgGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGgcggcg                                                                                              >  1:77884/1‑148 (MQ=255)                                      gcggcgGATATCCGCGCTAACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGgcggcg                                                                                              >  5:66961/1‑148 (MQ=255)                                       cggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGa                                                                                             <  4:170907/148‑1 (MQ=255)                                        ggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGc                                                                                                                                    <  4:299901/109‑1 (MQ=255)                                        ggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAAc                                                                                           >  8:192395/1‑149 (MQ=255)                                            gATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                                       <  1:148775/149‑1 (MQ=255)                                             atatCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgtg                                                                                      <  3:194939/149‑1 (MQ=255)                                             atatCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgtg                                                                                      <  2:41032/149‑1 (MQ=255)                                             atatCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgtg                                                                                      <  5:179791/149‑1 (MQ=255)                                                 ccGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTggg                                                                                    >  5:101131/1‑147 (MQ=255)                                                   gcgcTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGCACTTGTGGGAAgc                                                                                >  5:101883/1‑149 (MQ=255)                                                    cgcTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAgcg                                                                               <  7:192395/149‑1 (MQ=255)                                                           cTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTccc                                                                        >  2:64948/1‑149 (MQ=255)                                                           cTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTccc                                                                        >  5:307563/1‑149 (MQ=255)                                                             gCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGt                                                                      >  1:163196/1‑149 (MQ=255)                                                                caagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCgg                                                                   <  2:77884/149‑1 (MQ=255)                                                                 aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAAc                                                                                           <  1:63538/124‑1 (MQ=255)                                                                 aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAAc                                                                                           >  2:63538/1‑124 (MQ=255)                                                                 aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCgg                                                                   >  7:67908/1‑148 (MQ=255)                                                                 aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACGTGTGGGAAGCGCTCTCCCGTCggg                                                                  >  8:237445/1‑149 (MQ=255)                                                                   gcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                                >  1:170359/1‑149 (MQ=255)                                                                        tCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCctc                                                           >  6:259602/1‑149 (MQ=255)                                                                            cGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTcgg                                                       <  8:67908/149‑1 (MQ=255)                                                                            cGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCGG‑GC‑GGGCATTCAGCAG‑CGGTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAgcg                                                                               >  8:221842/1‑125 (MQ=255)                                                                              ttGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTcggcg                                                     <  6:247995/149‑1 (MQ=255)                                                                               tGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTccc                                                                        <  2:28471/129‑1 (MQ=255)                                                                               tGCTCACGTCATGCCCATTGAAGTCGCGCACGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTccc                                                                        >  1:28471/1‑129 (MQ=255)                                                                                gCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTcggcgg                                                    <  6:307563/148‑1 (MQ=255)                                                                                  tCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                                       <  6:43443/111‑1 (MQ=255)                                                                                  tCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTgg                                                                                     >  5:150624/1‑113 (MQ=255)                                                                                  tCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTgg                                                                                     <  6:150624/113‑1 (MQ=255)                                                                                   cACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTg                                                 >  1:127775/1‑148 (MQ=255)                                                                                           gCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGATTGATTGAGGCGAACTTGTGGGAAGCTCTCTACCGTCGGGATGCCTCTCGGCGGGTGCACTGGTcc                                        >  2:133994/1‑149 (MQ=25)                                                                                             ccATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGACGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTccac                                      >  6:27455/1‑149 (MQ=255)                                                                                               aTTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTccacca                                    >  7:76184/1‑149 (MQ=255)                                                                                                   aaGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                                <  2:302350/117‑1 (MQ=255)                                                                                                   aaGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                                >  1:302350/1‑117 (MQ=255)                                                                                                      tCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGTGGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGc                              >  1:88694/1‑149 (MQ=255)                                                                                                        gcgcAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGAt                           <  2:170359/149‑1 (MQ=255)                                                                                                        gcgcAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGAt                           >  8:172668/1‑149 (MQ=255)                                                                                                            aGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaa                        >  2:320431/1‑148 (MQ=255)                                                                                                             ggCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaa                      <  2:163196/149‑1 (MQ=255)                                                                                                             ggCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaa                      <  5:259602/149‑1 (MQ=255)                                                                                                              gCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTTGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGATGGCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa                     <  7:230388/149‑1 (MQ=255)                                                                                                                  ccgccAGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTAGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCAt                                   <  7:263183/129‑1 (MQ=37)                                                                                                                  ccgccAGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTAGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCAt                                   >  8:263183/5‑133 (MQ=37)                                                                                                                   gCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCa                <  1:320431/149‑1 (MQ=255)                                                                                                                      cAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACt              >  3:187829/1‑148 (MQ=255)                                                                                                                       aGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCa            <  2:88694/149‑1 (MQ=255)                                                                                                                           ctgctCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGc        <  6:101131/149‑1 (MQ=255)                                                                                                                           atgctCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTTTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGc        <  8:294567/148‑1 (MQ=255)                                                                                                                            tgctCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCg       >  5:94412/1‑149 (MQ=255)                                                                                                                                ctGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGGGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATc   >  2:94040/1‑149 (MQ=255)                                                                                                                                ctGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaa                       >  7:143362/1‑129 (MQ=255)                                                                                                                                ctGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaa                       <  8:143362/129‑1 (MQ=255)                                                                                                                                 tGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCg  <  8:76184/149‑1 (MQ=255)                                                                                                                                 |                                                                                                                                                        TGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAG‑CGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCACTCATTGCGTATCG  >  NC_002947/2998720‑2998996

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