Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,799,180	A→G	5.6%	intergenic (‑194/+22)	PP_2452 ← / ← ansB	hypothetical protein/glutaminase‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,799,180	0	A	G	5.6%	203.2 / 2.5	71	intergenic (‑194/+22)	PP_2452/ansB	hypothetical protein/glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

AAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTG  >  NC_002947/2799033‑2799322                                                                                                                                                    |                                                                                                                                                aagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCAc                                                                                                                                                >  6:71142/1‑149 (MQ=255) aagtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCAGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCaa                                                                                                                                                >  7:65520/1‑148 (MQ=255)  agtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCAcc                                                                                                                                               <  4:289056/149‑1 (MQ=255)   gtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGGCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCg                                                                                                                                              <  4:126607/149‑1 (MQ=255)   gtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCg                                                                                                                                              >  8:65645/1‑149 (MQ=255)         tCGCAACAGCTTGTCAGGAGAAATCCACCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAACCGATGTAGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACTAGGCCCTCCGAGCCGGCGCACCGGCCCTg                                                                                                                                        >  3:230385/1‑149 (MQ=255)            cAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTATACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGCCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                     <  1:326041/149‑1 (MQ=255)            cAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTGAgg                                                                                                                                     <  5:280512/149‑1 (MQ=255)             aaCAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAggg                                                                                                                                    <  8:47737/149‑1 (MQ=255)               cAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGcc                                                                                                                                           >  7:68375/1‑140 (MQ=255)               cAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGcc                                                                                                                                           <  8:68375/140‑1 (MQ=255)                 gCTTGTCCCGAGAAATCCAGCGTTGGCCTGTAAACCATGACGCCGGTCAGAAAGCTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACCCGGTTGAGGAATTCGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGaa                                                                                                                                <  4:92832/149‑1 (MQ=255)                  cTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAg                                                                                                                               >  8:167225/1‑149 (MQ=255)                    tGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCa                                                                                                                             >  3:57528/1‑149 (MQ=255)                        aGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCa                                                                                                                        <  8:88751/149‑1 (MQ=255)                        aGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCa                                                                                                                        <  7:295398/149‑1 (MQ=255)                         ggAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCAGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCGTGAGGGGGAAGCCC‑TCAg                                                                                                                       >  6:262874/1‑149 (MQ=255)                          gagaAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAg                                                                                                                       >  4:182548/1‑148 (MQ=255)                            gaAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGCTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTAc                                                                                                                    >  7:222739/1‑149 (MQ=255)                              aaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTc                                                                                                                  >  1:197866/1‑149 (MQ=255)                                tCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGTGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTccc                                                                                                                <  6:162146/149‑1 (MQ=255)                                           ccTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGc                                                                                                     >  7:139534/1‑149 (MQ=255)                                             tGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTg                                                                                                   >  8:194514/1‑149 (MQ=255)                                                 aaCCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAg                                                                                                <  2:197866/148‑1 (MQ=255)                                                        aGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTg                                                                                         <  8:139534/148‑1 (MQ=255)                                                            cGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGCTcag                                                                                                      >  2:112372/1‑129 (MQ=255)                                                            cGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTc                                                                                    >  2:313235/1‑149 (MQ=255)                                                               tCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAg                                                                                                <  2:43891/134‑1 (MQ=255)                                                               tCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAg                                                                                                >  1:43891/1‑134 (MQ=255)                                                                taGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                <  5:307927/148‑1 (MQ=255)                                                                 aGAAACGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                <  5:71142/148‑1 (MQ=255)                                                                  gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTgg                                                                                >  6:307927/1‑147 (MQ=255)                                                                  gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGt                                                                              >  4:307560/1‑149 (MQ=255)                                                                      ggTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                          <  2:149934/149‑1 (MQ=255)                                                                      ggTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                          <  7:140878/149‑1 (MQ=255)                                                                      ggTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTg                                                                          <  6:136437/149‑1 (MQ=255)                                                                        tAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTgg                                                                         <  7:65645/148‑1 (MQ=255)                                                                         aaaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGt                                                                        >  6:161661/1‑148 (MQ=255)                                                                          aaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTg                                                                                 <  6:162938/138‑1 (MQ=255)                                                                          aaCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTg                                                                                 >  5:162938/1‑138 (MQ=255)                                                                             gATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCg                                                                   >  1:51843/1‑149 (MQ=255)                                                                              aTGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGc                                                                  <  1:18656/149‑1 (MQ=255)                                                                                 tCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTa                                                                <  7:167225/148‑1 (MQ=255)                                                                                         ttGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGc                                                         >  2:117181/1‑147 (MQ=255)                                                                                                 tctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCcgc                                               >  6:6128/1‑149 (MQ=255)                                                                                                     tgtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCtt                                           >  7:135189/1‑149 (MQ=255)                                                                                                      gtCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCttt                                          >  7:260040/1‑149 (MQ=255)                                                                                                         tACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTCCCATCGCCAGGATCCGCGCTTTTTg                                       >  2:72082/1‑149 (MQ=255)                                                                                                         tACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATcc                                                 >  4:236995/1‑139 (MQ=255)                                                                                                         tACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATcc                                                 <  3:236995/139‑1 (MQ=255)                                                                                                          aCGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACAGGACCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATACGCGCTTTTTGc                                      >  1:278166/1‑149 (MQ=255)                                                                                                           cgcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGc                                      <  1:117181/148‑1 (MQ=255)                                                                                                            gctgTTGAGGAATTGGCAACGAGGCCCTCTGAGCCGGCGCACCGTCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGTGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                                    <  2:278166/146‑1 (MQ=255)                                                                                                            gcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                                    <  1:313235/149‑1 (MQ=255)                                                                                                                  gAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCa                               <  8:65520/148‑1 (MQ=255)                                                                                                                     gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTc                           >  5:265228/1‑149 (MQ=255)                                                                                                                     gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTc                           <  3:307560/149‑1 (MQ=255)                                                                                                                       aTTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                         <  5:161661/149‑1 (MQ=255)                                                                                                                       aTTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGt                         <  2:51843/149‑1 (MQ=255)                                                                                                                           gCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGc                     >  1:134805/1‑149 (MQ=255)                                                                                                                            cAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCTCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGcc                    <  8:222739/149‑1 (MQ=255)                                                                                                                                  ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                       <  6:197725/140‑1 (MQ=255)                                                                                                                                  ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                       >  5:197725/1‑140 (MQ=255)                                                                                                                                  ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc              <  5:126844/149‑1 (MQ=255)                                                                                                                                    cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc              >  3:45292/1‑147 (MQ=255)                                                                                                                                     ccTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAg           >  8:215894/1‑149 (MQ=255)                                                                                                                                      cTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAg           >  1:217856/1‑149 (MQ=255)                                                                                                                                         aggtccGGCGCACCGGCGCTGAGGGGGAAGCAC‑TCAGTACTCCCAGATGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCTGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAGCCAGTc         <  4:230385/143‑1 (MQ=255)                                                                                                                                         agggCCGGTGCGCCGGCTCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGt       >  3:309335/4‑149 (MQ=255)                                                                                                                                              cGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg  <  2:217856/149‑1 (MQ=255)                                                                                                                                              cGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg  <  4:57528/149‑1 (MQ=255)                                                                                                                                              cGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTg  <  5:6128/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                AAGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCAC‑TCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTG  >  NC_002947/2799033‑2799322

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