Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,184,120	A→G	5.9%	D86G (GAC→GGC)	PP_4564 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,184,120	0	A	G	5.9%	241.2 / 4.1	85	D86G (GAC→GGC)	PP_4564	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (38/41);  new base G (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

TCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTT  >  NC_002947/5183978‑5184267                                                                                                                                               |                                                                                                                                                     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                                                                                                                            >  6:439505/1‑147 (MQ=255)                      aTGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTg                                                                                                                            >  2:92227/1‑147 (MQ=255)                       tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                                 >  7:129294/1‑141 (MQ=255)                       tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAgg                                                                                                                                 <  8:129294/141‑1 (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                                                                                                                                             >  2:294258/1‑126 (MQ=255)                           ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                             <  1:294258/126‑1 (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                                                                                   <  4:178382/149‑1 (MQ=255)                                                             cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                                   <  7:23099/149‑1 (MQ=255)                                                             cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                                   <  5:109355/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)                                                                                                                                             cgacgaCCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCt   <  8:19058/149‑1 (MQ=255)                                                                                                                                             cgacgaCCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCt   >  7:255062/1‑149 (MQ=255)                                                                                                                                              gacgacCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGTTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCtt  <  6:345412/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     TCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTGCCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGA‑CCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTT  >  NC_002947/5183978‑5184267

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