Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,845,783	T→C	7.0%	L25P (CTC→CCC)	pfeS‑II →	histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,845,783	0	T	C	7.0%	246.7 / 6.1	86	L25P (CTC→CCC)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base C (2/4);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAT  >  NC_002947/1845642‑1845930                                                                                                                                                   |                                                                                                                                                    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                                                                                                                                     <  4:230899/149‑1 (MQ=255)          aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCCCAGGTATAACTGAGGc                                                                                                                                     <  3:333760/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                                                                                                    <  6:297639/147‑1 (MQ=255)                                             aGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCg                                                                                                    >  2:45949/1‑147 (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                                                                                    <  6:330815/149‑1 (MQ=255)                                                           ggAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGa                                                                                    >  7:270173/1‑149 (MQ=255)                                                               tGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCtt                                                                                >  4:280870/1‑149 (MQ=255)                                                                    cATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                            <  4:285676/100‑1 (MQ=255)                                                                    cATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                            >  3:285676/1‑100 (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                                       >  7:85072/1‑149 (MQ=255)                                                                                                          cATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTg                                       >  1:279945/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)                                                                                                                                                   |                                                                                                                                                    ACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTC‑‑GACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATC‑‑‑GGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAT  >  NC_002947/1845642‑1845930

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