Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,102,518	T→G	6.5%	Y117D (TAC→GAC)	PP_2712 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,102,518	0	T	G	6.5%	196.0 / 4.6	62	Y117D (TAC→GAC)	PP_2712	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base G (2/2);  total (30/32)
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.65e-01

CTGAACCCTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGATCATC  >  NC_002947/3102372‑3102662                                                                                                                                                    |                                                                                                                                                 cTGAACCCTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTa                                                                                                                                                 >  2:55510/1‑148 (MQ=255)    aaCCCTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTaca                                                                                                                                               <  6:43345/147‑1 (MQ=255)      cccTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAg                                                                                                                                           >  2:56373/1‑149 (MQ=255)      cccTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAg                                                                                                                                           >  4:359718/1‑149 (MQ=255)           gATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCt                                                                                                                                        <  4:142899/147‑1 (MQ=255)           gATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTg                                                                                                                                       <  8:324934/148‑1 (MQ=255)              cTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCa                                                                                                                                   >  6:33461/1‑149 (MQ=255)                      gCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTAt                                                                                                                           <  3:359718/149‑1 (MQ=255)                        atcatcAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGa                                                                                                                         <  5:207685/149‑1 (MQ=255)                        atcatcAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGa                                                                                                                         >  5:229987/1‑149 (MQ=255)                         tcatcaACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGAc                                                                                                                        <  3:230159/149‑1 (MQ=255)                           atcaACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGAc                                                                                                                        >  3:92791/1‑147 (MQ=255)                            tcaACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAg                                                                                                                     <  1:55510/149‑1 (MQ=255)                              aaCAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTgg                                                                                                                                      >  6:238056/1‑130 (MQ=255)                              aaCAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTgg                                                                                                                                      <  5:238056/130‑1 (MQ=255)                              aaCAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGt                                                                                                                   <  8:275062/149‑1 (MQ=255)                              aaCAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGt                                                                                                                   >  1:337279/1‑149 (MQ=255)                               aCAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACCTCCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGt                                                                                                                   <  1:287003/149‑1 (MQ=255)                                 agcagTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACAAGCCTGGGCACGCTGTATGACCAGGTGc                                                                                                                 >  7:107115/1‑148 (MQ=255)                                                  gACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGtt                                                                                               <  7:355913/149‑1 (MQ=255)                                                   aCGGGTGGCCTCGGTGCACTTGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGttt                                                                                              <  5:33461/149‑1 (MQ=255)                                                        tGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGt                                                                                          >  6:47456/1‑148 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGc                                                                                                        <  2:298250/130‑1 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGc                                                                                                        >  1:298250/1‑130 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGc                                                                                     >  8:207562/1‑149 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTGACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGAt                                                                                       <  7:89113/147‑1 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTGACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGAt                                                                                       <  7:94173/147‑1 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTGACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGAt                                                                                       >  8:89113/1‑147 (MQ=255)                                                            cTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTGACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGAt                                                                                       >  8:94173/1‑147 (MQ=255)                                                                gTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCgag                                                                                 <  6:230679/149‑1 (MQ=255)                                                                 tGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCgaga                                                                                <  2:337279/149‑1 (MQ=255)                                                                   cACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGAt                                                                               <  8:107115/148‑1 (MQ=255)                                                                        ggCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGcc                                                                           <  7:207562/147‑1 (MQ=255)                                                                         gCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGa                                                                         >  1:75552/1‑148 (MQ=255)                                                                          cgaGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGAc                                                                        >  6:213883/3‑148 (MQ=255)                                                                                cAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGAtgc                                                                 >  3:121930/1‑149 (MQ=255)                                                                                 aaGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGAtgct                                                                >  5:282282/1‑149 (MQ=255)                                                                                  aGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGAtgctg                                                               <  4:92791/149‑1 (MQ=255)                                                                                          aaaGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGc                                                       >  7:207675/1‑149 (MQ=255)                                                                                            aGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCt                                                     <  2:75552/149‑1 (MQ=255)                                                                                             gCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCg                                                                                                  >  6:185524/1‑103 (MQ=255)                                                                                             gCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCg                                                                                                  <  5:185524/103‑1 (MQ=255)                                                                                                 gCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGc                                                  <  8:207675/147‑1 (MQ=255)                                                                                                 gCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGc                                                  <  4:121930/147‑1 (MQ=255)                                                                                                   cTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGAt                                              <  5:213883/149‑1 (MQ=255)                                                                                                      gAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATc                                             <  5:47456/147‑1 (MQ=255)                                                                                                              cATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGa                                               <  6:218994/137‑1 (MQ=255)                                                                                                              cATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGa                                               >  5:218994/1‑137 (MQ=255)                                                                                                              cATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCAGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACg                                   >  6:256892/1‑149 (MQ=255)                                                                                                               aTGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTg                                                                                                          <  6:370475/77‑1 (MQ=255)                                                                                                               aTGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTg                                                                                                          >  5:370475/1‑77 (MQ=255)                                                                                                                    aTTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATc                              >  1:362094/1‑148 (MQ=255)                                                                                                                        gACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAg                         >  2:212128/1‑149 (MQ=255)                                                                                                                          cc‑ccGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGc                        >  6:106100/1‑148 (MQ=255)                                                                                                                                      tGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGc                                                  <  7:262949/111‑1 (MQ=255)                                                                                                                                      tGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGc                                                  >  8:262949/1‑111 (MQ=255)                                                                                                                                       gCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTggg           >  2:181402/1‑149 (MQ=255)                                                                                                                                        cAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTgg            <  1:56373/147‑1 (MQ=255)                                                                                                                                             ttCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGat      >  5:206261/1‑148 (MQ=255)                                                                                                                                              tCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGatca    >  5:264121/1‑149 (MQ=255)                                                                                                                                                ggCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGatcatc  <  6:229987/149‑1 (MQ=255)                                                                                                                                                ggCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGatcatc  <  1:212128/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CTGAACCCTTGATCTGGGAGGGCATCATCAACAGCAGTTGGGCGGGGGCGACGGGTGGCCTCGGTGCACTGGGCAAGGTCAAGGCGGCGAAAGCTAGCCTGGAGCAGGCCATGAAATTGGACC‑CCGCCGCCCTGCAAGGTTCGGCTTACACCAGCCTGGGCACGCTGTATGACCAGGTGCCAGGCTGGCCTATCGGTTTCGGTGATGCCGAGATGGCCGACAAGATGCTGCGCAAAGCCTTGCAGATCAACCCCAACGGTATCGACAGCAACTATTTCTGGGCCGATCATC  >  NC_002947/3102372‑3102662

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