Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,133,981	T→A	6.5%	Q350L (CAG→CTG)	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,133,981	0	T	A	6.5%	180.6 / 2.7	61	Q350L (CAG→CTG)	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base A (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

ACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGTC  >  NC_002947/4133834‑4134129                                                                                                                                                     |                                                                                                                                                     accacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGcc‑tgcctg                                                                                                                                                     >  2:343432/1‑149 (MQ=255)   cacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGc                                                                                                                                                   >  7:177307/1‑149 (MQ=255)    acTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCg                                                                                                                                                  <  3:102168/149‑1 (MQ=255)     cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGt                                                                                                                                                 <  2:157373/149‑1 (MQ=255)        gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTg                                                                                                                                               <  6:176125/148‑1 (MQ=255)         ggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCa                                                                                                                                             <  8:177307/149‑1 (MQ=255)             aaGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTc                                                                                                                                         >  1:317000/1‑149 (MQ=255)              aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCg                                                                                                                                        <  2:317000/149‑1 (MQ=255)                 aTGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGgcg                                                                                                                                     >  6:290290/1‑149 (MQ=255)                       aCGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAAt                                                                                                                               >  3:327719/1‑149 (MQ=255)                           ccTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCa                                                                                                                           <  4:327719/149‑1 (MQ=255)                            cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑AGGCAGGCGTTGCat                                                                                                                                            <  1:80065/131‑2 (MQ=255)                            cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑AGGCAGGCGTTGCat                                                                                                                                            >  2:80065/1‑130 (MQ=255)                               aCCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                                         <  5:290290/147‑1 (MQ=255)                                     tGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                                 >  1:367515/1‑149 (MQ=255)                                      gCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑AGGCAGGCGTTGCat                                                                                                                                            >  4:252146/1‑120 (MQ=255)                                      gCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑AGGCAGGCGTTGCat                                                                                                                                            <  3:252146/121‑2 (MQ=255)                                        cAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGt                                                                                                              >  7:254838/1‑149 (MQ=255)                                          gCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCTAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                            <  1:198671/149‑1 (MQ=255)                                          gCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                            <  4:75413/149‑1 (MQ=255)                                                        tGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGAc                                                                                                           >  2:140456/1‑136 (MQ=255)                                                        tGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGAc                                                                                                           <  1:140456/136‑1 (MQ=255)                                                          cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAg                                                                                              <  8:254838/147‑1 (MQ=255)                                                           gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                           >  4:83882/1‑149 (MQ=255)                                                             cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGa                                                                                         >  1:175261/1‑149 (MQ=255)                                                              aGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGAc                                                                                        <  3:324271/149‑1 (MQ=255)                                                                        gcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                                 <  7:271447/114‑1 (MQ=255)                                                                        gcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTggg                                                                                                                 >  8:271447/1‑114 (MQ=255)                                                                        gcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATc                                                                              <  4:217432/149‑1 (MQ=255)                                                                           aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                           >  1:13809/1‑149 (MQ=255)                                                                           aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                           >  6:378668/1‑149 (MQ=255)                                                                           aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                           >  1:199849/1‑149 (MQ=255)                                                                            accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                          <  3:83882/149‑1 (MQ=255)                                                                             ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                         <  1:148475/149‑1 (MQ=255)                                                                                cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                                       >  4:136774/1‑148 (MQ=255)                                                                                 ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCa                                                                      >  4:136650/1‑148 (MQ=255)                                                                                 ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCa                                                                      >  5:100219/1‑148 (MQ=255)                                                                                  cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    >  7:102014/1‑149 (MQ=255)                                                                                  cTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    <  1:343432/149‑1 (MQ=255)                                                                                     aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAccc                                                                 >  7:111469/1‑149 (MQ=255)                                                                                      cgcgGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    >  6:45873/1‑145 (MQ=255)                                                                                      cgcgGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    <  5:45873/145‑1 (MQ=255)                                                                                       gcgGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCtt                                                               >  7:60833/1‑149 (MQ=255)                                                                                         gggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCg                                                             <  8:102014/149‑1 (MQ=255)                                                                                          ggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                            >  7:95589/1‑149 (MQ=255)                                                                                                  aaGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    >  2:152836/1‑133 (MQ=255)                                                                                                  aaGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                    <  1:152836/133‑1 (MQ=255)                                                                                                   aGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                                    >  5:292773/1‑148 (MQ=255)                                                                                                      ggATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                                <  2:175261/149‑1 (MQ=255)                                                                                                         tGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                            <  5:71279/134‑1 (MQ=255)                                                                                                         tGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                            >  6:71279/1‑134 (MQ=255)                                                                                                               gggCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCt                                       <  8:111469/149‑1 (MQ=255)                                                                                                                  ccTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgc                                    <  4:129530/149‑1 (MQ=255)                                                                                                                   cTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgcc                                   <  2:13809/149‑1 (MQ=255)                                                                                                                      cAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgcc                                <  8:218964/149‑1 (MQ=255)                                                                                                                        gCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                              <  2:307361/149‑1 (MQ=255)                                                                                                                             acCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGTCTCCCGCCATCGGCggg                                                                        <  3:241452/102‑1 (MQ=255)                                                                                                                             acCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                        >  4:241452/1‑102 (MQ=255)                                                                                                                               cgcgCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                       <  6:100219/149‑1 (MQ=255)                                                                                                                                ccgcTCATGC‑ACGCCGGGCCGGGGGTTGCTGTGCGGCGCGCAATTCCCCCCGGCGGGGCGTGACGACCGGCACTACGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTGTTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc                      <  2:367515/148‑1 (MQ=255)                                                                                                                                          aaCGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGtt            >  6:65303/1‑149 (MQ=255)                                                                                                                                          aaCGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGtt            <  7:13409/149‑1 (MQ=255)                                                                                                                                             gcc‑tgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTa         >  4:266625/1‑149 (MQ=255)                                                                                                                                              cc‑tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTAc        <  3:136774/149‑1 (MQ=255)                                                                                                                                                 tgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcg      >  5:241114/1‑149 (MQ=255)                                                                                                                                                   cctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACcgcg    <  6:241114/149‑1 (MQ=255)                                                                                                                                                     tgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGCGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGCGCAAACCCTTCGCCGTATTCGCAGACGTAAGCCTTGCCGCCACTGGCCAGCCACTCCAGTTGTACAGAGTc  >  3:285245/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCC‑TGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGTC  >  NC_002947/4133834‑4134129

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