Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,228,100	C→G	6.0%	R313R (CGC→CGG)	PP_2828 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,228,100	0	C	G	6.0%	210.3 / 3.2	66	R313R (CGC→CGG)	PP_2828	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/29);  new base G (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

TGATCGCTATTGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCTGCTGTTCAGC  >  NC_002947/3227952‑3228247                                                                                                                                                     |                                                                                                                                                     tGATCGCTATTGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTgcgc                                                                                                                                                      <  2:153996/149‑1 (MQ=255)     cGCTATTGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGccc                                                                                                                                                  >  2:25179/1‑149 (MQ=255)       cTATTGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGTCGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCg                                                                                                                                                 >  1:286476/1‑148 (MQ=255)          ttGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCa                                                                                                                                               <  1:244931/147‑1 (MQ=255)              tgaCCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTg                                                                                                                                          >  3:126608/1‑148 (MQ=255)                 ccTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑Gc                                                                                                                                     >  6:51879/1‑149 (MQ=255)                    gTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGc                                                                                                                                  >  5:239048/1‑149 (MQ=255)                        ccGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGACACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTgg                                                                                                                               <  1:297782/148‑1 (MQ=255)                         cGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGcc                                                                                                                             >  1:114786/1‑149 (MQ=255)                          gCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCg                                                                                                                            >  6:138246/1‑149 (MQ=255)                          gCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCg                                                                                                                            <  1:25179/149‑1 (MQ=255)                             cttgcttgAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGgc                                                                                                                         <  5:49545/149‑1 (MQ=255)                              ttgcttgAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGgcg                                                                                                                        >  3:356661/1‑149 (MQ=255)                              ttgcttgAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGgcg                                                                                                                        >  5:373594/1‑149 (MQ=255)                                gcttgAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGgc                                                                                                                         <  6:244804/146‑1 (MQ=255)                                gcttgAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGgc                                                                                                                         >  5:244804/1‑146 (MQ=255)                                     aaGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAgc                                                                                                                  <  4:356661/148‑1 (MQ=255)                                       gagCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCt                                                                                                               >  3:389294/1‑149 (MQ=255)                                         gcagcagcTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGa                                                                                                             <  1:151286/149‑1 (MQ=255)                                         gcagcagcTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGa                                                                                                             <  6:239048/149‑1 (MQ=255)                                          cagcagcTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAg                                                                                                            >  1:56205/1‑149 (MQ=255)                                             cagcTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCt                                                                                                         <  2:56205/149‑1 (MQ=255)                                               gcTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTgg                                                                                                       <  7:15474/149‑1 (MQ=255)                                                   gACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTggcagg                                                                                                   >  7:126699/1‑149 (MQ=255)                                                   gACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTggca                                                                                                     <  2:286476/147‑1 (MQ=255)                                                      gcgATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGCCTggcaggg                                                                                                  >  4:61172/1‑147 (MQ=255)                                                        gATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGgcagggcagg                                                                                              <  3:112461/149‑1 (MQ=255)                                                            ctaggTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCggg                                                                                          >  8:322267/4‑149 (MQ=255)                                                             tGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGGACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGt                                                                                         >  3:97567/1‑149 (MQ=255)                                                               ggTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGt                                                                                         <  1:3080/147‑1 (MQ=255)                                                                     cACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGgcagc                                                                                   >  2:212489/1‑147 (MQ=255)                                                                      aCTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCagcaag                                                                                >  6:267726/1‑149 (MQ=255)                                                                      aCTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCagcaa                                                                                 <  5:138246/148‑1 (MQ=255)                                                                        ttCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCagcaag                                                                                >  8:200498/1‑147 (MQ=255)                                                                              aCGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGACGCCGCTGGCCGGGCGCCCAGCGCTGAGGCTggcaggg                                                                                                  <  1:72696/124‑1 (MQ=255)                                                                              aCGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACCGCACAGCAGCATGAACTGCGCGCCCCCAGTCTGGA‑GCCGCTGGCCGGGCCCCCAGCGCTCAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGa                                                                        >  4:343580/1‑149 (MQ=255)                                                                              aCGGCTGCGGATGCTTTGCAGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTggcaggg                                                                                                  >  2:72696/1‑123 (MQ=255)                                                                                ggCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGgcagc                                                                                   >  7:107953/1‑136 (MQ=255)                                                                                ggCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGgcagc                                                                                   <  8:107953/136‑1 (MQ=255)                                                                                     cGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACgccgtg                                                                 <  7:208995/149‑1 (MQ=255)                                                                                     cGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACgccgtg                                                                 <  1:267520/149‑1 (MQ=255)                                                                                          gCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGtt                                                            >  6:198877/1‑149 (MQ=255)                                                                                           cTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTg                                                           <  1:25544/149‑1 (MQ=255)                                                                                            tttGCCGCTGTCCGTATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCgg                                                                                           <  5:3934/116‑1 (MQ=255)                                                                                            tttGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCgg                                                                                           >  6:3934/1‑116 (MQ=255)                                                                                             ttGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGc                                                          >  8:144253/1‑148 (MQ=255)                                                                                                   cTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTgg                                                   <  5:267726/149‑1 (MQ=255)                                                                                                   cTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTgg                                                   <  4:389294/149‑1 (MQ=255)                                                                                                    tGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGc                                                  <  2:76906/149‑1 (MQ=255)                                                                                                         gAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGc                                             <  1:212489/149‑1 (MQ=255)                                                                                                         gAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGc                                             >  3:264766/1‑149 (MQ=255)                                                                                                               gtcgACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACgg                                       >  6:243659/1‑149 (MQ=255)                                                                                                                  gACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCa                                     >  4:215070/1‑148 (MQ=255)                                                                                                                   aCTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCAt                                    >  6:145071/1‑148 (MQ=255)                                                                                                                     tCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCTGGTGGCAGCAAGCAGCGTGACGCAGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCaa                                 >  8:234082/1‑149 (MQ=255)                                                                                                                       gCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAAc                                <  7:144253/148‑1 (MQ=255)                                                                                                                           gCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAgg                            <  3:215070/148‑1 (MQ=255)                                                                                                                              aCGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGccc                         >  4:171804/1‑148 (MQ=255)                                                                                                                                   acaGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGAAAGTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTgg                   >  1:134644/1‑149 (MQ=255)                                                                                                                                      gCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGGCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGgtgt                <  3:343580/149‑1 (MQ=255)                                                                                                                                      gCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGgtgt                <  4:264766/149‑1 (MQ=255)                                                                                                                                             gAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCctgct         >  1:79135/1‑149 (MQ=255)                                                                                                                                               aCTGCGGGCGCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGct                                     >  2:165066/1‑118 (MQ=255)                                                                                                                                               aCTGCGGGCGCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGct                                     <  1:165066/119‑2 (MQ=255)                                                                                                                                               aCTGCGGGCGCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCa                                     >  2:298307/1‑119 (MQ=255)                                                                                                                                               aCTGCGGGCGCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCa                                     <  1:298307/119‑1 (MQ=255)                                                                                                                                                cTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCTGCTGtt      <  5:198877/149‑1 (MQ=255)                                                                                                                                                   cgcgcCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCTGCTGTTCAg   >  6:62245/1‑149 (MQ=255)                                                                                                                                                    gcgcCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCTGCTGTTCAGc  <  2:114786/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGATCGCTATTGATGACCTGTTGCCGCACTTGCTTGAAGAGCAGCAGCTTGACGCGATGCTGGGTTGTCACTTCGAGACGGCTGCGGATGCTTTGCCGCTGTCCGAATTCGTCGACTCGCTAGCCACGGCACAGCTGCATGAACTGCGCGCCCGCAGTCTGGA‑GCCGCTGGCCGGGCGCCCAGCGCTGAGGCTGGCAGGGCAGGCGGGTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCTGCTGTTCAGC  >  NC_002947/3227952‑3228247

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