Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,746	C→T	8.3%	A79V (GCT→GTT) ‡	PP_4562 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,746	0	C	T	8.3%	172.5 / 4.3	60	A79V (GCT→GTT) ‡	PP_4562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base T (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

TGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAG  >  NC_002947/5182602‑5182892                                                                                                                                                 |                                                                                                                                                   tGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGc                                                                                                                                                    <  6:166526/145‑1 (MQ=255) tGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGAAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGc                                                                                                                                                    >  5:166526/1‑145 (MQ=255)   gACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCAc                                                                                                                                                <  1:61135/147‑1 (MQ=255)      aacaaACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCtt                                                                                                                                             >  2:178336/1‑147 (MQ=255)      aacaaACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGa                                                                                                                                           <  8:43788/149‑1 (MQ=255)           acgcccccTGTGGTGGTTCGCAGGGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGTGAGCCTTGAGGaaa                                                                                                                                      <  8:238260/143‑1 (MQ=255)            cGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAAt                                                                                                                                     <  4:266762/149‑1 (MQ=255)                 ccTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGTGAGCCTTGAGGaaa                                                                                                                                      >  1:179416/1‑143 (MQ=255)                 ccTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGTGAGCCTTGAGGaaa                                                                                                                                      <  2:179416/143‑1 (MQ=255)                     tggtggTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTgcg                                                                                                                            <  4:172926/149‑1 (MQ=255)                     tggtggTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTg                                                                                                                              <  5:131757/147‑1 (MQ=255)                             cGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTg                                                                                                                    >  5:108687/1‑149 (MQ=255)                                aGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGc                                                                                                                  >  1:149780/1‑148 (MQ=255)                                  tGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGTGAGCCTTGAGGaaa                                                                                                                                      >  4:134550/1‑126 (MQ=255)                                  tGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGTGAGCCTTGAGGaaa                                                                                                                                      <  3:134550/126‑1 (MQ=255)                                       gTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGcgccg                                                                                                          >  3:138606/1‑149 (MQ=255)                                        tGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGcgccgc                                                                                                         <  1:178336/149‑1 (MQ=255)                                             ctgctgGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTccac                                                                                                    <  8:144107/149‑1 (MQ=255)                                              tgctgGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTccac                                                                                                    >  5:357662/1‑148 (MQ=255)                                                ctgGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGCTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGACTGGCAAGCGCCGCTccacca                                                                                                  >  8:1508/1‑148 (MQ=255)                                                   ggCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCa                                                                                               <  4:138606/148‑1 (MQ=255)                                                    gCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGTAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGc                                                                                             >  1:344296/1‑149 (MQ=255)                                                     cAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGcc                                                                                            >  4:53713/1‑149 (MQ=255)                                                           ccAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCt                                                                                           <  5:364165/144‑1 (MQ=255)                                                           ccAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCt                                                                                           >  6:364165/1‑144 (MQ=255)                                                           ccAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGccc                                                                                      >  4:285064/1‑149 (MQ=255)                                                                  gACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAAt                                                                                                                                     >  6:208021/1‑95 (MQ=255)                                                                  gACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAAt                                                                                                                                     <  5:208021/95‑1 (MQ=255)                                                                    cTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCg                                                                             <  3:285064/149‑1 (MQ=255)                                                                         tCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAGGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCgg                                                                        >  3:114772/1‑149 (MQ=255)                                                                         tCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCgg                                                                        >  4:33877/1‑149 (MQ=255)                                                                           aaCCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGa                                                                      <  3:33877/149‑1 (MQ=255)                                                                            aCCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGAc                                                                     <  5:364867/149‑1 (MQ=255)                                                                              ctgctgctgCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATgg                                                                                                                                   >  1:246285/1‑85 (MQ=255)                                                                              ctgctgctgCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATgg                                                                                                                                   <  2:246285/85‑1 (MQ=255)                                                                                      gcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCtgctgc                                                           >  5:90426/1‑149 (MQ=255)                                                                                       cacaCGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCg                                                                         >  6:70406/1‑134 (MQ=255)                                                                                       cacaCGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCg                                                                         <  5:70406/134‑1 (MQ=255)                                                                                                gAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCg                                                 >  3:255204/1‑149 (MQ=255)                                                                                                        tGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTg                                         >  1:180527/1‑149 (MQ=255)                                                                                                          aGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGTCTCGCTGCTGCGCGACACCCGTCCGGGTGGc                                       <  7:1508/149‑1 (MQ=255)                                                                                                              gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTg                                                                                                                    <  8:196440/68‑1 (MQ=255)                                                                                                              gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTg                                                                                                                    >  7:196440/1‑68 (MQ=255)                                                                                                               cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGGACGGGTGTCGcgca                                  <  3:53713/149‑1 (MQ=255)                                                                                                                             gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAc                           >  5:257434/1‑142 (MQ=255)                                                                                                                             gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAc                           <  6:257434/142‑1 (MQ=255)                                                                                                                             gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTcc                       <  1:86999/146‑1 (MQ=255)                                                                                                                             gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTcc                       >  2:86999/1‑146 (MQ=255)                                                                                                                             gCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGc                     >  6:265566/1‑148 (MQ=255)                                                                                                                               ccaccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAAc                  >  2:300612/1‑149 (MQ=255)                                                                                                                                caccaATTTCCTCAGGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACt                 <  4:114772/149‑1 (MQ=255)                                                                                                                                caccaATTTCCTCAAGGCTCAGCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCACCt                 >  3:348475/1‑149 (MQ=255)                                                                                                                                caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACt                 >  3:233625/1‑149 (MQ=255)                                                                                                                                caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACt                 >  7:63368/1‑149 (MQ=255)                                                                                                                                caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACt                 >  6:70185/1‑149 (MQ=255)                                                                                                                                  ccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTc               >  4:284942/1‑149 (MQ=255)                                                                                                                                          cTCAAGGCTCACCTTGAGTAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGc       <  2:344296/149‑1 (MQ=255)                                                                                                                                          cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACa          >  7:309985/1‑146 (MQ=255)                                                                                                                                          cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACa          <  8:309985/146‑1 (MQ=255)                                                                                                                                          cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGc       <  2:180527/149‑1 (MQ=255)                                                                                                                                              cggCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCACGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTaa   <  4:255204/148‑1 (MQ=255)                                                                                                                                               ggCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg  <  6:90426/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   TGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAG  >  NC_002947/5182602‑5182892

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