Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,240,567	C→G	8.1%	H43Q (CAC→CAG)	PP_4618 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,240,567	0	C	G	8.1%	194.9 / 6.6	62	H43Q (CAC→CAG)	PP_4618	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base G (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

TGTGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGG  >  NC_002947/5240423‑5240712                                                                                                                                                 |                                                                                                                                                  tgtgtgGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGAcc                                                                                                                                               <  5:187738/149‑1 (MQ=255)     tgGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCGCTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCa                                                                                                                                           >  8:226448/1‑149 (MQ=255)           tACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCCGCGCTGCAGGTCGTGCAgcgca                                                                                                                                      <  2:255947/148‑2 (MQ=255)             ccccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                       >  8:300436/1‑145 (MQ=255)               ccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                  >  5:285232/1‑148 (MQ=255)                cATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                  >  6:174197/1‑147 (MQ=255)                cATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATaa                                                                                                                                <  4:26958/149‑1 (MQ=255)                  tGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                              <  1:145482/149‑1 (MQ=255)                   gAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgcg                                                                                                                             >  3:127868/1‑149 (MQ=255)                      gCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCgag                                                                                                                                    <  7:309297/139‑1 (MQ=255)                      gCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCgag                                                                                                                                    >  8:309297/1‑139 (MQ=255)                        cACGTCAAGCCTGCCTGGCCTGCCTGAAACGCGAGCCTGTCGACATCCTGGAAGCCGCACTGTGGATTCCTCCCGAGCACGACCCTGTTGTCCACCCCCCTGCCAGCCTTCATCCGCTGCACGACCTCCAGCGCGAGATAACCGCCAAc                                                                                                                        >  8:130039/1‑149 (MQ=255)                         aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAcc                                                                                                                       <  4:127868/149‑1 (MQ=255)                                gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCTCTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                <  7:226448/149‑1 (MQ=255)                                gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                <  6:285232/149‑1 (MQ=255)                                   tgcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCt                                                                                                             >  6:43627/1‑149 (MQ=255)                                    gcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCACCGCGAGATAAGCGCCAACCTGCCGATGCtt                                                                                                            >  3:65999/1‑149 (MQ=255)                                     cctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCt                                                                                                             <  2:133749/147‑1 (MQ=255)                                     cctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                           <  5:174197/149‑1 (MQ=255)                                     cctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                           <  7:300436/149‑1 (MQ=255)                                      ctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                           >  7:81044/1‑148 (MQ=255)                                         gcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCAc                                                                                                        <  5:119948/148‑1 (MQ=255)                                          cctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCaa                                                                                                                         <  1:149916/130‑1 (MQ=255)                                          cctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCaa                                                                                                                         >  2:149916/1‑130 (MQ=255)                                            tgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCt                                                                                                             >  5:68927/1‑140 (MQ=255)                                            tgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCt                                                                                                             <  6:68927/140‑1 (MQ=255)                                              cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATg                                                                                                    <  5:9949/147‑1 (MQ=255)                                                tgtaaCGCGAGCCTGTCGACCTGCTGGAAGCCGCAGTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGaa                                                                                                <  7:130039/146‑1 (MQ=255)                                                tggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                <  8:322070/133‑1 (MQ=255)                                                tggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                >  7:322070/1‑133 (MQ=255)                                                  gAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGaa                                                                                                <  2:96734/147‑1 (MQ=255)                                                      gcgAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgc                                                                                                                                       >  8:108185/1‑104 (MQ=255)                                                      gcgAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgc                                                                                                                                       <  7:108185/104‑1 (MQ=255)                                                                  aCCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgctg                                                                              >  5:264766/1‑149 (MQ=255)                                                                   ccTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACtgctgc                                                                             <  5:43627/149‑1 (MQ=255)                                                                    ctgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcg                                                                            >  7:34669/1‑149 (MQ=255)                                                                      gctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATg                                                                                                               <  8:137433/112‑1 (MQ=255)                                                                      gctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATg                                                                                                               >  7:137433/1‑112 (MQ=255)                                                                       ctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCa                                                                         <  8:81044/149‑1 (MQ=255)                                                                          gAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGccc                                                                                         >  2:24622/1‑130 (MQ=255)                                                                          gAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGccc                                                                                         <  1:24622/130‑1 (MQ=255)                                                                               cGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATg                                                                 <  6:264766/149‑1 (MQ=255)                                                                                             ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTca                                                    <  7:258281/148‑1 (MQ=255)                                                                                                 cGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCa                                                                                        >  7:91618/1‑108 (MQ=255)                                                                                                 cGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCa                                                                                        <  8:91618/108‑1 (MQ=255)                                                                                                       cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGaa                                          >  7:174585/1‑148 (MQ=255)                                                                                                             tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGc                                                         <  5:233726/127‑1 (MQ=255)                                                                                                             tGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGc                                                         >  6:233726/1‑127 (MQ=255)                                                                                                               gtgtCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAcc                                   <  8:34669/147‑1 (MQ=255)                                                                                                                    gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                           <  5:229202/118‑1 (MQ=255)                                                                                                                    gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                           >  6:229202/1‑118 (MQ=255)                                                                                                                       cccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                         >  7:21400/1‑149 (MQ=255)                                                                                                                        ccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGcccc                        >  6:49191/1‑149 (MQ=255)                                                                                                                                aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcggc                >  2:69367/1‑149 (MQ=255)                                                                                                                                aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcggc                >  1:199408/1‑149 (MQ=255)                                                                                                                                  ccTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga              <  4:83785/149‑1 (MQ=255)                                                                                                                                       tcgcgctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCAAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg         >  1:323937/3‑149 (MQ=255)                                                                                                                                       ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg         <  4:65999/149‑1 (MQ=255)                                                                                                                                        tgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATgg        <  8:21400/149‑1 (MQ=255)                                                                                                                                           gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGAc      >  6:87066/1‑148 (MQ=255)                                                                                                                                               cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg  <  5:49191/148‑1 (MQ=255)                                                                                                                                               cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg  <  2:199408/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TGTGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGG  >  NC_002947/5240423‑5240712

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