Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,886,764	C→G	5.5%	W84C (TGG→TGC)	PP_5161 ←	lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,886,764	0	C	G	5.5%	240.4 / 3.8	72	W84C (TGG→TGC)	PP_5161	lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base G (2/2);  total (40/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 = 5.89e-01

CGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCG  >  NC_002947/5886619‑5886910                                                                                                                                                  |                                                                                                                                                   cgcTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                 <  4:90703/149‑1 (MQ=255)   cTGTACATGCCGTGGCCGGGCAGGTCAAAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACATCGTCGCGATTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               <  2:54375/149‑1 (MQ=255)   cTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               <  2:200669/149‑1 (MQ=255)      gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCGTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               <  3:185116/146‑1 (MQ=255)      gACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCGTGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               >  4:185116/1‑146 (MQ=255)       ccAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgcctgc                                                                                                                                           <  7:92938/148‑1 (MQ=255)        cAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               <  1:12636/144‑1 (MQ=255)        cAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAtgc                                                                                                                                               >  2:12636/1‑144 (MQ=255)              gTGGCCGGGCAGGTCACAGCTGCTCACGGCGAAACCCTGCTTGAGCGCCCACTCGATCACATCTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCACGAACAGCGTCGCGCTTGGCTGCTCGGGCAGGCATGGCTGCCCGAGCa                                                                                                                                    <  7:225476/149‑1 (MQ=255)               tGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAg                                                                                                                                   <  7:54029/149‑1 (MQ=255)                           tCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGc                                                                                                                       >  1:330543/1‑149 (MQ=255)                                          gCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGCAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTAcc                                                                                                        >  4:230686/1‑149 (MQ=255)                                                 ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                 <  4:99065/101‑1 (MQ=255)                                                 ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCAt                                                                                                                                                 >  3:99065/1‑101 (MQ=255)                                                 ccTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccagccag                                                                                                 >  8:149567/1‑149 (MQ=255)                                                    gCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtg                                                                                              >  5:315605/1‑149 (MQ=255)                                                     cTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCACATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                             <  4:146638/149‑1 (MQ=255)                                                      ttGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAAcc                                                                                                                           >  1:129142/1‑118 (MQ=255)                                                      ttGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAAcc                                                                                                                           <  2:129142/118‑1 (MQ=255)                                                           tGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCg                                                                                       >  3:77518/1‑149 (MQ=255)                                                               cACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTTCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAg                                                                                   >  3:201988/1‑149 (MQ=255)                                                                aCTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAgg                                                                                  >  7:327311/1‑149 (MQ=255)                                                                  tCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGTAGGtt                                                                                >  3:304590/1‑149 (MQ=255)                                                                    gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                             >  2:257877/1‑134 (MQ=255)                                                                    gATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgt                                                                                             <  1:257877/134‑1 (MQ=255)                                                                      tCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAgg                                                                            <  6:80680/149‑1 (MQ=255)                                                                            gTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCa                                                                         >  3:300586/1‑146 (MQ=255)                                                                             tCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCAGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAg                                                                      >  6:250331/1‑148 (MQ=255)                                                                               ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                           >  2:324378/1‑125 (MQ=255)                                                                               ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCtgtgt                                                                                           <  1:324378/125‑1 (MQ=255)                                                                               ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGTAGGTTCAGGCCATAGaa                                                                    <  4:304590/148‑1 (MQ=255)                                                                               ggTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGaa                                                                    <  4:300586/148‑1 (MQ=255)                                                                                 tAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                             >  5:351559/1‑137 (MQ=255)                                                                                 tAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAg                                                                             <  6:351559/137‑1 (MQ=255)                                                                                  agagGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACccagc                                                                                                    >  2:87565/1‑113 (MQ=255)                                                                                    agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTg                                                              <  2:330543/149‑1 (MQ=255)                                                                                    agGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTg                                                              <  8:327311/149‑1 (MQ=255)                                                                                     ggCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                             <  6:315605/149‑1 (MQ=255)                                                                                      gCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                             >  4:88985/1‑148 (MQ=255)                                                                                           tgtgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTgg                                                             <  8:106232/143‑1 (MQ=255)                                                                                           tgtgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAgaac                                                       >  7:106232/1‑146 (MQ=255)                                                                                             tgGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTTCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCt                                                      <  4:201988/148‑1 (MQ=255)                                                                                                tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTa                                                           <  5:323172/140‑1 (MQ=255)                                                                                                tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTa                                                           >  6:323172/1‑140 (MQ=255)                                                                                                tCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTgcgc                                                  >  8:159314/1‑149 (MQ=255)                                                                                                    aaTAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTtgt                                              >  6:194358/1‑149 (MQ=255)                                                                                                                gcagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggc                                  >  7:69933/1‑149 (MQ=255)                                                                                                                 cagGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                 >  7:203544/1‑149 (MQ=255)                                                                                                                   ggAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggctgg                               >  1:334456/1‑149 (MQ=255)                                                                                                                    gAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGACCCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCAGGCAGGTTCAGGCCATAGAAGAGCTGGTAGTCCTGCGCTTGTGCCGACAGtggatg                                >  6:350787/1‑147 (MQ=255)                                                                                                                       cAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAg                           <  5:194358/149‑1 (MQ=255)                                                                                                                         gCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAgtg                         <  7:159314/149‑1 (MQ=255)                                                                                                                                  ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                 >  6:66264/1‑132 (MQ=255)                                                                                                                                  ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtggct                                 <  5:66264/132‑1 (MQ=255)                                                                                                                                  ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgc                >  5:267439/1‑149 (MQ=255)                                                                                                                                  ttGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgc                >  1:19614/1‑149 (MQ=255)                                                                                                                                    ggCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGgcgc              >  3:135392/1‑149 (MQ=255)                                                                                                                                       tgctCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTACCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAg           >  1:244109/1‑149 (MQ=255)                                                                                                                                        gctCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAgg          >  5:83516/1‑149 (MQ=255)                                                                                                                                         ctCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGc         >  3:244206/1‑149 (MQ=255)                                                                                                                                            gggCAGGCATGGCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCggg      >  8:185055/1‑149 (MQ=255)                                                                                                                                            gggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCgg       <  4:244206/148‑1 (MQ=255)                                                                                                                                             ggCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGt     <  6:83516/149‑1 (MQ=255)                                                                                                                                              gCAGGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                   >  7:111030/1‑118 (MQ=255)                                                                                                                                              gCAGGCATGGCTGCCCGAGCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGtgg                                   <  8:111030/118‑1 (MQ=255)                                                                                                                                              gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGtt                                                           <  5:347364/94‑2 (MQ=255)                                                                                                                                              gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGtt                                                           >  6:347364/1‑93 (MQ=255)                                                                                                                                              gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGAGAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTg    >  5:336872/1‑149 (MQ=255)                                                                                                                                              gCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTg    <  8:203544/149‑1 (MQ=255)                                                                                                                                               cAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgc   <  8:69933/149‑1 (MQ=255)                                                                                                                                                 gCCATTCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGCTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCAGACAGTCGCTGGCGAGTGGTGAACGGCGCCAGGCGCGTgcg  >  4:17629/1‑148 (MQ=255)                                                                                                                                                 gCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTgcg  >  5:127889/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGCTGGACAGGCCGTGGCCGGGCAGGTCACAGCTGATCACGGCGAAACCCTGCTTGAGTGCCCACTCGATCACATGTCGGTAGAGGCCCATGTGGTCGTAATAGCCGTGCAGCAGGAACAGCGTCGCGCTTGGCTGCTCGGGCAGCCATGCCTGCCCGACCAGGTCGAACCCCGCTGCCTGAAAGCTACCCAGCCAGCTGTGTGCCGGCAGGTTCAGGCCATAGAAGCGCTGGTAGTCCTGCGCTTGTGCCGACAGTGGCTGGCGAGTGGTGAACGGCGCCAGGCGGGTGCG  >  NC_002947/5886619‑5886910

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