Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,797,544	A→T	7.6%	C165S (TGC→AGC)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,797,544	0	A	T	7.6%	184.6 / 6.4	66	C165S (TGC→AGC)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base T (3/2);  total (36/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 = 6.48e-01

CGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTG  >  NC_002947/3797396‑3797692                                                                                                                                                     |                                                                                                                                                     cGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGCCCCCGCCGTCGCCGATGCCTGGATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCGCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAg                                                                                                                                                     <  2:215856/150‑1 (MQ=255)  gCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAgg                                                                                                                                                    >  2:544627/1‑150 (MQ=255)        cgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGATCAGGCTGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCg                                                                                                                                              <  1:1553092/150‑1 (MQ=255)         gaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCgg                                                                                                                                             >  2:684204/1‑150 (MQ=255)           tGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGTTCCACGGTGTTGCCTCAGGCCAGCGGCa                                                                                                                                           <  2:218150/150‑1 (MQ=255)            ggcgggcgCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaa                                                                                                                                          <  1:771813/150‑2 (MQ=255)            ggCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCat                                                                                                                                          <  1:942711/150‑1 (MQ=255)                             gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTgg                                                                                                                         <  1:344117/150‑1 (MQ=255)                                    gCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTAt                                                                                                                  <  2:1443277/150‑1 (MQ=255)                                      tgatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCa                                                                                                                >  2:1863415/1‑150 (MQ=255)                                       gatgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAc                                                                                                               >  1:998062/1‑150 (MQ=255)                                           ccATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                           <  2:1238701/149‑1 (MQ=255)                                           acacTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                           <  1:1006833/150‑1 (MQ=255)                                           aCATTGACCACGCCGTCGCCGATGCCTCCATCGTGCCGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                           <  2:1842783/150‑1 (MQ=255)                                            cATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCgggg                                                                                                          >  2:1501833/1‑150 (MQ=255)                                            cATTGACCACGCCGTCGCCGATGCCGGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCgggg                                                                                                          <  2:606219/150‑1 (MQ=255)                                             aTTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGc                                                                                                         >  2:430401/1‑150 (MQ=255)                                             aTTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGc                                                                                                         >  2:420916/1‑150 (MQ=255)                                             aTTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGGGGCCAGTATCACCggggg                                                                                                         >  2:219010/1‑149 (MQ=255)                                              ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCg                                                                                                        >  1:596232/1‑150 (MQ=255)                                                 cccACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTc                                                                                                     <  2:1434910/149‑1 (MQ=255)                                                   cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCa                                                                                                   >  2:622024/1‑150 (MQ=255)                                                     cGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACg                                                                                                 <  1:1501833/150‑1 (MQ=255)                                                         gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCa                                                                                                                      <  2:237207/125‑1 (MQ=255)                                                         gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCa                                                                                                                      >  1:237207/1‑125 (MQ=255)                                                                     tGCNTCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAacac                                                                                 >  2:1201350/1‑150 (MQ=255)                                                                        aTCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAAcaccac                                                                              <  1:236295/150‑1 (MQ=255)                                                                               agcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAg                                                                       <  1:1640615/150‑1 (MQ=255)                                                                                gcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                      >  2:1621329/1‑150 (MQ=255)                                                                                gcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                      >  1:899559/1‑150 (MQ=255)                                                                                  aCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgc                                                                    <  1:684204/150‑1 (MQ=255)                                                                                  aCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCAATGCCCAACACCACGCCGCAGGgc                                                                    <  2:596232/150‑1 (MQ=255)                                                                                   ccTAGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcg                                                                   >  1:1073483/1‑150 (MQ=255)                                                                                       gCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGCTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgcg                                                               >  2:800511/1‑150 (MQ=255)                                                                                            tCAGCCGATGGACCGCCGCGCTCAGCTCGGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAggg                                                          <  1:1369509/150‑1 (MQ=255)                                                                                             cAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCATGGc                                                         <  2:659128/150‑1 (MQ=255)                                                                                                 cGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                     <  1:622024/150‑1 (MQ=255)                                                                                                      gACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGCTGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGaa                                                >  1:257291/1‑150 (MQ=255)                                                                                                      gACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGaa                                                >  2:423671/1‑150 (MQ=255)                                                                                                        ccgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGc                                              >  1:438622/1‑150 (MQ=255)                                                                                                         cgccgGGGTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCt                                             >  2:118781/1‑150 (MQ=255)                                                                                                         cgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCt                                             >  2:1415298/1‑150 (MQ=255)                                                                                                         cgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCt                                             >  2:1788903/1‑150 (MQ=255)                                                                                                          gccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTg                                            >  2:86737/1‑150 (MQ=255)                                                                                                           ccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCCGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGc                                           >  1:364257/1‑150 (MQ=255)                                                                                                                cTCAGCGCCGAGGCCTTGAGCACCACGGTGTGGCCGCAGGCCAGCGGCATGGAAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGGGCGATGCCCAACACCCCGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGa                                      <  2:343200/150‑1 (MQ=255)                                                                                                                   aGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTg                                                    <  2:1729139/133‑1 (MQ=255)                                                                                                                   aGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTg                                                    >  1:1729139/1‑133 (MQ=255)                                                                                                                     cTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTc                                 >  2:1885789/1‑150 (MQ=255)                                                                                                                        cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGATACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGa                              >  1:1788997/1‑150 (MQ=255)                                                                                                                         gAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                             >  2:1103405/1‑150 (MQ=255)                                                                                                                         gAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                             <  1:430401/150‑1 (MQ=255)                                                                                                                               ttGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATc                       >  1:1492298/1‑150 (MQ=255)                                                                                                                                 gAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCAc                     <  1:544627/150‑1 (MQ=255)                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGGTGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGGTCACTTCAcc               >  1:1180757/1‑150 (MQ=255)                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAAc                                    <  2:1442536/129‑1 (MQ=255)                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAAc                                    >  1:1442536/1‑129 (MQ=255)                                                                                                                                             catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                     <  1:1028954/104‑1 (MQ=255)                                                                                                                                             catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                     >  2:1028954/3‑106 (MQ=255)                                                                                                                                             catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGAt         >  1:458345/3‑150 (MQ=255)                                                                                                                                              ttGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc        <  2:708296/150‑1 (MQ=255)                                                                                                                                              atGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc        >  2:501411/2‑150 (MQ=255)                                                                                                                                                 ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGc     <  1:1243580/150‑1 (MQ=255)                                                                                                                                                 ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGc     >  1:1102483/1‑150 (MQ=255)                                                                                                                                                  cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCg    >  2:1700712/1‑150 (MQ=255)                                                                                                                                                    cAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTg  <  1:1621329/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTG  >  NC_002947/3797396‑3797692

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