Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,549,906	A→T	5.4%	E102D (GAA→GAT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,549,906	0	A	T	5.4%	270.3 / 2.7	87	E102D (GAA→GAT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base T (2/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCT  >  NC_002947/1549759‑1550051                                                                                                                                                     |                                                                                                                                                  cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAg                                                                                                                                                  <  4:21309/149‑1 (MQ=255) cGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAg                                                                                                                                                  <  4:14834/149‑1 (MQ=255)       cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCt                                                                                                                                              >  5:19727/1‑147 (MQ=255)       cGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCt                                                                                                                                              <  6:19727/147‑1 (MQ=255)         tCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGaaa                                                                                                                                          <  5:310385/149‑1 (MQ=255)           gcaaaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCTTCGTTGACG‑ATGGCCTGagcc                                                                                                                                         <  8:280716/145‑4 (MQ=255)            ccAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGcc                                                                                                                                       >  6:355162/1‑149 (MQ=255)             cAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCg                                                                                                                                      >  5:91990/1‑149 (MQ=255)              aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAGGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGcc                                                                                                                                               >  6:142582/1‑139 (MQ=255)              aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCACGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGcc                                                                                                                                               <  5:142582/139‑1 (MQ=255)              aaaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGt                                                                                                                                     >  2:50027/1‑149 (MQ=255)               aaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTc                                                                                                                                    <  5:347300/149‑1 (MQ=255)               aaGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTc                                                                                                                                    <  5:39659/149‑1 (MQ=255)                   aaTCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGACGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTg                                                                                                                                <  1:221299/149‑1 (MQ=255)                       gAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCg                                                                                                                             <  6:364134/148‑1 (MQ=255)                         gCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCa                                                                                                                          <  2:13459/149‑1 (MQ=255)                            gAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATg                                                                                                                        >  6:181375/1‑148 (MQ=255)                             aaaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCACCG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGa                                                                                                                       <  1:50027/148‑1 (MQ=255)                              aaGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAAc                                                                                                                     <  1:271350/149‑1 (MQ=255)                                   gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAt                                                                                                                >  1:407702/1‑149 (MQ=255)                                   gCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGAt                                                                                                                >  4:61207/1‑149 (MQ=255)                                    cTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                               <  1:264182/149‑1 (MQ=255)                                        cGAACACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGAc                                                                                                            <  6:91990/148‑1 (MQ=255)                                         gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                                 >  8:5070/1‑142 (MQ=255)                                         gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGa                                                                                                                 <  7:5070/142‑1 (MQ=255)                                         gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                               >  3:177312/1‑144 (MQ=255)                                         gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                               <  4:177312/144‑1 (MQ=255)                                         gAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                                          >  7:267779/1‑149 (MQ=255)                                          aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                                          >  1:169079/1‑148 (MQ=255)                                          aaaaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCt                                                                                                          >  1:115501/1‑148 (MQ=255)                                            aaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGTTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                                                                                                       <  3:61207/149‑1 (MQ=255)                                            aaCATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGaa                                                                                                       <  8:267779/149‑1 (MQ=255)                                                  ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACGCAAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAgcgc                                                                                                   <  1:307890/148‑1 (MQ=255)                                                  ggCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGGCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGcgcg                                                                                                  <  1:112354/148‑1 (MQ=255)                                                   gcgcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCgg                                                                                                 >  8:314499/1‑148 (MQ=255)                                                     gcCCAGCTGGTCATGGTAGTTGCTTCCTAGGCCAACGACGGTGCCGGTGACGGCAACACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTAt                                                                                              <  8:327713/149‑1 (MQ=255)                                                     gcCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTAt                                                                                              >  7:7426/1‑149 (MQ=255)                                                           cTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                                        <  5:385460/149‑1 (MQ=255)                                                              gTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                                     >  2:10031/1‑149 (MQ=255)                                                              gTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGc                                                                                     >  7:267728/1‑149 (MQ=255)                                                                cAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCa                                                                                   >  7:131599/1‑149 (MQ=255)                                                                 aacgAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCAc                                                                                  <  7:257276/146‑1 (MQ=255)                                                                  aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACa                                                                                         >  7:280184/1‑141 (MQ=255)                                                                  aGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACa                                                                                         <  8:280184/141‑1 (MQ=255)                                                                       gTTGCTTCCAAGGCCACCGACGCGGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCCGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCCTGAACCCGATGGACCGTAAGCGCGGTATCGACAAGGCCACtgctgc                                                                            <  6:421148/149‑1 (MQ=255)                                                                                      aaCGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCCCTGCTGCCGTCGTTGCCGAAcg                                                             >  7:175686/1‑148 (MQ=255)                                                                                      aaCGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACt                                                             <  2:169079/149‑1 (MQ=255)                                                                                       acgacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACt                                                             >  8:315350/1‑148 (MQ=255)                                                                                         gacgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaa                                                          <  7:315350/149‑1 (MQ=255)                                                                                           cgCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaag                                                         >  5:83046/1‑148 (MQ=255)                                                                                               gCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCTCGGTATCGACAAGGCCActg                                                                                <  3:73643/121‑1 (MQ=255)                                                                                               gCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCActg                                                                                >  4:73643/1‑121 (MQ=255)                                                                                                 cGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGt                                                  <  8:267728/149‑1 (MQ=255)                                                                                                 cGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGt                                                  <  7:314499/149‑1 (MQ=255)                                                                                                    tGACGGCAGCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGAGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCa                                               <  1:201767/149‑1 (MQ=255)                                                                                                      aCGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                              <  2:115501/148‑1 (MQ=255)                                                                                                       cGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAg                                             <  5:181375/148‑1 (MQ=255)                                                                                                        ggCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGc                                            >  1:352856/1‑148 (MQ=255)                                                                                                          caccaccacCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCAt                                         <  2:352856/149‑1 (MQ=255)                                                                                                          caccaccacCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCAt                                         >  5:161008/1‑149 (MQ=255)                                                                                                               ccaccGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGAATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                                    <  3:17746/149‑1 (MQ=255)                                                                                                                caccGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGcc                                    >  3:155955/1‑148 (MQ=255)                                                                                                                  ccGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGa                                  >  3:396907/1‑148 (MQ=255)                                                                                                                  ccGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGAc                                 <  8:131599/149‑1 (MQ=255)                                                                                                                     cTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTc                               >  6:326732/1‑148 (MQ=255)                                                                                                                        ccGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAg                           <  4:155955/149‑1 (MQ=255)                                                                                                                             cTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCa                       <  1:10031/148‑1 (MQ=255)                                                                                                                                gCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                                        >  7:192185/1‑112 (MQ=255)                                                                                                                                gCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTgaaga                                                        <  8:192185/112‑1 (MQ=255)                                                                                                                                gCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGc                   >  6:138729/1‑149 (MQ=255)                                                                                                                                  tCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATg                                                                                                               >  6:46899/1‑55 (MQ=255)                                                                                                                                  tCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGcc                  >  6:183814/1‑148 (MQ=255)                                                                                                                                    aGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc                 <  5:183814/147‑1 (MQ=255)                                                                                                                                    aGGCAATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGCCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc                 <  5:355162/147‑1 (MQ=255)                                                                                                                                     ggCCTTCGTTGACG‑ATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                              >  7:169671/1‑117 (MQ=255)                                                                                                                                     ggCCTTCGTTGACG‑ATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCaa                                              <  8:169671/117‑1 (MQ=255)                                                                                                                                      gCCTTCGTTGACG‑ATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCagg              >  4:215367/1‑148 (MQ=255)                                                                                                                                      gCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                                        <  8:15914/74‑1 (MQ=255)                                                                                                                                      gCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACaa                                                                                        >  7:15914/1‑74 (MQ=255)                                                                                                                                      gCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACGGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc                 <  6:1463/145‑1 (MQ=255)                                                                                                                                      gCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACGGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGccc                 >  5:1463/1‑145 (MQ=255)                                                                                                                                           cgttgaCG‑ATGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCaggtagg          <  3:215367/142‑1 (MQ=255)                                                                                                                                            gTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAggtaggta        <  6:83046/148‑1 (MQ=255)                                                                                                                                            gTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAc       >  8:134707/1‑149 (MQ=255)                                                                                                                                             tCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAc       >  8:242768/1‑148 (MQ=255)                                                                                                                                             tCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCAATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTAc       >  4:297191/1‑148 (MQ=255)                                                                                                                                                 cG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCAtct  >  1:312137/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGTTTCCGTCGCCAAAGAAATCGAGCTGAAAGACGCTTTCGAAAACATGGGCGCCCAGCTGGTCAAGGAAGTTGCTTCCAAGGCCAACGACGCTGCCGGTGACGGCACCACCACCGCTACCGTTCTGGCTCAGGCCATCGTCAACG‑AAGGCCTGAAAGCCGTCGCTGCCGGCATGAACCCGATGGACCTGAAGCGCGGTATCGACAAGGCCACTGCTGCCGTCGTTGCCGAACTGAAGAACCTGTCCAAGCCATGCGCCGACTCCAAGGCCATCGCCCAGGTAGGTACCATCT  >  NC_002947/1549759‑1550051

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