Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,780,949	T→C	5.5%	W83R (TGG→CGG)	PP_2434 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,780,949	0	T	C	5.5%	218.1 / 2.4	73	W83R (TGG→CGG)	PP_2434	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/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 = 9.25e-01

GCCGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCG  >  NC_002947/2780802‑2781093                                                                                                                                                    |                                                                                                                                                  gccgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTg                                                                                                                                                  >  2:264743/1‑149 (MQ=255) gccgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTg                                                                                                                                                  >  7:44867/1‑149 (MQ=255)   cgccACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGTAGCCCGGCTGGt                                                                                                                                                <  3:148491/149‑1 (MQ=255)     ccaccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCt                                                                                                                                              >  5:258928/1‑149 (MQ=255)        ccaGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCt                                                                                                                                           <  5:245110/149‑1 (MQ=255)          aGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCa                                                                                                                                         <  6:85685/149‑1 (MQ=255)                    gggCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGtt                                                                                                                                >  4:259258/1‑148 (MQ=255)                      gCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGAc                                                                                                                             <  2:93713/149‑1 (MQ=255)                             tGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACg                                                                                                                      >  6:191606/1‑149 (MQ=255)                              ggAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGc                                                                                                                     >  1:179446/1‑149 (MQ=255)                               gtattgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTCGCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCt                                                                                                                    <  8:157270/145‑1 (MQ=255)                                  gtgtTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTggg                                                                                                                 <  1:264743/149‑1 (MQ=255)                                     ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                              >  8:174368/1‑149 (MQ=255)                                     ttGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCaa                                                                                                              >  2:373773/1‑149 (MQ=255)                                        gTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑g                                                                                                          >  5:348602/1‑149 (MQ=255)                                            aCCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCttt                                                                                                      >  8:276071/1‑149 (MQ=255)                                             ccTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTc                                                                                                     >  5:290821/1‑149 (MQ=255)                                              cTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCt                                                                                                    >  2:222875/1‑149 (MQ=255)                                                ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                         >  5:60084/1‑142 (MQ=255)                                                ggAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑Gc                                                                                                         <  6:60084/142‑1 (MQ=255)                                                 gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTgg                                                                                                  <  6:258928/148‑1 (MQ=255)                                                 gagaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGa                                                                                                 <  8:49708/149‑1 (MQ=255)                                                  agaTCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGAt                                                                                                <  3:259258/149‑1 (MQ=255)                                                     tCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGAc                                                                                             >  7:69598/1‑149 (MQ=255)                                                         ccccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGtt                                                                                         <  8:44867/149‑1 (MQ=255)                                                          cccGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTc                                                                                        >  5:100783/1‑149 (MQ=255)                                                            cGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCgg                                                                                      >  6:380628/1‑149 (MQ=255)                                                             gAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGCACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTGTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCGTTCTGGATGCCTGTTCgg                                                                                      >  6:315739/1‑148 (MQ=255)                                                                  tttCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATtgc                                                                                >  4:45158/1‑149 (MQ=255)                                                                          ccaccaTGCGCTTGCGTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGCGTTGCCCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCg                                                                        >  8:370510/1‑149 (MQ=255)                                                                           caccaTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGATTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCgg                                                                       >  2:317280/1‑149 (MQ=255)                                                                               aTGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGcc                                                                   <  6:100783/149‑1 (MQ=255)                                                                                   gcttgcttCGCATCGAACTCCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTATGCAGCCCGGCTGGTATGCTCAGCCTGACTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCTGCATTGCTGCTGCCGGTGCCAAgc                                                               >  3:273301/1‑149 (MQ=255)                                                                                        cttCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGAt                                                          <  3:45158/149‑1 (MQ=255)                                                                                            gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                       <  7:276071/148‑1 (MQ=255)                                                                                            gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCAt                                                       <  7:174368/148‑1 (MQ=255)                                                                                            gCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATg                                                      <  1:222875/149‑1 (MQ=255)                                                                                             cATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGc                                                     <  1:373773/149‑1 (MQ=255)                                                                                               tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTTCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                   <  6:290821/149‑1 (MQ=255)                                                                                               tCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGccc                                                   <  8:69598/149‑1 (MQ=255)                                                                                                 gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCa                                                  <  2:179446/148‑1 (MQ=255)                                                                                                 gAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAg                                                 <  8:331241/149‑1 (MQ=255)                                                                                                    cTGCACGACACGGTTTCCCAGGCGCAGTTGCTCGCTTTGCAGCTCGGCTTGTCTGCTCAGCCTGAGTTGACCCGAACGCTTGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCTGGTGCCAAGCGCGATCATGCCCAGcac                                              <  4:281345/149‑1 (MQ=255)                                                                                                       cccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                           >  3:347040/1‑149 (MQ=255)                                                                                                        ccGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                           >  6:290548/1‑148 (MQ=255)                                                                                                               ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTAcc                                   <  5:191606/149‑1 (MQ=255)                                                                                                               ggTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTAAACATTGCTTCAGCCGGTGCCAAGCGCGATAATTCCCAGCACCACCAATTAcc                                   >  1:323860/1‑149 (MQ=255)                                                                                                                gTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                  <  1:162701/149‑1 (MQ=255)                                                                                                                     ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                              >  2:196133/1‑148 (MQ=255)                                                                                                                     ccAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTc                              >  6:37502/1‑148 (MQ=255)                                                                                                                      cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                      <  8:284940/139‑1 (MQ=255)                                                                                                                      cAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAAtt                                      >  7:284940/1‑139 (MQ=255)                                                                                                                             aGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACgg                     <  4:347040/149‑1 (MQ=255)                                                                                                                               ttGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATCTTCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGc                    >  1:309389/1‑149 (MQ=255)                                                                                                                                 gCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                 >  5:267289/1‑149 (MQ=255)                                                                                                                                  cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCACCCCACGGCAccc                >  6:293457/1‑149 (MQ=255)                                                                                                                                  cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                 <  5:37502/148‑1 (MQ=255)                                                                                                                                  cccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCAcc                 <  1:374106/148‑1 (MQ=255)                                                                                                                                   ccGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCAACAATTACCTGTTCAACCAACGGCAcc                 >  2:374106/1‑147 (MQ=255)                                                                                                                                    cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                  <  5:146063/129‑1 (MQ=255)                                                                                                                                    cGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCt                                  >  6:146063/1‑129 (MQ=255)                                                                                                                                     gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg               >  2:103345/1‑147 (MQ=255)                                                                                                                                     gCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCg               >  8:337623/1‑147 (MQ=255)                                                                                                                                       tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                               <  6:232367/110‑1 (MQ=255)                                                                                                                                       tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                               <  6:227004/110‑1 (MQ=255)                                                                                                                                       tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                               >  5:227004/4‑113 (MQ=255)                                                                                                                                       tgagCAGACCAGCCGGGCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCAGTGCCAAGCGCGATCATGCCCAGca                                               >  5:232367/4‑113 (MQ=255)                                                                                                                                          gCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGgcg         >  1:111752/1‑148 (MQ=255)                                                                                                                                           cAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCa       <  3:89375/149‑1 (MQ=255)                                                                                                                                             gCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCaga     <  5:290548/149‑1 (MQ=255)                                                                                                                                              cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                           >  1:370033/1‑110 (MQ=255)                                                                                                                                              cccGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGcaccac                                           <  2:370033/110‑1 (MQ=255)                                                                                                                                                cGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAgcg  <  1:196133/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCCGCCACCAGCCCTCCTGGGGCGATGCTGGAAGTGTTGGTGCACCTGGAGATCGACCCCGAGGATTTCCCTACCACCATGCGCTTGCTTCGCATCGAACTGCCCGACACGGTTTCCCAGGCGCAGTTGCCCGCTTTGCAGCCCGGCTGGTCTGCTCAGCCTGAGTTGACCCGAACGCTGGGCAATC‑GCTTTCTGGATGACTGTTCGGCATTGCTGCTGCCGGTGCCAAGCGCGATCATGCCCAGCACCACCAATTACCTGTTCAACCCACGGCACCCGCAGGCGCAGAGCG  >  NC_002947/2780802‑2781093

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