Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,184,117	C→T	5.9%	A85V (GCC→GTC)	PP_4564 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,184,117	0	C	T	5.9%	256.4 / 3.0	85	A85V (GCC→GTC)	PP_4564	lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base T (3/2);  total (42/43)
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

GTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAG  >  NC_002947/5183969‑5184262                                                                                                                                                     |                                                                                                                                                   gTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGc                                                                                                                                                    <  1:14552/149‑1 (MQ=255) gTGCCGTGGGCGCTGAGCGACTTCCTGTACATGAACCGCGAGCTAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGc                                                                                                                                                    <  4:344629/149‑1 (MQ=255)          tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCa                                                                                                                                           <  5:287601/149‑1 (MQ=255)          tCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCa                                                                                                                                           <  5:120388/149‑1 (MQ=255)            gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTATCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAg                                                                                                                                          <  1:307518/148‑1 (MQ=255)            gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAgg                                                                                                                                         >  2:410300/1‑149 (MQ=255)             cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        >  6:115422/1‑149 (MQ=255)                aGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                       <  5:103607/147‑1 (MQ=255)                aGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                                      <  3:21023/148‑1 (MQ=255)                 gCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCAcc                                                                                                                                    <  5:199509/149‑1 (MQ=255)                    aCTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCTGCTACGTGACATGGTCGCCGACGACCACGTCACCgcc                                                                                                                                 >  1:333904/1‑147 (MQ=255)                    aCTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑c                                                                                                                                <  2:54186/149‑1 (MQ=255)                        ccTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCa                                                                                                                              <  4:63635/147‑1 (MQ=255)                            gACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGccc                                                                                                                         <  1:112659/148‑1 (MQ=255)                             aCATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                       <  1:111365/149‑1 (MQ=255)                              cATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGCCCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                       >  8:311079/1‑148 (MQ=255)                               aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                       >  2:92227/1‑147 (MQ=255)                               aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                       >  6:439505/1‑147 (MQ=255)                                tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                            <  8:129294/141‑1 (MQ=255)                                tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                            >  7:129294/1‑141 (MQ=255)                                 gAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGACTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgcag                                                                                                                    >  5:95098/1‑148 (MQ=255)                                  aaCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgcagc                                                                                                                   >  8:23099/1‑148 (MQ=255)                                   aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCaa                                                                                                                 <  5:115422/149‑1 (MQ=255)                                    ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        <  1:294258/126‑1 (MQ=255)                                    ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                        >  2:294258/1‑126 (MQ=255)                                    ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCGACGTGACATGGGCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCaaa                                                                                                                <  6:117039/149‑1 (MQ=255)                                     cgcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAg                                                                                                               >  7:169995/1‑149 (MQ=255)                                       cgtgctaaGGGCGCTTCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTGGCGGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAATGCa                                                                                                             <  2:333904/143‑1 (MQ=255)                                                     tgctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                               >  6:377288/1‑149 (MQ=255)                                                       ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTg                                                                                              <  2:14498/148‑1 (MQ=255)                                                       ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                             <  6:159571/149‑1 (MQ=255)                                                        tgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                             <  1:290632/149‑1 (MQ=255)                                                         ggATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                           >  2:385769/1‑149 (MQ=255)                                                         ggATGATCGGTTGCGCGAGCACCTGACCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                           >  5:203104/1‑149 (MQ=255)                                                          gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgca                                                                                                                     >  2:307170/1‑122 (MQ=255)                                                          gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTgca                                                                                                                     <  1:307170/122‑1 (MQ=255)                                                          gatgatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAgg                                                                                           <  1:385769/148‑1 (MQ=255)                                                                 ggTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGCCGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCCAGGTCCGTGAGGCAATTGGt                                                                                   >  8:348885/1‑149 (MQ=255)                                                                 ggTGGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGt                                                                                   <  4:300735/149‑1 (MQ=255)                                                                  gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGGCCGGGTAACTAGAGTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTGACCGA‑CCAGGCCGTGCAGCAAAGCACCCGCGAGGTCGGTGAGGCAATTGGTc                                                                                  >  7:272786/1‑149 (MQ=255)                                                                  gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                  >  1:188588/1‑149 (MQ=255)                                                                  gTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGACGCAATTGGTc                                                                                  >  7:270559/1‑149 (MQ=255)                                                                      cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                              <  7:23099/149‑1 (MQ=255)                                                                      cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                              <  5:109355/149‑1 (MQ=255)                                                                      cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                              <  4:178382/149‑1 (MQ=255)                                                                       gcgAGCACCTGGCCTGGCATTGCAAGACTCCACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGGGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTg                                                                             <  4:242945/149‑1 (MQ=255)                                                                              cctggcctggCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                      <  2:301517/149‑1 (MQ=255)                                                                               ctggcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                      >  8:84277/1‑148 (MQ=255)                                                                                  gcctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGACTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGaaa                                                                   <  6:95098/148‑1 (MQ=255)                                                                                    ctggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCAACGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCa                                                                >  2:189837/1‑149 (MQ=255)                                                                                     tggcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                               >  7:19058/1‑149 (MQ=255)                                                                                       gcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAccc                                                             >  5:127037/1‑149 (MQ=255)                                                                                       gcATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAcc                                                              <  1:410300/148‑1 (MQ=255)                                                                                            gCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCg                                                        >  8:371213/1‑149 (MQ=255)                                                                                             cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGTCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                        <  1:44089/132‑1 (MQ=255)                                                                                             cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGTCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGa                                                                        >  2:44089/1‑132 (MQ=255)                                                                                             cAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGGCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCgg                                                       <  7:128058/149‑1 (MQ=255)                                                                                              aaGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGt                                                                                               >  6:167241/1‑108 (MQ=255)                                                                                                  cTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAcc                                                  >  5:366026/1‑149 (MQ=255)                                                                                                     aaCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGaa                                               >  3:331454/1‑149 (MQ=255)                                                                                                      aCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                              <  6:234554/149‑1 (MQ=255)                                                                                                      aCTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                              <  6:240598/149‑1 (MQ=255)                                                                                                        tGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCAAGGAAATCACCCCTTCGGCCACCGAActg                                            >  1:61530/1‑149 (MQ=255)                                                                                                           ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGACctgctg                                         >  2:87434/1‑149 (MQ=255)                                                                                                           ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAActgctg                                         >  3:21963/1‑149 (MQ=255)                                                                                                            cGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACtgctgc                                        >  1:328344/1‑149 (MQ=255)                                                                                                            cGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAAACACCCCTTCGGCCACCGAACtgctgc                                        >  1:65327/1‑149 (MQ=255)                                                                                                                   cTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                                 <  4:21963/149‑1 (MQ=255)                                                                                                                    tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTCCTGCGCGGCAt                                 >  4:315705/1‑148 (MQ=255)                                                                                                                       aCTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCACGAGc                             <  5:377288/149‑1 (MQ=255)                                                                                                                           gCTCGACCGGCTACGTGACAGGGTCGCCGACGCCCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGa                           <  8:307155/147‑1 (MQ=255)                                                                                                                             tCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAc                       >  7:375136/1‑149 (MQ=255)                                                                                                                               gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAcc                      <  5:364971/148‑1 (MQ=255)                                                                                                                               gACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAAcc                      >  6:364971/1‑148 (MQ=255)                                                                                                                                  ctggtcggTGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGt                  >  1:433575/8‑149 (MQ=39)                                                                                                                                     cttcGTTACATGGTCGACGACGACATGGTCAGCGA‑GCAGGGCCTGCAGCATACCAGCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGATATGACCGCTTCCGCCACGGAACTGCTGCGCGGCATGAGCGACAAGCAGGTAgc               <  8:270559/146‑1 (MQ=255)                                                                                                                                     cTACGTGACATGGTCGCCGACGTCCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCTATTGTCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgc               <  8:272786/149‑1 (MQ=255)                                                                                                                                        cgggACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa            <  2:179340/146‑1 (MQ=255)                                                                                                                                        cGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg              <  5:439505/147‑1 (MQ=255)                                                                                                                                         gTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCACCAAAGCCCCCGCGACGTCCGTGAGGCAATTCGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTCCGCGGCATCAGAGACAACCAGGAAgcc              >  4:361579/1‑145 (MQ=255)                                                                                                                                          tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg              >  1:185646/1‑145 (MQ=255)                                                                                                                                          tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg                >  7:299233/1‑143 (MQ=255)                                                                                                                                          tGACCTGGTCGTCGGCGACCATGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg                <  8:299233/143‑1 (MQ=255)                                                                                                                                          tGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcg              >  3:182411/1‑145 (MQ=255)                                                                                                                                          tGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGa          <  6:366026/149‑1 (MQ=255)                                                                                                                                            aCATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGAt         <  8:375136/148‑1 (MQ=255)                                                                                                                                                  tCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAg  >  4:132335/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GTGCCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAG  >  NC_002947/5183969‑5184262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑