Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,139,637	C→A	5.2%	Q162H (CAG→CAT)	PP_5481 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,139,637	0	C	A	5.2%	254.9 / 3.2	77	Q162H (CAG→CAT)	PP_5481	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base A (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCTCGTC  >  NC_002947/2139494‑2139783                                                                                                                                                |                                                                                                                                                   gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACt                                                                                                                                                   <  3:110049/145‑1 (MQ=255) gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACt                                                                                                                                                   >  4:110049/1‑145 (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                                                                                                                                      >  6:136516/1‑149 (MQ=255)          atatCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccca                                                                                                                                      <  5:136516/149‑1 (MQ=255)             tCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCccacccac                                                                                                                                     >  1:53494/1‑147 (MQ=255)                 cATTGGCACGGCTATCACCTCCATGCGGGGGGCCACGGCGGGGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCGGAAGCCTTGCACCGGCATTGGCTACGACGCCACTGCCGGCCCACCCACCAgcgc                                                                                                                               <  8:125934/149‑1 (MQ=255)                      gCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTa                                                                                                                          <  3:362758/149‑1 (MQ=255)                       cACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGc                                                                                                                                                 >  2:191782/1‑125 (MQ=255)                       cACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGc                                                                                                                                                 <  1:191782/125‑1 (MQ=255)                        aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                              <  4:206055/127‑1 (MQ=255)                        aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                              >  3:206055/1‑127 (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                                                                      >  8:35196/1‑148 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                                      >  1:121522/1‑148 (MQ=255)                                                                           gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACaccacc                                                                      >  6:315040/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:98992/1‑149 (MQ=255)                                                                                     tGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAg                                                           >  2:94075/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                                                   <  8:267051/121‑1 (MQ=255)                                                                                                                         cACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                                   >  7:267051/1‑121 (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)                                                                                                                                            gCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGAc                                     >  2:261135/1‑116 (MQ=255)                                                                                                                                            gCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtcg    >  2:58284/1‑149 (MQ=255)                                                                                                                                            gCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACCCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGAc                                     <  1:261135/116‑1 (MQ=255)                                                                                                                                              cACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtcgtc  <  8:352356/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCTCGTC  >  NC_002947/2139494‑2139783

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