Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,887,077	A→G	6.8%	intergenic (‑62/‑68)	PP_5161 ← / → PP_5162	lysophospholipase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,887,077	0	A	G	6.8%	173.4 / 2.7	59	intergenic (‑62/‑68)	PP_5161/PP_5162	lysophospholipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGCCTACA  >  NC_002947/5886929‑5887223                                                                                                                                                     |                                                                                                                                                   ggCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCca                                                                                                                                                    <  3:61882/149‑1 (MQ=255) ggCAGGGGGCAGCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCca                                                                                                                                                    <  7:144611/149‑1 (MQ=255)   cAGGGGGCAGCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCcaca                                                                                                                                                  <  5:219654/149‑1 (MQ=255)    aGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAAGCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAg                                                                                                                                                 <  5:179298/149‑1 (MQ=255)     gggggCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAgg                                                                                                                                                <  5:255080/149‑1 (MQ=255)          cATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCg                                                                                                                                           <  8:224923/149‑1 (MQ=255)                gTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTCCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCa                                                                                                                                     <  3:69998/149‑1 (MQ=255)                     tCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCttt                                                                                                                                >  6:142883/1‑149 (MQ=255)                        acaTTGATGCTAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAg                                                                                                                             <  5:77926/149‑1 (MQ=255)                           ttGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                                          <  2:114707/149‑1 (MQ=255)                              aTGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAtt                                                                                                                        >  3:21880/1‑148 (MQ=255)                              aTGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAtt                                                                                                                        <  6:4854/148‑1 (MQ=255)                                gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                                          <  8:236362/144‑1 (MQ=255)                                gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGaaa                                                                                                                          >  7:236362/1‑144 (MQ=255)                                gCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAg                                                                                                                     >  1:245443/1‑149 (MQ=255)                                  aaGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGc                                                                                                                    <  3:26695/148‑1 (MQ=255)                                  aaGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACAAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCTGCTTTTAGAAATTGAGc                                                                                                                    >  1:137416/1‑148 (MQ=255)                                   aGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCaa                                                                                                                  <  4:21880/149‑1 (MQ=255)                                    gTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGc                                                                                                                    >  5:258003/1‑146 (MQ=255)                                      aTTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGa                                                                                                                >  6:148944/1‑148 (MQ=255)                                        tGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                                         >  1:68909/1‑121 (MQ=255)                                        tGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGgc                                                                                                                                         <  2:68909/121‑1 (MQ=255)                                         gATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAg                                                                                                             <  1:160121/148‑1 (MQ=255)                                             tGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCt                                                                                                        >  3:116451/1‑149 (MQ=255)                                              ggAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCtt                                                                                                       >  4:142039/1‑149 (MQ=255)                                               gAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAAGTGAGCAAGACAGTTGCtgc                                                                                                      >  4:98421/1‑147 (MQ=255)                                                     tCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAg                                                                                                >  4:96567/1‑149 (MQ=255)                                                       aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                                               <  3:142039/148‑1 (MQ=255)                                                       aGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAgg                                                                                               <  2:245443/148‑1 (MQ=255)                                                            tGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAcc                                                                                         >  1:11160/1‑149 (MQ=255)                                                            tGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCAcc                                                                                         >  3:203515/1‑149 (MQ=255)                                                                   cgcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCtt                                                                                  <  3:32767/149‑1 (MQ=255)                                                                    gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTc                                                                                 <  5:142883/149‑1 (MQ=255)                                                                    gcCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTc                                                                                 <  5:148944/149‑1 (MQ=255)                                                                                           gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGa                                                                                                                            >  1:270741/1‑83 (MQ=255)                                                                                           gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGa                                                                                                                            <  2:270741/83‑1 (MQ=255)                                                                                           gAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                                           <  4:116451/148‑1 (MQ=255)                                                                                            aacaacCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                                           <  1:115675/147‑1 (MQ=255)                                                                                            aacaacCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAAc                                                           >  2:115675/1‑147 (MQ=255)                                                                                                  cAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgtg                                                    >  6:257217/1‑148 (MQ=255)                                                                                                        cGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTATTAGAAATTGCGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGATTAGCCGAATCTTCAAACCGCTGTGCACGGt                                              >  8:142189/1‑148 (MQ=255)                                                                                                          aaCTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                                      >  8:142072/1‑90 (MQ=255)                                                                                                          aaCTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTc                                                                                                      <  7:142072/90‑1 (MQ=255)                                                                                                                    gggAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg                                 <  1:113950/149‑1 (MQ=255)                                                                                                                      gAGCGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg                                 <  5:47984/147‑1 (MQ=38)                                                                                                                      gAGCGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg                                 >  6:47984/1‑147 (MQ=38)                                                                                                                         cGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCgg                            >  3:174577/1‑149 (MQ=255)                                                                                                                         cGGCCTTGCGTCGCGAAAGGCCGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCgg                            >  3:51230/1‑149 (MQ=38)                                                                                                                          ggCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGgc                           >  7:230467/1‑149 (MQ=255)                                                                                                                             cTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGcc                        >  4:103779/1‑149 (MQ=255)                                                                                                                                   cgACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATc                  <  4:174577/149‑1 (MQ=255)                                                                                                                                    gACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGc                                                                        <  7:116126/94‑1 (MQ=255)                                                                                                                                    gACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGc                                                                        >  8:116126/1‑94 (MQ=255)                                                                                                                                       gCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg                                 <  1:5497/130‑1 (MQ=255)                                                                                                                                       gCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCttgttg                                 >  2:5497/1‑130 (MQ=255)                                                                                                                                         aaGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTgg            >  1:170150/1‑149 (MQ=255)                                                                                                                                           ggCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTgg            >  6:41020/1‑147 (MQ=255)                                                                                                                                             ccGCTCCCGCAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAgg        <  4:51230/149‑1 (MQ=255)                                                                                                                                              cGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGc       <  5:257217/149‑1 (MQ=255)                                                                                                                                                   ccACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGCCTa    >  1:241256/1‑147 (MQ=255)                                                                                                                                                   ccACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGCCTAc   <  2:170150/148‑1 (MQ=255)                                                                                                                                                   ccACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGCCTACa  <  5:41020/149‑1 (MQ=255)                                                                                                                                                     acaGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgt                                                     >  3:127475/1‑96 (MQ=255)                                                                                                                                                     acaGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCtgt                                                     <  4:127475/96‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCAGGGGGCATCGGGTACATCCACATTGATGCAAGTATTGATCTGGAATCTTCAGCTCTGGCCGGCGCCATGGCAAGCTTGCTCACCTCGAACAACCAGCATCGAACTAATTGCGGGAGCGGCCTTTCGCGACGCAAGGCCGCTCCCACAGGGATCGGCGCCAGCTTTTAGAAATTGAGCAAGACAGTTGCTTCCACAAGGGCCACCGACTCTTCATGAATAGCCGAATCATCAAACCGCTGTGCACGGTCATCGCCTTGTTGGTCGGCGCCGTAATCCTCTGGCAGGCCTACA  >  NC_002947/5886929‑5887223

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