Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,652,964	G→A	5.5%	H1269Y (CAT→TAT)	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,652,964	0	G	A	5.5%	225.0 / 2.1	73	H1269Y (CAT→TAT)	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base A (2/2);  total (38/35)
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.25e-01

GTGCCGCCAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGC  >  NC_002947/1652819‑1653111                                                                                                                                                  |                                                                                                                                                    gtgCCGCCAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCCGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCa                                                                                                                                                   <  6:165444/148‑1 (MQ=255)   gccgccAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAgg                                                                                                                                                 <  7:47854/148‑1 (MQ=255)     cgccAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCcg                                                                                                                                              >  8:257933/1‑149 (MQ=255)        cAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCgg                                                                                                                                           <  5:104043/149‑1 (MQ=255)             tGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCCCGCAGGTCACGGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCGACGTTGAGCCGGCCAGGCAGGCCGCGGGTGTc                                                                                                                                      >  8:226614/1‑149 (MQ=255)               cccTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACT‑GTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTcc                                                                                                                                     >  6:254822/1‑147 (MQ=255)                 cTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCgg                                                                                                                                           <  4:241542/140‑1 (MQ=255)                 cTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCgg                                                                                                                                           <  6:311558/140‑1 (MQ=255)                 cTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCgg                                                                                                                                           >  5:311558/1‑140 (MQ=255)                 cTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCgg                                                                                                                                           >  3:241542/1‑140 (MQ=255)                      tCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGcc                                                                                                                             >  7:117006/1‑149 (MQ=255)                      tCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGcc                                                                                                                             >  2:145547/1‑149 (MQ=255)                        cAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCGCGACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGccc                                                                                                                            <  6:286007/148‑1 (MQ=255)                         aGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCtt                                                                                                                          <  8:113489/149‑1 (MQ=255)                          gggTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTg                                                                                                                         >  2:300484/1‑149 (MQ=255)                           ggTGGCACTCTCTATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTgg                                                                                                                        <  1:300484/149‑1 (MQ=255)                           ggTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTgg                                                                                                                        >  5:200072/1‑149 (MQ=255)                           ggTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTgg                                                                                                                        >  6:90227/1‑149 (MQ=255)                            gTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGc                                                                                                                       <  8:117006/149‑1 (MQ=255)                             tGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGcc                                                                                                                      <  6:200072/149‑1 (MQ=255)                                  ctctCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGtt                                                                                                                 >  5:167338/1‑149 (MQ=255)                                   tctcGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGtt                                                                                                                 >  7:308273/1‑148 (MQ=255)                                   tctcGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGtt                                                                                                                 >  7:312947/1‑148 (MQ=255)                                    ctcGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAg                                                                                                                                  >  1:211751/1‑130 (MQ=255)                                    ctcGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAg                                                                                                                                  <  2:211751/130‑1 (MQ=255)                                             ggTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGCGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCGTGGCCATGTTGAAACCGgcgc                                                                                                      >  5:268987/1‑149 (MQ=255)                                              gTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCgcc                                                                                                     <  7:194278/149‑1 (MQ=255)                                                         gCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGt                                                                                          <  7:226614/149‑1 (MQ=255)                                                            ttGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATACAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTc                                                                                         >  3:252005/1‑147 (MQ=255)                                                            ttGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATACAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTc                                                                                         >  3:247320/1‑147 (MQ=255)                                                              gCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAgg                                                                                     >  6:34239/1‑149 (MQ=255)                                                               cAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGt                                                                                                                  <  4:175018/119‑1 (MQ=255)                                                               cAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGt                                                                                                                  >  3:175018/1‑119 (MQ=255)                                                                   tCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTgtc                                                                                <  5:90227/149‑1 (MQ=255)                                                                    cGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTgtcg                                                                               >  6:12037/1‑149 (MQ=255)                                                                          gcTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTgtcgcgtcgc                                                                         >  1:948/1‑149 (MQ=255)                                                                                cgcgcAGGTCACTGTTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAACCCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCctgg                                                                   <  8:308273/149‑1 (MQ=255)                                                                                cgcgcAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCctgg                                                                   <  8:312947/149‑1 (MQ=255)                                                                                 gcgcAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCctggc                                                                  <  1:145547/149‑1 (MQ=255)                                                                                  cgcAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCctggct                                                                 <  7:257933/149‑1 (MQ=255)                                                                                          aCTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAgg                                                         >  5:96030/1‑149 (MQ=255)                                                                                            tGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGaa                                                             >  8:200555/1‑143 (MQ=255)                                                                                            tGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGaa                                                             <  7:200555/143‑1 (MQ=255)                                                                                            tGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCt                                                       >  8:260518/1‑149 (MQ=255)                                                                                              gTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGc                                                        <  8:248536/146‑1 (MQ=255)                                                                                              gTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGc                                                        >  7:248536/1‑146 (MQ=255)                                                                                               tATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTc                                                     >  8:193923/1‑148 (MQ=255)                                                                                                 ttCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAAc                                                  >  5:10910/1‑149 (MQ=255)                                                                                                  accAGGTGAGCTTATCGCGCTTGTCGAGCTCTACCATGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGCCGCAGGCCGAGATGTCGAGGTTGTCGCGTCGCACct                                                                     >  6:329715/2‑129 (MQ=255)                                                                                                                   cgcTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTc                                  >  4:48021/1‑147 (MQ=255)                                                                                                                    gcTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGat                               >  7:121313/1‑148 (MQ=255)                                                                                                                    gcTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAg                               <  6:167338/149‑1 (MQ=255)                                                                                                                     cTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAgg                              >  7:134162/1‑149 (MQ=255)                                                                                                                        gTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGc                             <  5:254822/147‑1 (MQ=255)                                                                                                                         tCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTc                          >  6:367858/1‑149 (MQ=255)                                                                                                                            aGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAg                       >  5:199666/1‑149 (MQ=255)                                                                                                                            aGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAg                       <  3:48021/149‑1 (MQ=255)                                                                                                                            aGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAg                       >  3:245249/1‑149 (MQ=255)                                                                                                                              ctctACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAggg                     >  3:20566/1‑149 (MQ=255)                                                                                                                                ctACCTTGAGCCGGCCATACAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTg                   <  4:252005/149‑1 (MQ=255)                                                                                                                                atACCTTGAGCCTGCCATACAGGCCGCGGGTGGCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTg                   <  4:247320/148‑1 (MQ=255)                                                                                                                                 tACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGa                  >  6:278710/1‑149 (MQ=255)                                                                                                                                  aCCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGAt                 >  1:197541/1‑149 (MQ=255)                                                                                                                                   ccTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGAt                 <  5:34239/148‑1 (MQ=255)                                                                                                                                    cTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGa                        <  7:84904/140‑1 (MQ=255)                                                                                                                                    cTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGa                        >  8:84904/1‑140 (MQ=255)                                                                                                                                       gAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTg              <  8:91701/147‑1 (MQ=255)                                                                                                                                        aGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTc           <  5:278710/149‑1 (MQ=255)                                                                                                                                           cGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGc        <  6:96030/149‑1 (MQ=255)                                                                                                                                           cGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGc        <  4:20566/149‑1 (MQ=255)                                                                                                                                               cATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTTATCTGTTCGGCCGAg    <  8:121313/149‑1 (MQ=255)                                                                                                                                               cATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAg    <  2:197541/149‑1 (MQ=255)                                                                                                                                                 tGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGc  <  4:245249/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTGCCGCCAATATGCCCTGCATCCAGGGTGGCACTCTCGATGCGGGTATCGCCGCGGCTTTGCAGGTCGATGCGCTCGGCGCGCAGGTCACTGGTATTCCAGGTGAGGTTATCGCGCTTGTCGAGCTCTACCTTGAGCCGGCCATGCAGGCCGCGGGTGTCCAAGGCGCCCTTGGCCATGTTGAAACCGGCGCCGGCCGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGC  >  NC_002947/1652819‑1653111

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