Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,139,631	A→T	5.5%	R164R (CGT→CGA)	PP_5481 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,139,631	0	A	T	5.5%	211.4 / 2.2	72	R164R (CGT→CGA)	PP_5481	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base T (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

GGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAG  >  NC_002947/2139494‑2139777                                                                                                                                          |                                                                                                                                                   gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACt                                                                                                                                             >  4:110049/1‑145 (MQ=255) gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACt                                                                                                                                             <  3:110049/145‑1 (MQ=255)      tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGccca                                                                                                                                    >  5:283366/1‑149 (MQ=255)        ggATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccac                                                                                                                                   <  7:213220/148‑1 (MQ=255)         gATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccac                                                                                                                                   >  1:269478/1‑147 (MQ=255)         gATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccac                                                                                                                                   <  2:269478/147‑1 (MQ=255)          atatCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccca                                                                                                                                <  5:136516/149‑1 (MQ=255)          atatCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccca                                                                                                                                >  6:136516/1‑149 (MQ=255)             tCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCccacccac                                                                                                                               >  1:53494/1‑147 (MQ=255)                 cATTGGCACGGCTATCACCTCCATGCGGGGGGCCACGGCGGGGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCGGAAGCCTTGCACCGGCATTGGCTACGACGCCACTGCCGGCCCACCCACCAgcgc                                                                                                                         <  8:125934/149‑1 (MQ=255)                      gCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTa                                                                                                                    <  3:362758/149‑1 (MQ=255)                       cACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGc                                                                                                                                           <  1:191782/125‑1 (MQ=255)                       cACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGc                                                                                                                                           >  2:191782/1‑125 (MQ=255)                        aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                        >  3:206055/1‑127 (MQ=255)                        aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                        <  4:206055/127‑1 (MQ=255)                            cGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAt                                                                                                              <  6:283366/149‑1 (MQ=255)                             gATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtg                                                                                                             >  1:118071/1‑149 (MQ=255)                                cACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                                          >  3:256329/1‑149 (MQ=255)                                 aCCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                                          >  7:106277/1‑148 (MQ=255)                                   cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                                          >  3:61325/1‑146 (MQ=255)                                   cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgag                                                                                                       <  3:244848/149‑1 (MQ=255)                                        tGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTa                                                                                                  <  2:118071/149‑1 (MQ=255)                                        tGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTa                                                                                                  <  2:53494/149‑1 (MQ=255)                                           tggtggCCACTGCGGTGCTCTACGTTGTTTTCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCg                                                                                                                                         >  1:32616/1‑107 (MQ=39)                                           tggtggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCg                                                                                                                                         <  2:32616/107‑1 (MQ=39)                                             gtggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGc                                                                                               <  5:113554/147‑1 (MQ=255)                                                  cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATg                                                                                        >  5:350172/1‑149 (MQ=255)                                                  cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATg                                                                                        >  3:67216/1‑149 (MQ=255)                                                     tGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      >  4:104435/1‑148 (MQ=255)                                                          tGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGcac                                                                                <  4:256329/149‑1 (MQ=255)                                                            ctctACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGcacca                                                                              <  5:85267/149‑1 (MQ=255)                                                                 cGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCAcgccg                                                                         >  1:288598/1‑149 (MQ=255)                                                                 cGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCAcgccg                                                                         <  4:61325/149‑1 (MQ=255)                                                                 cGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCAcgc                                                                           <  3:351430/147‑1 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                                >  6:315040/1‑148 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                                >  8:35196/1‑148 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                                >  1:121522/1‑148 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccg                                                               >  8:183221/1‑149 (MQ=255)                                                                                cTTGCTGGTCAGCGGCCCCGGCATGCGGGCCTGAAGCCTTACACCGGCATTGGCGACGACGCCACTGCAGGCCCACCCACCAGCGCCACTAGAACATGGAGGAGAAGTAGCCAACGATCCCCAGCACAACACCGCACACCAcccagcc                                                           >  8:232133/1‑148 (MQ=255)                                                                                    cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACcgagccgagc                                                       <  4:67216/148‑1 (MQ=255)                                                                                    cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCa                                                      <  2:288598/149‑1 (MQ=255)                                                                                     tGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAg                                                     >  2:94075/1‑149 (MQ=255)                                                                                     tGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAg                                                     >  2:98992/1‑149 (MQ=255)                                                                                        tCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAgg                                                    >  6:226140/1‑147 (MQ=255)                                                                                          aGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAgg                                                <  4:334537/149‑1 (MQ=255)                                                                                          aGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                                 <  7:232133/148‑1 (MQ=255)                                                                                           gCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAggg                                               >  1:98566/1‑149 (MQ=255)                                                                                             ggCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                             <  3:104435/149‑1 (MQ=255)                                                                                                   ggCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCca                                         >  1:179869/1‑147 (MQ=255)                                                                                                     cATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTTGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGc                                     >  1:240442/1‑149 (MQ=255)                                                                                                            gCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGa                                >  2:95761/1‑147 (MQ=255)                                                                                                            gCCTGAAGCCTTGCACCGGCATTGGAGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACa                              <  5:315040/149‑1 (MQ=255)                                                                                                             ccTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAg                             <  2:362124/149‑1 (MQ=255)                                                                                                              cTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGc                            <  5:226140/149‑1 (MQ=255)                                                                                                                 aaGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTgg                         >  6:223216/1‑149 (MQ=255)                                                                                                                  aGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGt                        <  2:179869/149‑1 (MQ=255)                                                                                                                  aGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGt                        <  2:46553/149‑1 (MQ=255)                                                                                                                    ccTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCg                      >  2:263632/1‑149 (MQ=255)                                                                                                                    ccTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCg                      >  7:42838/1‑149 (MQ=255)                                                                                                                      ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGt                    >  3:105906/1‑149 (MQ=255)                                                                                                                      ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGt                    >  4:235151/1‑149 (MQ=255)                                                                                                                         cACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                             >  7:267051/1‑121 (MQ=255)                                                                                                                         cACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                             <  8:267051/121‑1 (MQ=255)                                                                                                                         cACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGa                  <  6:350172/148‑1 (MQ=255)                                                                                                                          aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATc                <  1:224829/149‑1 (MQ=255)                                                                                                                               cATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTg           >  3:264676/1‑149 (MQ=255)                                                                                                                                aTTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTg           >  7:352356/1‑148 (MQ=255)                                                                                                                                   ggCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGaa        <  7:35196/148‑1 (MQ=255)                                                                                                                                    gCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgag                                                                                                       <  1:370498/52‑1 (MQ=255)                                                                                                                                    gCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgag                                                                                                       >  2:370498/1‑52 (MQ=255)                                                                                                                                    gCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAg       <  2:121522/148‑1 (MQ=255)                                                                                                                                    gCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAg       <  8:106277/148‑1 (MQ=255)                                                                                                                                        cgacgCCACTGCCGGCCCACCCACCAGCGCCAGTTGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAg  <  2:240442/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                   GGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAG  >  NC_002947/2139494‑2139777

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