Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,845,787	T→G	6.6%	S26R (AGT→AGG)	pfeS‑II →	histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,845,787	0	T	G	6.6%	279.4 / 9.6	91	S26R (AGT→AGG)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base G (2/4);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGC  >  NC_002947/1845642‑1845932                                                                                                                                                       |                                                                                                                                                  aCGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGta                                                                                                                                                 >  3:230899/1‑148 (MQ=255)  cGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTtata                                                                                                                                               <  1:110390/149‑1 (MQ=255)   gAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATAc                                                                                                                                              <  7:345856/149‑1 (MQ=255)          aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGTATAACTGAGGc                                                                                                                                       <  3:333760/149‑1 (MQ=255)          aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGTATAACTGAGGc                                                                                                                                       <  4:230899/149‑1 (MQ=255)           gATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGTATAACTGAGGcc                                                                                                                                      <  2:143551/149‑1 (MQ=255)             tCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCgg                                                                                                                                    <  1:226946/149‑1 (MQ=255)             tCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGTATAACTGAGGCCga                                                                                                                                    <  3:246446/149‑2 (MQ=255)               cGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGGGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGc                                                                                                                                   <  5:131482/148‑1 (MQ=255)                gCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGcac                                                                                                                                 <  8:59005/149‑1 (MQ=255)                  aCCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCGTCAGTTATACCTGGGGCCGGCACAt                                                                                                                               <  6:367651/149‑1 (MQ=255)                  aCCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                               <  7:101356/149‑1 (MQ=255)                       gtgGGGCAAGGGTCACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATAc                                                                                                                                              >  8:352335/1‑129 (MQ=255)                       gtgGGGCAAGGGTCACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATAc                                                                                                                                              <  7:352335/129‑1 (MQ=255)                        tgGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGaaa                                                                                                                          <  3:372200/148‑1 (MQ=255)                        tgGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAAc                                                                                                                         <  8:318553/149‑1 (MQ=255)                           ggCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCa                                                                                                                      <  6:173699/149‑1 (MQ=255)                            gCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAAcc                                                                                                                        >  4:238821/1‑146 (MQ=255)                            gCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAAcc                                                                                                                        >  8:77945/1‑146 (MQ=255)                               aGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAAt                                                                                                                  <  7:69836/149‑1 (MQ=255)                                gggTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCCAGTCCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGaa                                                                                                                   <  1:130753/149‑1 (MQ=255)                                        tGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCCGCCTGCTGATGATCGGAGGGAGAAGTAATACCTGGGGCCGGAACATGAAAACCCAGAATGCCtt                                                                                                             >  1:324155/1‑148 (MQ=255)                                         gCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCctct                                                                                                        >  3:162187/1‑149 (MQ=255)                                          cTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCtctc                                                                                                       <  3:172360/149‑1 (MQ=255)                                          cTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCtctc                                                                                                       <  5:282812/149‑1 (MQ=255)                                             aGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCg                                                                                                      >  2:45949/1‑147 (MQ=255)                                             aGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCg                                                                                                      <  6:297639/147‑1 (MQ=255)                                             aGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACATGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCg                                                                                                      >  1:148687/1‑147 (MQ=255)                                                  cAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCgg                                                                                               <  5:229372/149‑1 (MQ=255)                                                   aGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGc                                                                                              >  3:32773/1‑149 (MQ=255)                                                      gAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACg                                                                                           >  4:58095/1‑149 (MQ=255)                                                           ggAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGa                                                                                      >  7:270173/1‑149 (MQ=255)                                                           ggAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGa                                                                                      <  6:330815/149‑1 (MQ=255)                                                               tGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCtt                                                                                  >  4:280870/1‑149 (MQ=255)                                                                    cATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                              >  3:285676/1‑100 (MQ=255)                                                                    cATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                              <  4:285676/100‑1 (MQ=255)                                                                             gACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                        <  7:77945/147‑1 (MQ=255)                                                                             gACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATg                                                                       <  4:188397/148‑1 (MQ=255)                                                                               ccGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccg                                                                    >  8:222065/1‑149 (MQ=255)                                                                                cGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgc                                                                   <  6:355575/149‑1 (MQ=255)                                                                                cGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgc                                                                   >  6:189580/1‑149 (MQ=255)                                                                                  ccATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgcc                                                                  >  2:77271/1‑148 (MQ=255)                                                                                        gCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCt                                                            <  3:238821/148‑1 (MQ=255)                                                                                        gCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTg                                                           <  2:148687/149‑1 (MQ=255)                                                                                         cTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTg                                                           >  4:15478/1‑148 (MQ=255)                                                                                          tGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGaa                                                         <  8:163563/149‑1 (MQ=255)                                                                                           gTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGaa                                                         >  6:7601/1‑148 (MQ=255)                                                                                                 gAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCt                                                  >  5:97536/1‑149 (MQ=255)                                                                                                  taGCTGGCCATCTTGCTGGTGGGCTTCCGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTg                                                 <  2:324155/148‑1 (MQ=255)                                                                                                       ggCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGACGAAGCtt                                                           >  8:180235/1‑133 (MQ=255)                                                                                                         ccATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                          >  7:288683/1‑149 (MQ=255)                                                                                                          cATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTg                                         >  1:279945/1‑149 (MQ=255)                                                                                                          cATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTg                                         >  7:85072/1‑149 (MQ=255)                                                                                                             cTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcg                                      >  4:351666/1‑149 (MQ=255)                                                                                                             cTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcg                                      >  3:339143/1‑149 (MQ=255)                                                                                                              ttGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcgg                                     >  2:105394/1‑149 (MQ=255)                                                                                                                gCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcga                                     <  2:279945/147‑2 (MQ=255)                                                                                                                  tggtggGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGgcgggcgg                                 >  3:310244/1‑149 (MQ=255)                                                                                                                     tggGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGcgggcggg                                >  5:85132/1‑147 (MQ=255)                                                                                                                     tggGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCa                              <  7:222065/149‑1 (MQ=255)                                                                                                                      gggCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAg                             >  6:156744/1‑149 (MQ=255)                                                                                                                      gggCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAg                             <  3:58095/149‑1 (MQ=255)                                                                                                                        gCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAggg                           <  5:7601/149‑1 (MQ=255)                                                                                                                          ttCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCt                         >  1:31581/1‑149 (MQ=255)                                                                                                                           tCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTg                        >  8:354370/1‑149 (MQ=255)                                                                                                                               cctgcTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAACCCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGAc                     <  7:139947/148‑1 (MQ=255)                                                                                                                               cctgcTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTgg                       <  2:82697/146‑1 (MQ=255)                                                                                                                                ctgctgATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGAc                     >  3:239967/1‑147 (MQ=255)                                                                                                                                 tgctgATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAg                  >  3:47317/1‑149 (MQ=255)                                                                                                                                  gctgATGCTC‑‑‑GGCCTCAGTTATACCTTGGGCCGGAAGATGGGAACCCGGAATGCGTTCCTCTCGGAGCCGGCACGAAGGACCTTGCGGGGATATGCCGCCGATGCTGGACAGGCCTGGAGCAGAGGCGGCCGGGCAGGGCTGGAGCAgt                 >  8:325874/1‑149 (MQ=255)                                                                                                                                   ctgATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtg                <  4:47317/149‑1 (MQ=255)                                                                                                                                   ctgATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCACCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtg                <  1:45949/149‑1 (MQ=255)                                                                                                                                     gatgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtggg              >  2:1118/1‑149 (MQ=255)                                                                                                                                     gatgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtggg              >  4:267936/1‑149 (MQ=255)                                                                                                                                     gatgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtggg              >  2:7384/1‑149 (MQ=255)                                                                                                                                     gatgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTAGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCActggg              >  5:249601/1‑149 (MQ=255)                                                                                                                                      tgtgcc‑‑‑GGCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtgggt             >  3:195454/6‑149 (MQ=255)                                                                                                                                      atgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtgggt             <  8:85072/149‑1 (MQ=255)                                                                                                                                       tgatC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtgggtg            >  6:259627/1‑149 (MQ=255)                                                                                                                                         atC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                          <  4:83623/117‑1 (MQ=255)                                                                                                                                         atC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                          >  3:83623/1‑117 (MQ=255)                                                                                                                                          tC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcg         >  1:87019/1‑149 (MQ=255)                                                                                                                                           c‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcgg        >  2:84661/1‑149 (MQ=255)                                                                                                                                                 ccTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCa     <  1:77271/149‑1 (MQ=255)                                                                                                                                                 ccTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCa     >  5:149677/1‑149 (MQ=255)                                                                                                                                                  cTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAt    <  6:85132/149‑1 (MQ=255)                                                                                                                                                    cAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGc  <  8:288683/149‑1 (MQ=255)                                                                                                                                                     cgTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGc  <  4:32773/147‑1 (MQ=255)                                                                                                                                                     aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGcgggcggg                                <  5:340090/118‑1 (MQ=255)                                                                                                                                                     aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGcgggcggg                                >  6:340090/1‑118 (MQ=255)                                                                                                                                                     aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGc  <  7:354370/148‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGC  >  NC_002947/1845642‑1845932

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