Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,525,497	C→T	6.2%	A159A (GCC→GCT)	fadBA →	type II 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,525,497	0	C	T	6.2%	181.0 / 2.1	64	A159A (GCC→GCT)	fadBA	type II 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATCAACGTCAACCTGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTACGTGCTTC  >  NC_002947/2525351‑2525646                                                                                                                                                   |                                                                                                                                                       atcaACGTCAACCTGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCAATATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGgccgcc                                                                                                                                                     <  1:814767/150‑1 (MQ=255)  tcaACGTCAACCTGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCt                                                                                                                                                    >  1:767590/1‑150 (MQ=255)           aCCTGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCt                                                                                                                                          <  2:1701050/150‑1 (MQ=255)            ccTGATCGGCAGCTTCAGCCGGTTGCGTCTGGCTGCGGCGTGCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCACCAACACGGCCTCCAGTGCCGCCTACGACCGCCAGATGGGCCAGGCCGGCGTcg‑ccggccc                                                                                                                                         >  2:213571/1‑146 (MQ=255)             cTGATCGGCAGCTTCAACCTGTTGAGTCAGGGTTCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTcc                                                                                                                                        >  1:988939/1‑150 (MQ=255)              tGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCa                                                                                                                                       >  2:1277026/1‑150 (MQ=255)               gATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCaa                                                                                                                                      <  1:82944/150‑1 (MQ=255)                   ggCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCTGCCTACG‑CCGCCTCCAAGGGt                                                                                                                                  >  2:1599606/1‑150 (MQ=255)                   ggCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCTGCCTACG‑CCGCCTCCAAGGGt                                                                                                                                  >  2:1599269/1‑150 (MQ=255)                   ggCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGt                                                                                                                                  <  1:1388076/150‑1 (MQ=255)                     cAGCTTCAACCTGTTGCGTGTGGCTGCGGCGGACATGGCCGAAGGGGCTGCCGATTAGAGCGGCGAGCGTGGGGTCATCATCAACTCGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGc                                                                                                                                <  2:426136/150‑1 (MQ=255)                       gCTTCAACCTGTTGCGTCTGGCTGCGGCGGACATGGCCGAAGTGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCa                                                                                                                              <  2:780536/150‑1 (MQ=255)                        cTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTGCG‑CCGCCTCCAAGGGTGCCAt                                                                                                                             >  2:1086699/1‑150 (MQ=255)                                    gCGTCTGGCTGCGGCGGCCATGGCCAAAGGGGCTGCCGATGATAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGAc                                                                                                                 >  2:1263213/1‑150 (MQ=255)                                          ggCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTgcc                                                                                                           <  1:899163/150‑1 (MQ=255)                                              gcggcgGCCATGGCCGAAGGGGCTGCCGATGAGGGCGGCGAGCGTGGGGTCATCATCAACACGGAGCTCAACGTGGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTgccggcc                                                                                                       <  1:329688/150‑1 (MQ=255)                                                     ccATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTAGCCAGGCCGCCTACG‑CCGCCTCCAAGGCTGCCATTGCCAGCCTGACCTTGCCAGCcgcgcgcg                                                                                                >  2:244461/1‑150 (MQ=255)                                                       aTGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGaa                                                                                              >  1:817105/1‑150 (MQ=255)                                                            cGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGc                                                                                         >  2:744644/1‑150 (MQ=255)                                                            cGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATAGGCCAGGCCGCCTACG‑CCCACATGAAATGTCCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGc                                                                                         <  1:330869/150‑1 (MQ=255)                                                             gAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCa                                                                                        >  2:1912766/1‑150 (MQ=255)                                                                 gggCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCt                                                                                    >  1:1865221/1‑150 (MQ=255)                                                                      gCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGc                                                                               >  2:331839/1‑150 (MQ=255)                                                                       ccGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTCCAGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCa                                                                              >  1:346720/1‑150 (MQ=255)                                                                         gATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCa                                                                                                                        >  1:1444219/1‑106 (MQ=255)                                                                         gATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCa                                                                                                                        <  2:1444219/106‑1 (MQ=255)                                                                           tGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCC‑ACGTTTTCCACCAAGGGTGCCATTGGCAGCCTGACCTTGACGGCTGCGCGGGAACTGGCACGCTTCGGCATCCg                                                                          >  1:34660/1‑150 (MQ=255)                                                                            gagaGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCgt                                                                         <  1:829322/150‑1 (MQ=255)                                                                               agCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTT‑GGCATCCGTGtga                                                                     <  2:1254497/150‑1 (MQ=255)                                                                                 cggcgAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGtgat                                                                    <  2:1865221/150‑1 (MQ=255)                                                                                      aGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCa                                                               <  1:1263213/150‑1 (MQ=255)                                                                                       gCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCAt                                                              <  2:590594/150‑1 (MQ=255)                                                                                          tGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGc                                                           >  1:688477/1‑150 (MQ=255)                                                                                            gggTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTc                                                         <  2:351075/150‑1 (MQ=255)                                                                                               tcatcatcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCgg                                                      <  1:660848/150‑1 (MQ=255)                                                                                               tcatcatcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCgg                                                      >  2:451982/1‑150 (MQ=255)                                                                                                  tcatcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTa                                                   >  2:1701158/1‑150 (MQ=255)                                                                                                    atcgACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATc                                                 >  2:1459930/1‑150 (MQ=255)                                                                                                     tcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCgg                                                      >  1:822830/1‑144 (MQ=255)                                                                                                     tcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCgg                                                      <  2:822830/144‑1 (MQ=255)                                                                                                     tcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCt                                                >  2:245699/1‑150 (MQ=255)                                                                                                     tcaACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCGACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCt                                                <  1:591942/150‑1 (MQ=255)                                                                                                      caACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGc                                                                                                                          >  2:1127459/1‑75 (MQ=255)                                                                                                      caACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGc                                                                                                                          <  1:1127459/75‑1 (MQ=255)                                                                                                       aaCACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCttt                                              >  1:220820/1‑150 (MQ=255)                                                                                                        acacGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTg                                             >  1:1528506/1‑150 (MQ=255)                                                                                                          acGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGa                                            >  2:1406290/1‑149 (MQ=255)                                                                                                          acGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGa                                            <  1:1406290/149‑1 (MQ=255)                                                                                                            ggCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAAc                                         <  1:1701158/150‑1 (MQ=255)                                                                                                                      gccgccTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTatgatg                               <  2:1659599/150‑1 (MQ=255)                                                                                                                       ccgccTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATgg                              >  1:1000262/1‑150 (MQ=255)                                                                                                                                aCGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGa                     <  1:1912766/150‑1 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTgg                                                                                          <  1:432718/77‑1 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTgg                                                                                          >  2:432718/1‑77 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATcc                                                                           >  2:1698033/1‑92 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATcc                                                                           <  1:1698033/92‑1 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTcc                                                        >  2:1417627/1‑111 (MQ=255)                                                                                                                                    ccAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTcc                                                        <  1:1417627/111‑1 (MQ=255)                                                                                                                                     cAGATTGGCCAGGACGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGAt                <  2:877048/150‑1 (MQ=255)                                                                                                                                        aTTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAg             <  2:767590/150‑1 (MQ=255)                                                                                                                                        aTTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCCAATGGCCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAg             <  1:1277026/150‑1 (MQ=255)                                                                                                                                            gCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTAc         >  1:305936/1‑150 (MQ=255)                                                                                                                                            gCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTAc         >  1:1014836/1‑150 (MQ=255)                                                                                                                                              cAGGCTGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTACGt       <  1:1599606/150‑1 (MQ=255)                                                                                                                                              cAGGCTGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTACGt       <  1:1599269/150‑1 (MQ=255)                                                                                                                                                   cgccTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTGCTTGCTTc  >  1:1021691/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                       ATCAACGTCAACCTGATCGGCAGCTTCAACCTGTTGCGTCTGGCTGCGGCGGCCATGGCCGAAGGGGCTGCCGATGAGAGCGGCGAGCGTGGGGTCATCATCAACACGGCCTCCATTGCCGCCTATGACGGCCAGATTGGCCAGGCCGCCTACG‑CCGCCTCCAAGGGTGCCATTGCCAGCCTGACCTTGCCGGCCGCGCGCGAACTGGCACGCTTCGGCATCCGTGTGATGACCATCGCTCCGGGTATCTTTGAAACCCCTATGATGGCCGGCATGAGCGATGAGGTACGTGCTTC  >  NC_002947/2525351‑2525646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑