Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	904,934	Δ1 bp	6.6%	coding (1193/2160 nt)	PP_0788 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	904,934	0	C	.	6.6%	192.1 / 6.6	61	coding (1193/2160 nt)	PP_0788	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base . (2/2);  total (31/30)
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.04e-01

GTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTG  >  NC_002947/904788‑905076                                                                                                                                                   |                                                                                                                                                gTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCg                                                                                                                                               >  6:296700/1‑149 (MQ=255)  ttCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCgg                                                                                                                                              >  6:111406/1‑149 (MQ=255)     aTCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGgcgc                                                                                                                                           >  5:342969/1‑149 (MQ=255)         acaACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCc‑gcg                                                                                                                                      >  2:25043/1‑149 (MQ=255)              gacgacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTgg                                                                                                                                 <  1:209408/149‑1 (MQ=255)              gacgacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTgg                                                                                                                                 >  2:209408/1‑149 (MQ=255)                cgacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCG‑CGGCGCCGGCGAtttcgt                                                                                                                               <  2:56548/149‑5 (MQ=255)                 gacAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTg                                                                                                                              <  5:296700/149‑1 (MQ=255)                  acAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTgg                                                                                                                             <  5:111406/149‑1 (MQ=255)                    aaCCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGTACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGGTCAACGAACTCTCCGGCGCC‑GCGATTGGCTGGGa                                                                                                                           >  8:245189/1‑149 (MQ=255)                    aaCCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTgg                                                                                                                             >  5:23735/1‑147 (MQ=255)                    aaCCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGa                                                                                                                           <  2:259152/149‑1 (MQ=255)                      ccTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTgg                                                                                                                                 <  5:56362/141‑1 (MQ=255)                      ccTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTgg                                                                                                                                 >  6:56362/1‑141 (MQ=255)                          ccGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGc                                                                                                                       <  3:150178/147‑1 (MQ=255)                             gCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAg                                                                                                                  >  8:48379/1‑149 (MQ=255)                             gCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAg                                                                                                                  <  7:48379/149‑1 (MQ=255)                                        ggCGCTGTTTACGCGCCTTCTTATTCCCCTGCTCAATGTGGCATTTCTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGGGCGCCGTTTGCTCAACGAAATCGCCGTCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGGCAGc                                                                                                       <  8:320190/149‑1 (MQ=255)                                         gcgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCg                                                                                                      >  2:306465/1‑149 (MQ=255)                                          cgcTGATTACTCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGTCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGt                                                                                                     <  6:102484/149‑1 (MQ=255)                                          cgcTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGt                                                                                                     <  6:342969/149‑1 (MQ=255)                                              gATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTAAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACa                                                                                                 >  4:215605/1‑149 (MQ=255)                                               aTTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTgg                                                                                                                             <  3:219134/120‑1 (MQ=255)                                               aTTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTAAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTgg                                                                                                                             >  4:219134/1‑120 (MQ=255)                                                ttACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGc                                                                                               >  1:251703/1‑149 (MQ=255)                                                    gcgcCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                           <  2:251703/149‑1 (MQ=255)                                                    gcgcCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGc                                                                                           <  6:23735/149‑1 (MQ=255)                                                        cTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGa                                                                                                   >  4:104524/1‑137 (MQ=255)                                                        cTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGa                                                                                                   <  3:104524/137‑1 (MQ=255)                                                          tCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGc                                                                                     >  2:338791/1‑149 (MQ=255)                                                          tCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGc                                                                                     >  5:359072/1‑149 (MQ=255)                                                                   ccTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGc                                                                            <  6:359072/149‑1 (MQ=255)                                                                    cTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCg                                                                           >  3:35572/1‑149 (MQ=255)                                                                    cTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCg                                                                           >  6:193624/1‑149 (MQ=255)                                                                        tCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATc                                                                       >  3:67049/1‑149 (MQ=255)                                                                          aaGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATc                                                                       <  5:193624/147‑1 (MQ=255)                                                                            ggTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGc                                                                                               >  6:18552/1‑121 (MQ=255)                                                                            ggTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGc                                                                                               <  5:18552/121‑1 (MQ=255)                                                                             gTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCa                                                                  <  4:35572/149‑1 (MQ=255)                                                                                   ttACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCtt                                                            <  1:338791/149‑1 (MQ=255)                                                                                       tGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGcg                                                         >  1:235237/1‑147 (MQ=255)                                                                                          aCAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTTTCGTCCAGCGCTTGCTCAAt                                                     <  3:215605/149‑1 (MQ=255)                                                                                            aaGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGa                                                   >  7:196483/1‑149 (MQ=255)                                                                                                         gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGAt                                      >  5:142248/1‑149 (MQ=255)                                                                                                         gCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGAt                                      <  8:196483/149‑1 (MQ=255)                                                                                                                  gCCACCCGGCGCGGCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCgg                                                                                   <  2:156737/95‑1 (MQ=255)                                                                                                                  gCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCgg                                                                                   >  1:156737/1‑95 (MQ=255)                                                                                                                     aCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCcga                                          >  6:321974/1‑133 (MQ=255)                                                                                                                     aCCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCcga                                          <  5:321974/133‑1 (MQ=255)                                                                                                                       ccGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGc                        >  1:240299/1‑149 (MQ=255)                                                                                                                        cGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGcc                       <  7:245189/149‑1 (MQ=255)                                                                                                                                tgTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCGGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAtttt                                               <  4:184939/116‑1 (MQ=255)                                                                                                                                     gCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgac           <  1:306465/148‑1 (MQ=255)                                                                                                                                       tCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGAGTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTcgacgcc        >  7:264275/1‑147 (MQ=255)                                                                                                                                         cagccaaTCG‑CGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGa                     <  2:149653/129‑1 (MQ=255)                                                                                                                                         cagccaaTCG‑CGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGa                     >  1:149653/6‑134 (MQ=255)                                                                                                                                         aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGc                                                                            >  1:45858/1‑79 (MQ=255)                                                                                                                                         aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGc                                                                            <  2:45858/79‑1 (MQ=255)                                                                                                                                         aaCGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGAt                                      >  2:173348/1‑117 (MQ=255)                                                                                                                                            ccaaTCG‑CGGCGCCGGCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCt   >  7:262607/3‑149 (MQ=255)                                                                                                                                             aaaTCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTg  <  6:142248/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GTTCATCCACAACGACGACAACCTGCCGGCCAGCCTGCGGGCGCTGATTACGCGCCTTCATATTCCCCTGCTCAAGGTGGCATTACTGGACAAGGGGCTGTTCAGCCGGAGTAGCCACCCGGCGCGCCGTTTGCTCAACGAAATCGCCGGCGCC‑GCGATTGGCTGGGAGTGCGGCAGCGAGTGCCAGCGTGACAGCTTGCACATGCGGGTAGAGCGTATCGTCCAGCGCTTGCTCAATGATTTTGCCGACGATGCCAGCCTGTTCGCCGAGCTGCTCGACGACTTCCTG  >  NC_002947/904788‑905076

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